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Обобщенный метод скаляризации динамических упругих полей 
в трансверсально-изотропных средах и его новые применения" 


И. П. Мирошниченко, В. П. Сизов?” 


' Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 


*Ростовский-на-Дону научно-исследовательский институт радиосвязи, г. Ростов-на-Дону, Российская Федерация 


Сепега| зса]ат17аНоп тео4 о? дупапис еа5@с Яе!4$ т гапзуегзаПу 150горе тед1а ап@ 165 пе\у 


аррЁсаНоп$ " 
Г.Р. МиозписвепКо', У. Р. Зо?" 


'Ров Зе Тесвиса! ОлиуегзИу, Возюу-оп-оп, Визз1ап Еедеганой 


2? Возюу З‹епийс ВезеагсЬ шзНние юг ВаФосоглитигсаноп, Возюу-оп-Ооп, Визчап Еедеганоп 


Введение. При исследовании тензорных упругих полей 
перемещений, напряжений и деформаций в слоистых 
конструкциях из различных материалов, включая 
трансверсально-изотропные композиты, успешно 
применяется эффективный метод скаляризации тен- 
зорных полей. Данные поля могут быть выражены 
через скалярные потенциалы, соответствующие ква- 
зипродольным, квазипоперечным и чисто поперечным 
волнам. Такая скаляризация возможна, если рассмат- 
риваемые объекты являются тензорами относительно 
подгруппы общих преобразований координат, когда 
локальный аффинный базис имеет один инвариантный 
вектор, который совпадает с осью материальной сим- 
метрии материала. При этом в известных работах рас- 
сматриваются конструкции, где этот вектор совпадает 
с нормалью к границе между слоями. Однако, для 
практики представляют интерес и другие случаи вза- 
имного расположения оси материальной симметрии 
материала и границы между слоями. 

Целью является дальнейшее развитие применения 
метода скаляризации в граничных задачах динамиче- 
ской теории упругости на случаи произвольного рас- 
положения оси материальной симметрии по отноше- 
нию к границе между слоями. 

Методы 
методический аппарат разработан на основе использо- 


исследования. Предлагаемый — научно- 
вания обобщенного метода скаляризации динамических 
упругих полей перемещений, напряжений и деформаций в 
трансверсально-изотропных средах. 

Результаты исследования. Получены новые расчетные 
соотношения для определения полей перемещений, 
напряжений и деформаций в трансверсально-изотропных 


средах на случаи произвольного расположения осей 


Гигодиспоп. Ап е йаеп цесЬмаие оЁ {епзог Не! зсаапта- 
@оп 15$ зассеззШу зе У\@Ие шуезИгайпо {епзог @азс 
Неа оё а15р]асетеп, зте5зез ап 4еогпайоп$ ш фе Пау- 
еге4 знас@агез оЁ 41егепЕ таепа]з, шсе гапзуегзаПу 
1з0тор1с сотрозйез. ТБезе Йе!4$ сап Бе ехргеззе4 Шгопи?В 
Фе зса]аг роепНа!$ сотгезропат? ю Фе апаз1-1оп2Надта!, 
Чиаз1-фгапзуегзе, ап ‘гапзуегзе-опу \’ауез. бис зса]апта- 
боп 1$ розз1е 1 Ше об]ес{ ипаег сопз1дегайоп аге 1епзот$ 
те!айпо №ю Фе заботойпр оЁ репега| соог@ а сопуег$1опз, 
УПТеп Ше 1оса| аЁЙпе Ъа$1$ Баз опе шуапап{ уес{ог фай сош- 
с14ез \Иф Фе таепа| зуплтеу ах1$ оЁ е таета1. АЕ 15, 
фе Кпо\уп рарег$ сопз1Аег зтасвагез УПеге 1115$ уесюг сош- 
с14ез \ИН бе поггла! ю Фе Боипдагу Бебмееп 1ауегз. Но\уеу- 
ег, олег сазез оР Фе шиша! атапоетерё оР фе тшаепа] 
зуттейу ах15 оЁР фе шаепа| ап Фе Боипдамез Бебмееп 
1ауегз аге оф ииегез{ оп Фе ргасНса1 э14е. 

Маепта[ апа Мефо4. Тве \могк об]еспуе 1$ Каийег 4еуе]- 
ортеп{ оЁ Фе зса!аптайоп тефо4 аррИсаноп ш фе Боипда- 
гу уаае ргоепз оР 1е аупапис е!азисйу ШФеогу Юг Фе 
сазез оЁ ап агойгагу аггапоеглеп( оЁ фе шаема| зуштейу 
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материальной симметрии материалов слоев по отно- 
шению к границам между слоями. 

Обсуждение и заключения. Предлагаемый научно- 
методический аппарат успешно использован при 
определении напряженно-деформированного состоя- 
ния в слоистых конструкциях, выполненных из транс- 
версально-изотропных материалов, и при анализе ре- 
зультатов диагностики состояния плоскослоистых и 
слоистых цилиндрических конструкций, находящихся 
в эксплуатации. 
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Введение. При исследовании тензорных упругих полей перемещений, напряжений и деформаций 
в слоистых конструкциях из различных материалов, включая трансверсально-изотропные композиты 
[1-3], успешно применяется эффективный метод скаляризации тензорных полей, предложенный в [4], 
когда отмеченные поля могут быть выражены через скалярные потенциалы, соответствующие квазипро- 
дольным, квазипоперечным и чисто поперечным волнам соответственно. 

Такая скаляризация возможна, если рассматриваемые объекты являются тензорами относительно 
подгруппы общих преобразований координат, когда локальный аффинный базис имеет один инвариант- 
ный вектор, который совпадает с осью материальной симметрии материала. В [5] рассматриваются кон- 
струкции, где этот вектор совпадает с нормалью к границе между слоями. Однако для практики пред- 
ставляют интерес и другие случаи взаимного расположения оси материальной симметрии материала и 
границы между слоями. 

Цель исследования. Целью работы является дальнейшее развитие применения метода скаляриза- 
ции в граничных задачах динамической теории упругости на случаи произвольного расположения оси 
материальной симметрии по отношению к границе между слоями. 

Постановка задачи. Первоначально в системе координат с допустимым репером находятся ком- 
поненты тензоров перемещений и напряжений относительно этого репера. Затем, зная все компоненты 
тензорных полей в данной системе координат, можно путем перехода к системе координат, связанной с 
границей между слоями найти нормальные и касательные компоненты перемещений и напряжений, со- 
ответствующие площадке, лежащей на границе. Эти компоненты далее используются для удовлетворе- 
ния граничным условиям задачи [5]. 

Более подробно изложим вышесказанное на примере решения задач для плоскослоистых кон- 
струкций. 

Рассмотрим случай, когда главная ось симметрии трансверсально-изотропного материала со- 
ставляет угол © по отношению к нормали поверхности плоской границы между слоями (рис. 1). При 


этом считаем, что поле не зависит от координаты у. 





ЕК 


Рис. 1. Расчетная схема 
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Расчетные соотношения 
В соответствии с [4]: 


17) 


(р) Г 1 
[ОА - [5% +572, У ‚9+5 уу +857 ж+ 8 (6—6); 
8 


(2) 
с, б.н +4 ‚5767 +4, бум, и, У,У,)ф+ (1) 
(Г) (Т) 
+(а, 2,У, +, ‚575иу, + СА ‚+4. У У У )и+ 


+2 ва, (У 55" УМУ + аи (616 515 УК) у. 


Компоненты тензорных полей перемещений И, и напряжений 5, в системе координат, связан- 





ной с анизотропией материала хх“ › могут быть записаны следующим образом: 
(Г) 


В (В (Г) г. (т) 
= &(Р+Ь, + (= В, ш+а)и; 
8 
(2) (В) (Т) й в 
к = ВР, ф- Ве (Т) ; 
2 (2) 
(р р В (г) т 0 0 ый 
о =[а+4,-(4:+4.)5 ф-+т5(а+а.+а:-а.5`)»; 
1 (2) (22 000 (Г) ] (0 (т) 
ео 4; 6 В+а, рев 4-4, ); 
(2) (7) 4 (т) (т) (Г) (Г) 


2 . 2 
бк =(Ч+а.В ф-т 5 (4-4. В)». 
Здесь потенциалы квазипродольных ф и квазипоперечных у волн должны удовлетворять волновому 


(2) (Г) 
уравнению с соответствующим волновым числом & ( & или & ) и имеют вид: 

















ф= фе Ех” е’ В ей - 
(Г) й (Г) к (3) 
ж= ие Ех е! В хк , 
ф 
гдеё +В’ ==”; (д+8”) =0. (4) 
и 
Входящие в (2) коэффициенты Д и 4, определены в [5] и имеют вид: 
(2) 
(р) - 
8 (С, +Си) | 
ее о) ВОИ , 
в’р-й” Си (2° (С С, -Си) 
и. 9 
р _& фр-й (2С„+С,)-(=-Й С, __& р 
о В 
й юр-й Си- (а -Й (С, -С,-Си) 8-й 
(Г) 
Е (т) ео бы, 
В в = (т) пи} (5) 
Е в’р-й" (2Си +Сз)-(2"-®)С 
М>) (2) (2) (2) (2 (В (2) (Р) (2) ([) (2) (2) 
2 о а, === а, -й(2а, О, +а, О. +а. 0); 
В (2) (> (р (1) (2 
о а, бар. +2а.(20-+2,); 4, =2а, О; 
= (т) 
а. (Г) т 0 (г) 2 т) 
5= а, =а, &(-Б)+а; 8-тВ ; 


8 
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(Г) 





(Г) г) (1) (Г) г) р? (т) 
4, =а, 81-0.) +2а, а+а;| &-т р, |; 
5 
(т) 
(Г) (Г) (г (72| (т) ] 
4; =2а, &+2а.| #2. у |; 44 =2а, В, д; (6) 
8 8 
. 1 . 1 . 
а, =Си-2а,; а, =5(С,-С.); а =5 (С, +С:)+Сь +Сы 44; 
1 
а. = ый -С›)+Си;а; =С,+С,-2(С,+2Си), (7) 


где С, — модули упругости материала, записанные по свернутому индексу [6-7]. 


В формулах (5) и (6) для этих коэффициентов: 
_  с0$0 : _ 510 . _ 8. 
— соз(о+ 0) 7 — соз(о+ 0) РЕ ый с0$( +0) ` 
В граничные условия динамических задач теории упругости входят компоненты перемещений 


О., 0, и напряжений б., б.., записанные в системе координат 2 , х ‚, связанной с границей (рис. 1). 


2х? 


Л К > 
‚ х` связаны между собой соотношениями [8]: 


Системы координат 2, х их 
Е=-х“ шах” соза: Хх =х^ со ах” зта (8) 
или: 

х^ =Хс0$“—2510;х’ =Х 510, +2 ©0850. (9) 

Используя соотношения (8) и (9), по формулам [5]: 

т т 
0, “а. —. А. - 
ох' ИВ: Г: Вы 


где коэффициенты с чертой относятся к системе координат 2 ‚ х (назовем ИХ «НОВЫМИ» координатами), 


(10) 


а без черты к системе координат х”, х^ («старые» координаты). 


Запишем компоненты полей перемещений (0., И, и напряжений б., б.,, входящих в граничные 


2х? 
условия через компоненты (2): 

О. = сова О, - зах; 

О, =эта О, +соз ох; (11) 


б.. = 605’ мо, -зт2ао к +51” акк; 
в 1. . 
б. =512а(с, -ск)+@- 29" 9). (12) 


В «новых» координатах потенциальные функции (3) имеют вид: 


ми (о (р 
и: Хзша-+2 со 9) 7 Хс0$9.-2 эта.) . 
ф = фе 5 ( де В( ) 


2 


(Г) (Г) 
1Е (уяпа+7 08а) 1 В (хс050.—2 зш 0) 


®= ие 


2 


а перемещения (11) запишутся следующим образом: 


Т Г. 
> КОЛО (Г) = т 0 ф- (Г) Е (т) _ 
И. =созо[ 6 (2+20,)ф+(-ВБ, о 2)" -зшто[ 6 О, ф+(-2, Ре 2)" (13) 
8 8 
(Г) (Г)(Т) 
Е с Е _ Фо ФЕВ _ 
И, =зто[КО+ ,)$+(-Б; + 8)" +созой В В, ф-Б, ми] (14) 
8 8 


Здесь волновые числа &, В, © определены относительно «старого» репера (рис. 2). 
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Рис. 2. Схема «взаимосвязи» волновых чисел 


Волновые числа &, В, которые являются проекциями вектора 2 на «старый» репер х”, х^ и 


волновые числа ©, в. которые являются проекциями вектора 2 на «новый репер» 7. ‚ х связаны меж- 


ду собой следующими соотношениями: 








_  с0$0 _  с0$0 с ор_ 900 _ 510 . 
зш(а +0) ы с0$(0-+ 0) -_ зт(о + 0) г с05(© + 0) 2 





бы +0) соз(а+0)’ 


Подставляя (15) в (14), получим выражения для перемещений в «новых» координатах: 





и (2 (о т (Г) 
0. 18. (+0965 р бе (-соз0 р, — 909% я 
05(% + 0) с0$(о + 9) 
_ (0 с0805то 2 @  зтдсоза 2 9 зш(а+0)соз@ Ша 
9, =1а,[- о +Ь,) + РБ +а , (16) 
со$(а, + 9) с0$(© + 9) с05(о + 9) с05(© + 0) 
где, в соответствии с (13) и (15): 
.ы к (г (Г) 
ф=фе*“е*; т= ме“ "ео". (17) 


Аналогично, используя (2), (14) и (15), из (12) получим соотношения для компонент напряжений, 


входящих в граничные условия задач: 























шим ШО 20 102 @) со50зт0 2 ., Ф Ф зшо 
б_ ={с05? ма+а,-(а+а,)— ^^ —е?]+зт2а(-а ча, е?+ т? о(а- а, ^^ 92 ф+ 
== 0—0 “за +6у** 1 © й оз(а+6** (а, “оа+в) "$ 
воза с050 ел е соз° 0 2) | 510 г 0 (7 со5’0 ел 
с0$(© + 0) Е * со5? (+0) 5 с0$(0+ 0) 5 “сока +в)” 
с050 ФФ  со’0 ® _ 
т 20, 1е.(а-а. й 82); 
со$(а + 9) с05`(@+0) 
1. О 0 Ф  со520 2 @ 50 © . 1 © соз0зт@ 9 — 
5. ={ эш2а[а,-—(а,-+а,)——— в? 21-а-2 51? о) (-а,-+а 21 ф+ 
2х — [4,- (а, “сова +6)* “сои +0) | ( х> 3 “со а+6) "9 
(т) (т) (т) 29-х? 0 : (т) (Г) 2 (Г) 
Е а Пт (18) 
2 со$(а, + 0) с05° (&+0) ^ с0$(+0) 2 с0$`(&+0) 


Таким образом, компоненты перемещений и напряжений, входящие в граничные условия (систе- 
ма координат 2, х связана с границей), определяются соотношениями (16)—(18), где проекции волново- 


го вектора & в системах координат 2, х их”, х^ связаны соотношениями (15). 


Используемые в [5] матрицы С, характеризующие волновые свойства слоев [9-10], и которые 
являются основными элементами при построении конкретных решений граничных задач в слоистых 


В р://уезииК.Чопзеа.ги 


конструкциях в данном случае, когда оси 2. ‚ ея образуют с компонентами репера г“, е^ угол @, имеют 


вид: 
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с- бо С) (19) 
С ба 
где 
(р 
и ме р 2" 0 оо 
и 12 п 12 
Са = Со, = а оо 
с- С С Си С 0 еп“ еве*, 
ыы (Е) 2 
С: Со Си С 18:7 
а ес Си = = с —_ | _ 
41 42 41 42 (т 
оо 0 ев 
(г) (р 
ВВ. (20) 
2 2) соз@созо (2 
Си =18,(Р+ 2,); 
со$(а, + 9) 
(г) 0 соза 
Сь=8, (-с080 р); 
со$(а, + 9) 
‚2 созбзто, @ зш0соза (0 
И уе, 
с0$(© + 9) со$(а, + 9) 
о | $1(о-+ 0) соз0 р, та, 
с05(© + 09) с0$(о + 9) 
0 Ш 9 со5? 0 (2) , 1 (2) (2 соз0зт0 (Е) , (2) (2) $1? 0 (2) 
С; = 05? [4 +а.-(а.+а,)—— 2] +эт 20(—4,+а С аят” (а 8? 
>. а, (аа сова) ® 2” “ озз(а+6)** ' 4 совы +6) 7 ) 
с050 ФФ ФФ Ф 9 со’ 510 ФТФ @ соз’0 
С; = с05' © га. (аа, +4, а. 2 &:) 18, 3 @4 2 &:)+ 
с0$(© + 9) с0$°(м-+0) ^ с05(о + 9) 2 со$`(@+0) 
ВЕ со50 ФФ Ф зш’0 
+5 9 га, (4-4. 2 в 
с05(0- 0) с05` (4+0) 
1 Фо ® соз?0 © ® зи?0 © . 1) < соз0зт0 2 
С, =-зт2а[а,-(а.+а - 21-(1-2 52? а)(-а.-+а ра 
о, 4, (а, оо +0)* оао) ) 9° “сов (а +6) =" 
1. 0 0) (7) 1) соз? 0—эт” 0 | эт 0 19 @ с05?0 
Сома ааа И 29) +а-2? а) аа =). (20 
2 с05(0-+ 9 с05` (9-0) со$(а-+0) 2 с05`(&-+0) ^ 


В случае, когда а«=0 элементы (21) матрицы С совпадают с соответствующими элементами 
этой матрицы для случая совпадения нормали к границе и направления оси симметрии материала, кото- 


д | 
рые представлены формулами (3.82) в [5]. Если © = 5 то получается случай, когда ось симметрии материала 


касательна к поверхности границы и формулы (21) совпадают с выражениями (3.89)—(3.91) в [5]. 

Зная выражения для матриц С, можно построить решения различных задач, используя научно- 
методический аппарат, описанный в [5]. 

Выводы. Предлагаемый научно-методический аппарат успешно использован при определении 
напряженно-деформированного состояния в слоистых конструкциях, выполненных из трансверсально- 
изотропных материалов, и при анализе результатов диагностики состояния плоскослоистых и слоистых 
цилиндрических конструкций, находящихся в эксплуатации. 
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Контактная задача для двухслойного цилиндра" 
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Введение. Актуальность исследования контактных задач 
для цилиндрических тел обусловлена необходимостью 
проведения инженерных расчетов на контактную проч- 
ность валов, стержней и трубопроводов. В настоящей ра- 
боте изучается новая контактная задача статической тео- 
рии упругости о взаимодействии жесткого бандажа с бес- 
конечным двухслойным цилиндром, состоящим из внут- 
реннего сплошного и внешнего полого цилиндров, между 
которыми выполняются условия гладкого контакта. 
Наружный цилиндрический бандаж посажен с натягом и 
имеет конечную длину. При помощи интегрального пре- 
образования Фурье задача сводится к интегральному 
уравнению относительно неизвестного контактного дав- 
ления. 

Материалы и методы. Рассматриваются разные комбина- 
ции линейно-упругих материалов составного цилиндра. 
Исследуется асимптотика функции-символа ядра инте- 
грального уравнения в нуле и бесконечности, играющая 
важную роль для использования аналитических методов 
решения. Для решения интегрального уравнения вводится 
основной безразмерный геометрический параметр и при- 
меняется сингулярный асимптотический метод. 
Результаты исследования. В соответствии со свойствами 
функции-символа предложена специальная легко факто- 
ризуемая аппроксимация этой функции, пригодная в ши- 
роком диапазоне изменения параметров задачи. При по- 
мощи метода Монте-Карло рассчитаны параметры этой 
аппроксимации. Получены асимптотические формулы как 
для контактных давлений, так и для их интегральной ха- 
рактеристики. Расчеты сделаны для разных материалов и 
относительных толщин цилиндрического слоя, в том чис- 
ле для тонкостенных слоев. 

Обсуждение и заключения. Полученные асимптотические 
решения эффективны для относительно широких банда- 
жей, когда размер области контакта превышает диаметр 
составного цилиндра. Важно, что используемый метод 
остается применимым и для случаев, когда внешний ци- 


* Работа выполнена по гранту РФФИ 18-01-00017. 
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линдрический слой можно рассматривать как цилиндри- 
ческую оболочку. Асимптотические решения можно ре- 
комендовать инженерам для анализа контактной прочно- 
сти упругих деталей цилиндрической формы с упругим 
покрытием из другого материала. 


Ключевые слова: теория упругости, контактные задачи, 
составной цилиндр, аппроксимация, асимптотика. 
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Введение. Динамическая контактная задача для преднапряженного упругого цилиндра, наполненного 


жидкостью, изучалась в работе [1]. Статические контактные задачи для однородных упругих тел цилиндриче- 
ской формы рассматривались в работах [2—6] при помощи регулярного и сингулярного асимптотических мето- 
дов. Было установлено [4], что для цилиндрических тел символы ядер интегральных уравнений контактных 
задач характеризуются более сложным асимптотическим поведением в нуле и бесконечности, чем, например, в 
контактных задачах для упругой полосы. Это потребовало применения усложненных аппроксимаций этих сим- 
волов легко факторизуемыми функциями при использовании сингулярного асимптотического метода. Особен- 
но усложняется аппроксимация для полых тонкостенных цилиндров [6]. Предложенная аппроксимация [6] эф- 
фективна даже для случаев, когда тонкостенный упругий цилиндр можно рассматривать как цилиндрическую 
оболочку [7]. Исследовалась контактная задача о взаимодействии упругого кольца с упругим цилиндром [8]. 
Износ упругого цилиндра анализировался в работе [9]. Цель настоящего исследования — получить решение 
контактной задачи для составного двухслойного упругого цилиндра на основе сингулярного асимптотического 


метода и эффективной аппроксимации символа ядра интегрального уравнения. 

Материалы и методы. В цилиндрических координатах и, 2 (при осевой симметрии) рассмотрим беско- 
нечный упругий составной цилиндр внешнего радиуса А, который состоит из внутреннего сплошного цилиндра 
радиуса А!<А с упругими параметрами у1, С, (коэффициент Пуассона и модуль сдвига) и внешнего цилиндри- 
ческого слоя с упругими параметрами у, С. Между слоем и внутренним цилиндром выполняются условия сколь- 
зящей заделки. Рассмотрим контактную задачу о взаимодействии описанного составного цилиндра с жестким бан- 
дажом по области |2<а. При заданном натяге бандажа 6 требуется определить контактные давления 

с, =-9(2) (25а). 

Используя интегральное преобразование Фурье для решения собственно смешанной (контактной) краевой 

задачи для уравнений Ламе упругого равновесия и вводя безразмерные обозначения (штрихи далее опускаем) 








Ат, 5' -5, Н -=, ик, = к я (1) 
получим следующее интегральное уравнение относительно 4(б): 
ДЕ =лб (|(<1), (= [дшдсовбидаи, (2) 
где символ ядра имеет вид ” | 
п --Ч чак, (8) 
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а= (455 А ‘ Аз Ан А 4, = А, А») (АА — А.А) = 
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= АА ХА Ав = АА ›1, 
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—А,5.455 [А (Аз Ава В Аз Аз ) =. Ар (АзАв = А Ав У, 


4, =—(А, Аи - А.А, А (АА, -А,Ак)- ААА] + 
+А, А [А (Аз А, -А,А:)-А (А.А, - АА], 

А =и\ -2(-у)Л, А, ==—иК. -2(-у)К, А, =и|, А. =—иК, 
А, = (3—2у)ш — (и? +4 -у))1, А, =—(3-2\иК, — (м? +4А-у))К,, 
Аз =и(-и), Ад =—и(К! +иК,), 
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А =иКГ —=21-у)", А’ =-иКК, -2А-У)К'", Аз =ий, Ад 
А; =и К =-2-\,), Ак =иш, 
А, = 2 и -Е+4А-УЕГ, 
А, =—-—2иК: — (и +4а-УЕУКУ, 
Аз =ик Ч и? , Ад =-иКИ К! и? КЕ, 
А, = 6—2)" —Рк+4-У Е, А = ви -шЬ), 
=Г(и), К, =К, (и), Г =Г (К), К” =К, (иК), п=0, 1. 


Здесь //(и), К‚(и) — модифицированные функции Бесселя. Безразмерный параметр ^, характеризует от- 
носительную ширину области контакта. 
Функция Г(и) в нуле и бесконечности ведет себя следующим образом: 


(У-1@+=-\, +5%,) +”, 
2(у-1[+10+=-у, +=%,)+А"=] | 
& = (У+0(\, -9-=(-—П(\, +1. (4) 





Низ (м) = 240) = 


гр 
Ци) = —+—+0(и”) (и—> +0), Р=1-2у. 
ии 


При А=0 значение Г(0) совпадает с известным для однородного сплошного цилиндра [4]. 

Для решения уравнения (2) применим сингулярный асимптотический метод [3,4], эффективный при 
достаточно малых значениях /.. С целью применения метода Винера-Хопфа [10] использовалась легко факто- 
ризуемая аппроксимация функции Г(и) (3) выражением 


ми? + В? Бр и? + А*С? В ур) 
: 5_ехр| == ‚ @ = ехр| — |. 
и +С ) [2 +10“ 10 
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При расчетах брали два случая: 1) железо внутри, цинк снаружи (&=2,126, у=0,27, у,=0,28); 2) 
алюминий внутри, цинк снаружи (&=0,779, у=0,27, у,=0,34). В таблице 1 даны значения параметров 
аппроксимации (5), ее относительная погрешность на действительной оси 09 (%), рассчитанные при 
использовании метода Монте-Карло при разных относительных толщинах внешнего слоя К. 

Таблица 1 
Параметры аппроксимации 





А В С 0 А В С 0 


К Железо внутри цинка Алюминий внутри цинка 


0,10 1,230 1,549 4,638 2,5 0,451 6,113 3,394 
0,30 1,291 1,563 4,042 2,5 0,888 5,977 3,029 
0,50 1,296 1,851 5,006 3,0 2,399 1,467 2,524 
0,70 1,258 2,057 6,376 3,0 0,720 7,945 2,512 
0,90 1,024 4,360 92,381 5,0 1,481 9,890 9,369 
0,99 3,540 24,178 5,347 7,0 1,178 5,535 | 235,946 10 
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Поскольку, как следует из (4), 
1-у%, 
э 

21 -—у)(1+\,) 
то в случае малых # аппроксимация, рассчитываемая по формуле (5), для тонких внешних слоев должна услож- 
няться путем увеличения числа входящих в нее параметров. 

Результаты исследования. В результате применения метода Винера-Хопфа главный член 
асимптотического решения интегрального уравнения (2) при малых /^. можно построить в форме 


_6| [1+6 Е 1 
4() т т |4 >. т ((х|51), (7) 


_ #5) + 5) 


И’ ($) = аи, р ег (/В5 +[5-1 О(АС,5) , 


О(Е,5) = Е ехр(-Ез)ег /(В — Е)5) + т ег (Ву). 


Здесь егКх) — интеграл вероятностей. 
Для интегральной характеристики решения 


Р= [ час (8) 


вит (0) = (6) 





Ця) = ори Эк, 0?Эат, 
ИЯ 





на основании формул (7) получим выражение 


Р 2 2 2 2 р: ь 
С 2(2)+= (2 о О эк. (0? эат, (9) 
2 т =) а ехр( вы Е И Сиб, 


Т(Е,5) = < тм =) а ехр(-— во - 
С \ехр(-Р 5) — 
( = ЕЕ еРСВ-Е)5). 


Как показывают расчеты, погрешность асимптотик (7), (9) при ^.<1 не превышает (5+0)%, где 0 — по- 
грешность аппроксимации (5). 
В таблице 2 приведены значения интегральной характеристики Рб’', рассчитанные по формулам (9) 

















при разных значениях Ки ^. 
Таблица 2 


Значения Рб" 


























= 2 1 0,5 0,25 2 1 0,5 0,25 
К Железо внутри цинка Алюминий внутри цинка 
0,10 3,25 5,77 190,9 21,1 3,22 5,71 10,8 20,9 
0,30 3,37 6,02 11,4 22,2 3,20 5,65 19,7 20,7 
0,50 3,67 6,65 12,7 24,8 3,09 5,51 10,5 20,4 
0,70 4,22 7,78 15,0 29,4 2,96 5,33 19,2 19,9 
0,90 5,60 10,3 19,6 38,2 2,93 5,28 9,99 19,4 
0,99 6,91 12,3 23,2 44,9 2,99 5,31 9,94 19,2 



































Заключение. Как ВИДНО ИЗ таблицы яя с уменьшением А, интегральная характеристика контактных 
давлений возрастает, что связано с увеличением площади области контакта. Для случая более прочного 
материала внутри цинка (железа) контактные давления больше, чем для алюминия внутри цинка. С утончением 
слоя цинка вокруг железа (при возрастании 19) контактные давления существенно возрастают. При утончении 
слоя цинка вокруг алюминия этого не наблюдается, поскольку модуль продольной упругости (а также модуль 


Пожарский Д. А. и др. Контактная задача для двухслойного цилиндра 


Розпагзки Ш. А. апа ше ошйегз. Сотасй рто Мет рог а &то-1ауегей суйтаег 





сдвига) у алюминия немного меньше, чем у цинка. Найденные асимптотики можно рекомендовать инженерам 
для анализа контактных прочностных характеристик деталей цилиндрической формы с покрытием. 
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Сшаз Квагтапда ', 1. В. Ашураз””* 


люд ОшуегзНу, Гллаа, Зууедеп 
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Надежность оптимизации дизайна с использованием оптимальных факторов безопасности для 


крупномасштабных задач” 


Харманда М. Г.', Антибас И. Р.?"" 
' Лундский университет, г. Лунд, Швеция 


? Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 


Введение. Модель, основанная на оптимизации надёжно- 
сти, (КВОО) уменьшает структурный вес в некритических 
регионах, обеспечивает не только улучшенную конструк- 
цию, но и более высокий уровень уверенности в дизайне. 
Материалы и методы. Классический подход КВРО мо- 
жет быть выполнен в двух отдельных пространствах: фи- 
зическом пространстве и нормированном пространстве. 
Поскольку в вышеупомянутых двух пространствах требу- 
ется очень много повторных исследований, расчётное 
время для такой оптимизации является большой пробле- 
мой. Эффективный метод, называемый Орйтит ба у 
Еасюг (ОЗЕ), разработан и успешно применяется к не- 
скольким инженерным приложениям. 

Результаты исследования. Численное приложение по 
крупномасштабной задаче при усталостной загрузке пока- 
зывает эффективность разработанного метода КВОО от- 
носительно детерминированной оптимизации дизайна 
(ро). Эффективность метода ОЗЁЕ также распространя- 
ется на несколько режимов отказоустойчивости для 
управления несколькими выходными параметрами, таки- 
ми как структурный объем и атрибут повреждения. 
Обсуждение и заключения. Упрощенная стратегия внед- 
рения структуры ОЗЕ состоит из единственной задачи по 
оптимизации оценки проектной точки и прямой оценки 
оптимального решения с учетом составов ОЗЕР. Он предо- 
ставляет разработчикам эффективные решения, которые 
должны быть экономичными, удовлетворяющими требуе- 
мому уровню надежности с сокращённым расчётным вре- 
менем. 


Ключевые слова: оптимизация на основе надежности, 
структурная надежность, факторы безопасности 


* Работа выполнена в рамках инициативной НИР. 
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1. шаодисйоп 

\УВеп Беегииизис Оез1еп ОрНпиханоп (ООО) те#о4$ аге зе4, деегтиизНс орИтит 4ез1еп$ аге изиаПу 
ризНеа №0 фе 4ез1еп сопзгате Боипаагу, 1еауште Ш@е ог по гоот Юг {юегапсез (ог ипсецашиез) ш 4ез1еп, тапи#ас- 
фиге, ап орегайпе ргосеззез. бо, дейегиии$Нс орИтит 4ез1еп$ оМаштей \уИпойе сопз1АегаНоп оЁ ипсецатйез тау 
1еа4 1ю ипгепае 4ез1ет$, фегеоге сапе ог КейабИиу-Вазе4 Оез1еп ОрйпихаНоп (ВВОО). Ап КВРО зошйоп 
а гедисе$ Фе згасвага[ мел Е шт ипстса| ге21юп$ 40е$ поё оШу ргоу!е ап пиргоуе4 4ез1еп Бой а[50 а ШеВег 1еуё1 
о# сопй4епсе ш Фе 4ез1ет. ТВе Базтс 14еа 15 №0 соире Ше тепабИиу апа1у$1$ уф орнпихайоп ргоепт$. Т№1$ соч- 
рИп 15 а сошр]ех {азК, ага\уте оп а ШБ сотрийпе Ите ап4 сопуегоепсе заб Шиу, \ысЬ зепоч$ у Шшти5 Из аррИса- 
Шу ш геа] ргоегаз. То оуегсоте езе ЧИсиез, зеуега] пе о4$ Вауе Бееп е!абогайе4 [1]. ТВезе те#о4$ сап Бе 
с1аз;е по мо саегопез: патетса|! ап4 зеп-питепса| тео4$. Ап еЁЯсепе питепса]| тео са|Шеа Нубиа 
Мео4 1$ Базе оп зиияНапеой$ зоаНоп ю Ше гепабИИу апа Фе орнпитаноп рго Мет. 1 Ваз зассеззРаПу гедисе4 
Фе сотршаНопа| Яте ргоЫет. ТЬе а4уащасе оЁ Те Вубм4 те#фо4 аПо\уз и$ {0 зайзКу а тедите4 гепабшу 1еуе] Юг 
Ч1егеп{ сазез (айс, Чупапис, ...), Би е уесюг оЁ уанпа ез Неге сотатз Бой деегии1$Нс ап гап4от уама ез. 
То оуегсоше Бой А!Яси@ез, ап еЁИслеп{ зеп!-питенса! тефо4 са!еа Орйтит За у Еасюг (ОЗЕ) тефо4 Баз 
Бееп ргорозе4 ю зо1уе ргоетз ш $айсз [2], ап4 а!50 ап еЁЯслепЕ айегпануе зеп!-питенса! тео4 саПе4 За 
РошЕ Мео4 (5Р) Ваз Бееп ргорозе4 №0 зо!уе ргоМет ш 4упапис$ [3]. Кесепйу, а Коби5Е Нубма Меод (ВНМ) 1$ 
а150 Чеуеоре4 ю оуегсоте фе НМ ЧИ сшвез ш огаег №0 зо!уе ти@ахла| Райече датазе апа!уз15 ргоетаз [4]. Тве 
ВНМ 1еа4$ 1 гоби$( зо[айоп сотрагие мБ Ве НМ, Би Фе сотрийпе Яте 15 5НП а Ы2 дгауЪаск. 


2. Вепа ИИу-Вазе4 Оез оп Орйпи7аНоп 

2.1 Эеу@ортепт 6 

Тве сотриайопа! со оЁ зедиепйа1 КВРО арргоасВе$ 15 шисВ В1оКег {вап е ООО ргоседиге. Зеуега| 4еуе]- 
ортепё5 ассе]егайе4 1е изе оР Фе ВВРО тоде!. Те Вейаб у шдех АрргоасЬ (ВТА) апа фе Реогтапсе Меазиге 
АрргоасН (РМА) Бауе Бееп ргорозе4 [5]. Мехь 1е зедиепна! орипитайоп ап гепабИИу аззеззтепте (ЗОВА) 1$ 4е- 
уеюре4 №ю ппргоуе е еЁЯслепсу оЁ ргоБа Шзис орнпиханоп [6]. ТВе ЗОКА шефо4 етрюуз а зшве-оор зиаезу 
УВ а зепа| оЁ сус|ез оЁ деегил15Яс орипихайоп ап4 гепабИиу аззеззтепе. ТВе та]ог АЁНсиЙу Пез ш е еуача- 
Чоп оЁ Фе ргоба Ис сопзгаии в, \ушсв 15 ргобауе[у ехрепз1уе ап еуеп Фуегоез Юг тапу арр|сайопз. 1 15 с1еаг 
аЕ еЕРог$ \уеге Чтесе4 1ю\уага$ Фе аеуеортепЕ оЁ еРНстепЕ {есби1аие$ ю реюпт йе геПпабиу апа[у$15. Неге, Ве 
тепаб у шаех 15 сотрше4 Иегайуе]у Фа{ [еа4$ 10 ап епоглои$ аточиЕ о? сотршег ите ш Фе бое 4ез1еп ргосез$. 


2.2 Ваяс ВВОО югишаЧо0п$ 


ТгадопаПу, Юг фе гепабПиу-Базе4 4ез1еп орйпитайоп ргоседиге, мо зрасез аге изе4: фе рБуз1тса| зрасе 
ап4 фе погта|7е4 зрасе [7,8]. ТВегеРоге, пе гепабИиу-базе4 4е$1еп орниитаНоп 15 реогте4 Бу пезипе фе ЮПо\- 
шо (№0 ргоЫ ев: 


1. Орипихайоп ргоет: 


шт Х(х) 
зиБдесею &,(х)<0, А=Ь,...К (1) 
ап4 В(х, и) >В, 


УПеге Хх) 1$ Фе обуеснуе РапсНоп, ©(х) < 0 аге Фе аззос1айе сопзаш, В(х,а) 15 Ше гепаБИИу ш4ех оЁ Фе $гис- 
{аге, апа Д, 1$ Фе фагоей гепаб ИУ. 


2. Вепабщу апа[уз15: бе тепабИиу шаех В(х,а) 15 Фе пшипит 41$апсе Бебуееп 1е шпи зе РапсНоп Н(и) ап4 Фе 
опрш, ее Е1еиге 15. Ты шдех 1$ дейегилтеа Бу зо[уш фе шшитханоп рго ет: 
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ти а(и) в 
зиБесею Н(и)=0 


уВеге (и) 15 фе Ч $апсе ш Фе погта|Нте4 гап4отл зрасе, э1уеп Бу 4 = Уи ‚ апа Н(а) 15 фе регоплапсе Рапсйоп 


(ог Ший з(ае ВшсНоп) ш Фе погта|те4 зрасе, 4ейпе4 засВ аё Н(и) < 0 парПез$ гайитге, зее Е1оиге 15. ш Фе рВуз1- 
са] зрасе, фе птазе оЁ Н(и) 1$ Ше Пой зе РапсНоп С(х,у), зее Еизге Та. Озше е с!а5$1са]! арргоасЬ, фе КВРО 
ргосез$ 15 сагие4 оц ш \о 5расез. ТВаф |еа4$ №0 а 126 сотршайопа! ите ргоет. ТБегеЮге, Беге 15 а гоп пее4 
тю деуеор еРНслет{ те#од$ [9]. 


М№ттайзе4 $ расе 












Еайиге Еайиге 
Кезоп Кезоп 
Дезеп Н(и)=0 


Рой 
А 


и, 


а) Ь) 
Е. 1. РВуз1са| ап4 погта|2е4 зрасез 


ш фе Бе 4 оРтейпабИиу-Базе4 4ез1еп орнпихайоп, ме @зпеи$Ь Бебмееп мо {урез оРГуапаез: 


1. ТЬе орбпитайоп уамаШез, \Б1сЬ аге деегили5Яс уапаез ю Бе а4лазе4 уф а меу/ ю орйпиите Фе 
$17108; Феу гергезеп{ е сопо| рагате{егз оЁ фе тесКатшса!| зузет (1.е. Читепз1оп$, таёепа|5, [оа4$, ес.), 
ап Фе ргоба зас тоде! (1.е. теапз ап4 запдаг4 деулайот$ оЁтапдот уапаб]ез); 

2. ТВе гапдот уанаШез, уурлсЬ гергезеп! Бе ипсецашиез ш Фе зузет. ЕасЬ о# фезе уапаез 15 14епиВеа 
Бу Ше буре оё а1561Бийоп 1а\ ап Фе аззос1ае4 рагатеегз. ТБезе уапаб]ез тау Бе фе сеотейлс Айтеп- 
$1018, Ше свагасет$с$ оРе шаета! ог Фе ежегпа] ]оа4$. 


3. ОрИишиш За{еёу Еасог (ОЗЕ) 


ТВе Рагйа| ЗаРеёу Еасфог$ (РЗЕ) ргезеще4 т [10] изе Ше са|бгайоп те о4$ {ай пее4 №ю ргорозе зоте соп- 
эта 5 Чите Фе саНбгайоп ргосезз ю шсгеазе фе еЁйслепсу ап Фе ассигасу. ТБе гезпе зоайоп уБеп и5те 
РУЕ тау по! гергезеп{ а 21офа| ог еуеп 1оса] орйтит. Ц тау зайзЁКу Ше гедите4 гейабИиу [еуе| Бесаизе о# Фе еЁЙ- 
с1епсу оР\е изе4 орйпиханоп а]еогИт. Ап еЁйсет ОЗЕ тео еззепнаПу 4ерепа$ оп Фе зай РасНоп оЁ фе орй- 
ташу сопаюопз оЁ е гепабИиу шдех рго ет (2). ТВ1$ тео4 ргоу!4ез Ве 4ез1епег а [еа5ё \/ В а 1оса| тепабИиу- 
Базе оритит Иво ад Чопа! сотрийпе со$. ТЬ15 тео4 ваз Бееп Баз1саПу 4еуе1оре4 Юг а погта| 51 иНоп 
сазе [11]. п 1$ \отК, И 1$ вепега|те4 10 Бе аррПе4 © а зте апа пиаре РаЙаге сазез. 


3.1 шее ЕаЙаге Моде (ЗЕМ) 
Тве ЗЕМ геНаб у ргоет сап Бе ууп@еп аз: 


ВМ = пи 4(м,) = и + +... 36: Н(иьи,,.ьи,)<0 (3) 


уВеге 4(и,) 15 Фе пшипит 415апсе Бебмееп Ше 4ез1оп ро ап4 Фе орйта! зомНоп. Ава Н (и,) <0 гергезет Фе 


аЦаге тоде. ТВе сотгезроп4 т? апа]уйса! ЮгтиаНоп изте ОЗЕ сап Бе утШеп аз ЮПо\з: 





=... п (4) 
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уВеге Ше $121 оЁ- 4ерепа$ оп Фе з12п оЁ\е 4епуайуе, 1.е., 


а в ОЕ 1=1,...й 


7 7 


Еогтиаяоп (4) ргоу14ез Чегет оритит уа[аез оЁ е погтаПте4 уама ез ай {йе 4ез1еп рошЁ ап4 аКше шю ас- 
соипё а зпе @Иате тоде. Ш [2], а зиаПаг ЮгтиЙавоп сап Бе Юипа Юг $еуега| 15а оп8. 


3.2 Мшар ЕаЙиге Моде (МЕМ) 


Озшя Фе зате ОЗЕ деуеортет Юг а зте РаЙиге тоде [2], ше МЕМ ргоет сап Бе \упщеп аз: 


Ви = ий (м, = + и За Н (ищи, )<0 (5) 


2^т 


Н, (и,) < 0 гергезепе фе @1егеп{ аИоге то4ез. ТЬе сотгезроп@ т? апа[уйса! ЮгтиаНоп изте ОЗЕ сап Бе упиеп 


аз ЮПо\м: 


(6) 








Еогищайоп (6) ргоу14ез Чегет оритит уа[аез оЁ 1е погтаНте4 уама ез ай {йе 4ез1еп рошЁ апа {акте ш® ас- 
соипЕ зеуега! @Паге то4ез. 


3.3 ОЗЕ вот 


Те Орйтит Заеёу Еасюг (ОЗЕ) а]еогИ бт сап Бе еазПу пар!етещеа ш @гее рипстра| 5%ерз (Е1е. 2). Те 
Нгзё %ер 1$ ю деегтте е 4ез1юп рошё соп$14еппе е 110$ асйуе сопзгашЕ аз а Пти зе ВапсНоп С(у). ТВе орН- 
шиайоп ргоет 1$ ю пишите фе об]еснуе Рапсйоп зиб]ес ю фе Пти з(айе ап Фе дебегилил15Яс сопзбгаш. Тре 
тезите зо оп 1$ {егте4 Фе 4ез1еп рошё. ТБе зесопа %ер 1$ ю сотрие Ше заЁебу Рсюг$ изше Фе едиаНоп$ (4) 
апа (6). Тве 14 чер 1$ ю са]сище Фе орйта| зоайоп шсште фе уаез оЁ Фе завейу асюгз ш Фе сотпршайоп 
оЁ фе уащез оРШе 4ез1еп уапа ез ап4 еп 4ееглте те орйтит 4ез1еп оР Фе затасвге. 







5етаройи: Х,,Уъ- 





Ета!шае ОБзесйуе -уипсйоп 


з 
з 
я. = 
58 
5 
за 
> 
А 







Орйитгайоп Еогтшайоп 


УЕ$ 


Сотрше О5Е 


Меру 2 &3 


Вееущиие Еипсйоп$ 


Орйтит Ро 





Ее. 2. ОЗЕ аюоги та 
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4. Митенса! Арр|сайоп 


Тве обзеснуе оЁ 1$ аррПсайоп 15 0 зВо\ фе адуащасе оР Фе ВВГО Бу ОЗЕ ге!айуе №0 Фе ООО узВеп 4еа]- 
ше м1 ЗЕМ апа МЕМ. ТБе сотгезроп4т» таепа| ргорег@ез аге: Уоипз”5 тодиз Е = 206.8СРа, Ро1550п’5 гайо 
у= 0.29 апа аепзиу р = 7820 Ке/и?. Тке епдигапсе Шиайз ог фе теуегзе4 1епзюп $655 ап ‘огзоп 565$ #-1 апа 1-1 
убайе4 айег 2.106 сус1ез аге едиа] © 252МРа апа 182МРа, гезресиуу. Еог тоге 4ейа$ або фе @иеие даа апа 
тшешфо4$ оё 111$ ехатр[е, Ше пиеге$е4 геадег сап геРег фо [12]. ТВе 1епе® ап Фе Вет оЁ Фе зваЧе4 рае аге: 
Г=0.14 т апа Н=0.1 т, гезреснуе[у (Е1о. За). ш Япие е|етеп{ апа]уз1$, ве р1ае 1$ зиррозе4 зиир|у Яхе4 оп 5 Юг 
е4еез ап4 1$ шо4деПеа Бу 32 е1юМ-по4е здиаге е]етеп{$ \фисЬ ргодасе по ощ-оЁр1апе 5гезз (Ее. 3Ъ). Неге, ве 
{сКпез$ оЁРеасВ @етепи Т, 1$ сопз14еге4 аз а гапдот уапае ап4 15 теап аз а дйегииизЯс уапаЫе {та{ [еа4$ 10 32 


т 
5 5.5 Б 5.5 г 7.5 8 8.5 


х103 


деегиаиитзИс уама ]ез ап 32 гап4от уапа ез. 











а Ь 
Ее. 3: 
а - Оипеп$101$ оЁ Фе реггайе4 р]ае ап4 
Ь - Тыскпез$ 15а Ноп Юг геза т орбит 4ез12п 


Тве об]еснуе оР Фе ОПО ргоседге 1$ №0 штитияе Фе уоате заб уесё 10 Ше Райгие датазе сопзгали 5. Неге, \уе 
сопз14ег а 1офа| зафебу асюг 5, =1.25 аррИе4 0 Ше датазе Р,‚„ап@ Базе оп Ше епотеегие ехрепепсе. ТВе 
ВВПО ргосе4иге саппоё сопёго] поё оу фе гепабИиу 1еуе!| Би зеуега| оиёри{ рагатеегз. ш [2], фе гепабПиу 1еуе1 
Баз Бееп сопгоПе4 соп$14етие а фагое{ гепабИиу [еуе], Бо\уеуег, Веге \ме зееК №0 сопго] е офег ошёриЕ рагате{ег$ 
зисН аз фе згасга| уоТате ап4 Ве датазе сгИепоп. 

Тве ООО апа ВВПО гези 5 ш {аЫе (1) зБо\и фа Фе ООО саппо! ргоу14е фе дез1епег \/ИВ а гедите4 гепаб!- 
пу 1еуе! мЮШе ВВПО Ъу ОЗЕ аПо\ соптоШпя фе заёейу 1еуе!5. Ассог4те {10 Фе ргоепа$ 3 ап 5, е уаше оЁ Фе 
21ора| заеёу Гасог 13 аррПе@ 10 Ше иррег Чатазе ший Ду =1 ю е5, =1.25. ТЫ$ мау Ше аПо\уаЫе Чатаре Ш 


Бе: 2, =0.8. Тне уап4агА деманопз аге соп$14еге4 аз ргорогйопа| о# {пе теап уашез: о; = 0.57и, , {=1,...,32. АЙег 
Вауше орйтите4 Фе знасеоге, фе гези шт? ООО уоше (ТаЫе 1) \уаз юиа4 ю Бе Ирро=105.64спт. Тне сотезропа- 
ш8 генаб у шдех уаз КЮюцпа 10 Ъе: Во =2.73. ТЬ15 гези ще уаше 4оез поё Ъе1опз №0 Ше 5апдаг4 збгасвага| епэ1- 


пеегиё пог Во = 2.73 [3—4.25] . Но\еуег, Юг фе зате оритит уоате, е КВРО Юг шшШар!е ЁаПиге то4дез 


(изше 5 апа 6) ргоу14ез 1е ез1епег \/ИВ а тоге тепае оритит згасвге уу В = 3.64. Еюмге$ 4а апа Ъ $Во\ 


АВОО” 
Фе пуегуа! оРШе датасе 915 1БаНоп оРаП зигасеаге @ФисКпез$ ог ООО апа ВВПРО 1 Ше зате уо!ате Грро = Гввро.. 
Тве гезите Чатазе а1зи1бийоп пиегуа! ор ВВОРО Бу ОЗЕ [0.52,0.8] 15 Бейег Фап Фе гезип» опе Бу ООО 
[0.4,0.8]. \УЬе Епзиге$ 46 апа с зВо\ Ше ищегуа[ оР Фе датазе ЧзиТиНоп оЁаП згасваге @сКпез$ ог ООО апа 
ВВОО ог Ше зате тахипии датасе уашез Орро=Рьвро =0.8. Тве ВВОО Ъу ОЗЕ рго\1Аез Ше 4ез1епег \и ап 


оритит з@гасваге \уйй ап шсгеазе 7% Биё тоге генаШе Бу 45% ге!ануе ю Ше гези из згисёаге Бу РОО (Ви» = 


3.78>Вььо =2.73). Неге, Ве гези т; датазе зил иНоп пщегуа| оЁВВРО Ъу ОЗЕ [0.56,0.84] 1$ а!зо Бецег ап Ше 


тез те опе Бу ООО [0.4,0.8]. Тва$, уВеп оМашише а БеНег датазе Ч1з61ийоп пиегуа|, Фе уо[ате 1$ гедисе4 ш 
Фе поп-спИса] згасбага] гез1оп5 Фа{ [еа45$ 0 есопопис ап4 гейае 4ез1еп5. 
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ТаЫе 1 
ООО ап4 ВВРО гезиИ$ фог ВМ 






































































































































а Реп Орйтит ЗоаНоп$ 
Роше ВВРО“” рро ВВОО” 
Уо[ите 81.9770 105.53 105.64 112.70 
Бутх 0.9945 =1 0.84 0.8 0.8 
Вы 0 3.64 2.73 3.78 
Т, 5.5784 8.914 7.7839 9.665 
Т> 7.7249 11.477 9.2979 12.434 
Тз 8.7796 8.340 10.1057 8.732 
Т. 8.7876 8.328 10.1055 8.718 
Т; 7.7435 11.571 9.2976 12.536 
Тв 5.604 9.001 7.7838 9.759 
Т, 6.8081 12.024 8.6421 13.129 
Тв» 6.0581 8.675 8.1291 9.265 
То 5.3822 3.715 7.5599 3.738 
Тю 5.3866 3.745 7.5597 3.770 
Ти 6.0644 8.678 8.1290 9.266 
Тр 6.8267 12.063 8.6422 13.169 
Тв 7.605 10.503 9.2265 11.214 
Та 4.7228 4.081 7.1685 4.157 
Тв 4.7194 4.061 7.1689 4.136 
Тб 7.608 10.507 9.2268 11.218 
Ти 7.6074 10.504 9.2267 11.214 
Тиз 4.7191 4.057 7.1688 4.131 
Тю 4.7231 4.082 7.1685 4.159 
То 7.6056 10.509 9.2264 11.220 
Ту 6.8256 12.057 8.6421 13.162 
Т22 6.0637 8.675 8.1290 9.262 
Т2з 5.3861 3.738 7.5597 3.762 
Та 5.3825 3.721 7.5599 3.745 
Т5 6.0588 8.680 8.1291 9.270 
Т5 6.8093 12.031 8.6421 13.137 
Ту 5.6029 8.992 7.7838 9.749 
Т28 7.7425 11.557 9.2976 12.521 
То 8.7871 8.311 10.1055 8.699 
Ту 8.7799 8.351 10.1057 8.744 
Тя 7.7257 11.491 9.2979 12.450 
Е Тз2 5.5794 8.922 7.7839 9.674 
3 
ы Те эёгасвага] гепаб Ищу 1еуе! 15 ппргоуе4 Бесаизе фе ОЗЕ-Базе4 зоайоп еззепнаПу 4ерепа$ оп е зеп1Нуйу 
8 эби4у з/шсЬ Чеегплше$ Ше го]е оЁеасБ рагатеег г@айуе 10 фе РаПиге ргобабИиу. Рог Фе сотрийпе ите, фе РОО 
2 ргоседиге пее4$ ю зо|уе мо зедиепйа] орбпитайноп ргоегаз (1) апа (2) ме Фе ВВОО Бу ОЗЕ сап геаП2е фе ор- 
Е егайоп ш ошШу опе зте орйпитаноп рго ет. 
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с) Оршиит ЗошНоп Бу АВРО?, ОЭшах=0.8, Уоите=112, В =3.78, ОЕ [0.53,0.8] 


Волге 4. Кезатс датазе 41561Баноп ищегуа! Бу ООО апа ВВОО ргоседигез 


ТаЫе 2 ргезет фе Чегет ОЗЁЕ гезиН$ соп$1Аетие бе Ппеаг ап4 попПпеаг @1зи1Бийоп$ (погтла], 1о2погтла|, 
ип огт, УеПЬо ап4 Ситбе! 415 ФиНопз), зее [2], заизРуте а гедите4 гепаб у 1еуе1 В =3 . ТВе гедите4 гена - 
Ну [еуе| сап Бе сопз14еге4 аз а отуеп 4айа ап пе а1гог т сопуегоез ю фе орйтит 4ез1еп а{ уетЙез Ше гедиге- 
тет (соптоПаЫе 4ез1етз: гепабИиу, уоите, Чатазе ...). Рог фе сошрийпз Яте сопзитрНоп, Юг а зше аИиге 
тоде, УТеп и$те ОРО ргоседиге, 5\о орНпитаНноп ргосеззез аге азе4 (Фе Игз{ 1$ ю Нпа Фе орнтит зоНоп апа 
Фе зесопа 1$ ю Япа Ше 4ез1еп рош®. ТВе рго ет Бесотез тасВ тоге сотр!ех Юг шаШр!е РаПиаге то4ез УТеге 
еасн Цаге то4е пее4$ а зерагае орйпитаНоп ргосез$ 10 Ипа Фе сотгезропате 4ез1еп рош!. Но\еуег, фе ВВОО Бу 
ОЗЕ пее4$ ошу а зше орипитайоп ргосез$ фо Йпа фе 4ез1еп рошё ап4 пех Фе орйтит зо[аноп 1$ апа1уйсаПу 
сотршеа изте ОЗЕ-ЗЕМ ог ОЗЕ-МЕМ Югтаайопз. ТВе КВОО Бу ОЪЕ 1$ Шеп саге оиё \уИКоиЕ аЧ4юопа] сот- 
рийп® Яте Бесачзе и Баз а зпе уапае уесюг фа{ 4ейпез Фе 4ез1оп роше. 


ТаЫе 2 
Глпеаг ап попПпеаг ВВОО гезий юг геаште4 тепабИиу ш4дех В=3 





















































Эфтисага| Рез1еп Орйтит КВОО зов оп$ 
Рагатейег Роше Мопта| Горпогта! Оп огт УУе!ЬиП Сите! 
Уоците 81.97 86.63 86.92 88.22 85.59 83.28 
Бах 0.99 =1 0.94 0.93 0.92 0.95 0. 98 
В 0 3.00 3.00 3.00 3.00 3.00 
Р: 50% =0.1% =0.1% =0.1% =0.1% =0.1% 
Т, 5.5784 5.8679 5.8890 5.9694 5.7954 5.6604 
Т› 7.7249 8.1947 8.2203 8.3520 8.0980 7.8541 
Тз 8.7796 9.3409 9.3684 9.5255 9.2330 8.9342 
Т. 8.7876 9.3488 9.3763 9.5334 9.2407 8.9428 
Т. 7.7435 8.2130 8.2387 8.3704 8.1160 7.8739 
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Те 5.6040 5.8933 5.9145 5.9949 5.8204 5.6873 
Т) 6.8081 7.1965 7.2204 7.3293 7.1097 6.9176 
Тв 6.0581 6.3822 6.4046 6.4950 6.3039 6.1518 
То 5.3822 5.6633 5.6836 5.7618 5.5934 5.4579 
Тю 5.3866 5.6677 5.6880 5.7662 5.5978 5.4625 
Ти 6.0644 6.3885 6.4109 6.5013 6.3101 6.1584 
Тр 6.8267 7.2148 7.2389 7.3476 7.1277 6.9374 
Тв 7.6050 8.0591 8.0848 8.2121 7.9634 7.7337 
Ти 4.7228 4.9483 4.9670 5.0288 4.8862 4.7866 
Тв 4.7194 4.9451 4.9638 5.0257 4.8831 4.7830 
Тб 7.6080 8.0621 8.0877 8.2151 7.9663 7.7369 
Ти 7.6074 8.0615 8.0872 8.2145 7.9657 7.7362 
Та 4.7191 4.9448 4.9635 5.0254 4.8828 4.7827 
То 4.7231 4.9486 4.9673 5.0291 4.8864 4.7869 
То 7.6056 8.0597 8.0854 8.2127 7.9640 7.7343 
Тр 6.8256 7.2137 7.2378 7.3466 7.1267 6.9362 
Т>2 6.0637 6.3878 6.4103 6.5007 6.3094 6.1577 
Т»з 5.3861 5.6672 5.6875 5.7657 5.5973 5.4620 
Та 5.3825 5.6636 5.6838 5.7620 5.5937 5.4583 
Т.5 6.0588 6.3828 6.4052 6.4956 6.3045 6.1525 
Т.5 6.8093 7.1976 7.2216 7.3304 7.1109 6.9189 
Ту 5.6029 5.8922 5.9135 5.9938 5.8193 5.6861 
Тв 7.7425 8.2121 8.2378 8.3695 8.1151 7.8729 
То 8.7871 9.3483 9.3759 9.5329 9.2402 8.9422 
Тзо 8.7799 9.3412 9.3687 9.5257 9.2332 8.9345 
Ти 7.7257 8.1954 8.2210 8.3527 8.0987 7.8550 
Тз2 5.5794 5.8689 5.8899 5.9704 5.7964 5.6615 





























Тве ОЗЕ тео 15 звоу\уп аз а Ч1зНпсйуе ю01 Юг КВОО ргоеплз.  зВо\у5 фе ЮПо\лпя а4уащаез: 

® — ТБе оМашеа гейабИиу-Базе4 орйтит 5о[айопз зВошА Бе тоге гепае ап фозе оМаштеа Бу ООО ргосе- 
Чиге Юг Ше зате орйтит уо[атез, 

® — ТЬе ОЗЕ ргоседиге пее4$ ошу а зше орнпихайоп ргосезз Юг фе 4ез1еп ро мо а@Ч1опа] сотрийпе 
ите Бесаизе й Ваз а зшее уапае уесюог а 4ейпез Фе 4ез1еп рошЕ Ве Фе ООО ргоседиге пее4$ мо 
орипихайоп ргосеззез. 

е  Эшсе фе шарог АЁЯсиМу Пез ш 1е еуалайоп оЁ Фе ргобаыИз$Нс сопзгаиз, уБсЬ 1$ рговауеу ехреп- 
уе ап еуеп Чуегоез ог тапу аррИсайопз, фе ОЗЕ ргоседиге рго\1Чез фе 4ез1епег \/ИВ ап апа[уйса1 
еуашаНоп \/ИВ зтаП сотрийпе Ите ге!ануе ю ОО ап оег ВВОО ргоседигез. 

е  АШтеПаБИИу шдех суашайотз Юг ВВРО-МЕМ з41е$ сап Бе апа[уйсаПу сагме4 опё Рог Чегет ргобаб!- 
15ис зи опз. ТВеге 15 по пее4 ю орнитаНоп ргосез$$ез. 


7. Соп $1015 


Тве ВВРО изше ОЗЕ Баз 5еуега| а4уатасез: зтаП питфег оЁ орбпитайоп уапаез, 2004 сопуегеепсе 5ваб!- 
Ну, зтаП сотрийпе ите, зайРасНоп оЁ Фе гедите4 гепабПиу 1еуе]5 апа 21офа! орйта ап4 тоге есопопис зоайоп 
Юг Ше зате геПабИиу шдех геапуе 10 Фе ООО ргосезз. ТЬе ОЗЕ-МЕМ Югт\щаЯоп сап Бе аррПе4 еазПу №ю ЧИегеп 
ЧзтичНоп 1а\уз. ш ее апа[уз1$, Фе зе оРе с[а551са] плео4 1еа4$ 10 ап ехиетеу 612В сотрийп? Ите ап4 а[зо 
оса] орйтта. 
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Механические свойства сервовитных пленок, формирующихся при трении в водных 


растворах карбоновых кислот" 


1 2 3 4 
В.Э. Бурлакова, Е. Г. Дроган’, А. И. Тюрин’, Т. С. Пирожкова”” 
1? Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 
34 НИИ «Нанотехнологии и наноматериалы» Тамбовского государственного университета им. Г. Р. Державина, Тамбов, 


Российская Федерация 


Месвашеса! ргорег#е$ оЁ зегуоуйе И!т$ Гогшед ш дигто Й1евоп адиеоц$ зоб оп$ о? сагБохуЦс ас1а$ ^^" 


У. Е. ВиНаКоуа', Е. С. Огозап?, А. 1. Туигию, Т. $. Рио7ВКоуа“" 


12 роп За Тесшиса] Ошуегзйу, Возюу-оп-Ооп, Киз1ап Еедегайоп 


34 Везеагсв шзниие оЁ Мапо{есНпо|021е5 апа Мапотаента15, С.В. ег2Ваут Татфбоу Зе ОшуегзИу, Татфоу, Визз1ап Еедеганоп 


Введение. В работе показано, каким образом природа орга- 
нической компоненты в систематическом ряду однооснов- 
ных карбоновых кислот влияет на трибологические харак- 
теристики пары трения «латунь — сталь» в водных раство- 
рах. Изучена зависимость физико-механических свойств 
антифрикционных пленок, формирующихся при трении, от 
природы смазочной композиции. Цели работы: изучить 
возможности использования карбоновых кислот как анти- 
фрикционных компонентов смазочного материала; оценить 
их влияние на механические свойства сервовитной пленки, 
формирующейся при трении латуни по стали. 

Материалы и методы. Проведены трибологические иссле- 
дования пары трения «латунь — сталь» на машине трения 
торцевого типа АЕ-5. Параметры шероховатости сервовит- 
ной пленки определялись с помощью оптической профи- 
лометрии. Микрогеометрия и структура объекта на нано- 
уровне исследовались с помощью атомно-силовой микро- 
скопии. Механические характеристики антифрикционной 
пленки изучали с помощью инструментального наноин- 
дентирования. 

Результаты исследования. Изучены трибологические ха- 
рактеристики трибосопряжения «латунь — сталь» и физи- 
ко-механические характеристики сервовитной пленки, 
формирующейся при трении в системе «латунь — водный 
раствор карбоновой кислоты — сталь». Установлено, что 
при увеличении длины углеводородного радикала коэффи- 
циент трения снижается. Обнаружены размерные эффекты 
в механических и трибологических свойствах сервовитной 
пленки, формирующейся на поверхности фрикционного 
взаимодействия в водных растворах карбоновых кислот. 
Обсуждение и заключения. Результаты исследования пока- 
зывают, что при фрикционном взаимодействии на поверх- 
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ности трения в водных растворах карбоновых кислот фор- 
мируется наноструктурная сервовитная пленка, резко сни- 
жающая коэффициент трения. Ее физико-механические и 
трибологические параметры зависят от состава модельной 
смазочной среды. Определено, что локальные физико- 
механические свойства зависят от способа получения сер- 
вовитного слоя, нагрузки и размера зоны деформирования. 
Полученные результаты могут быть использованы при раз- 
работке смазочных материалов. 


Ключевые слова: коэффициент трения, избирательный 
перенос, сервовитная пленка, размерные эффекты, шерохо- 
ватость поверхности. 
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280-288 
Введение. Стремительное развитие современного машиностроения расширяет возможности конструи- 


рования машин, работающих в экстремальных условиях. Таким образом, возрастает актуальность вопросов по- 
вышения ресурса узлов трения при растущих требованиях к механическим свойствам материалов. 

Условия трения с повышенной плотностью энергии предъявляют к машинам и оборудованию ряд спе- 
цифических требований в плане трибологических характеристик, надежности и срока службы. Критические 
условия работы трибосопряженных поверхностей влияют на передачу силы, потерю энергии и в конечном сче- 
те — на динамическое поведение всего механизма. Пространственные распределения неровностей топографии, 
как правило, негативно воздействуют на эксплуатационные возможности машин и механизмов. Иными слова- 
ми, несущая способность узлов трения напрямую зависит от основных параметров состояния поверхности. Они 
влияют на условия взаимодействия в области контакта, адгезию и деформацию, тем самым генерируют вибра- 
ции в процессе трения [1] и обусловливают неравномерный износ. 

Наибольшую антифрикционную эффективность в данном случае демонстрируют смазочные 
материалы [2—6], содержащие такие металлические присадки, как медь, алюминий, серебро. Данные материалы 
в процессе трения формируют на трибосопряженных поверхностях защитные металлоплакирующие пленки. 
Таким образом обеспечивается низкий коэффициент трения и умеренно низкий износ в широком диапазоне 
контактных давлений и скоростей скольжения [2, 3]. Металлоплакирование поверхностных слоев в процессе 
эксплуатации позволяет «залечивать» дефекты поверхности (трещины, поры, накопленные усталостные повре- 
ждения). Это значительно повышает антифрикционные характеристики и улучшает физико-механические свой- 
ства, что особо важно для надежной работы узлов трения и увеличения срока эксплуатации оборудования. 

Стоит отметить, что формирование пленок переноса возможно при трении в водно-спиртовых средах 
без заметного суммарного износа контактирующих поверхностей [7, 8]. Применение избирательного переноса в 
машиностроительной практике открывает реальную возможность для проектирования долговечных и эконо- 
мичных узлов трения в подвижных сопряжениях деталей и узлов машин и механизмов [8]. 

В связи с этим интересно рассмотреть влияние состава смазочной среды на физико-механические и 
трибологические параметры сервовитной пленки, формирующейся в трибосопряжении «латунь-сталь» (в каче- 
стве смазочной среды используются водные растворы карбоновых кислот). 

Материалы и методы. Эволюцию коэффициента трения системы «латунь 59 — водный раствор кар- 
боновой кислоты — сталь 40Х» исследовали в лаборатории «Гибридные функциональные материалы на основе 
графена» научно-образовательного центра (НОЦ) «Материалы» на торцевой машине трения типа АЕ-5. Опыты 
проводились при следующих условиях: скорость вращения — 180 об/мин; осевая нагрузка — 98 Н; время ис- 
пытания — 10 часов. В качестве органической компоненты смазочной композиции использовали однооснов- 


ные карбоновые кислоты предельного ряда с общей формулой К — СООН (К = С,Н›и+). 
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Механические характеристики (твердость — Я и модуль Юнга — Ё) сервовитных пленок определялись 
методом инструментального индентирования [9, 10]. При этом использовались: 

— наноиндентометр №Мапое51 Р1афогт 3 (Мсгтотиета!5, ИК) с подключенным функциональным блоком, поз- 
воляющим осуществлять нагрузки от 0,01 мН до 500 мН; 

— нанотрибоиндентометр 7/-950 ТиБот4етег — многофункциональный комплекс для механических испыта- 
ний материалов посредством динамического наноиндентирования. 

Методика наноиндентирования позволяет исследовать комплекс физико-механических свойств тонких 
приповерхностных слоев твердых тел и пленок толщиной до нескольких десятков нм [9-15]. Механические 
характеристики изучались с использованием алмазного индентора Берковича. Учитывая характерную толщину 
исследуемых пленок (от нескольких сотен нанометров до единиц микрометров), эксплуатационные физико- 
механические параметры измеряли на наноуровне при глубинах отпечатков индентора от 20 нм до единиц мик- 
рон. 

На наноуровне коэффициент трения исследовался на Гиротаетег Т1-950 при одновременном прило- 
жении к индентору нормальной и латеральной нагрузок [11, 16, 17]. Непрерывно регистрировались нормальная 
и латеральная составляющие сил (Ру и Ё/) и смещений (ймуи Й,), реализуемых при воспроизведении трибокон- 
тактов индентором. Регистрируемые данные анализировались, что позволяло судить о процессах трения и из- 
носа. Мерой оценки служил коэффициент трения (и. 

Для определения толщины сервовитной пленки и параметров шероховатости использовали оптический 
профилометр Сопюиг СТ-КТ с аналитическим программным обеспечением Из1оп 64, установленный в НОЦ 
«Материалы». Измерения проводились методом вертикальной сканирующей интерферометрии (75/) со скоро- 
стью сканирования 0,1 мкм/сек при повторяемости АМ5 0,01 нм. 

Топография поверхности сервовитной пленки изучалась с помощью атомно-силового микроскопа 
(АСМ) РНУТЕ Сотрас! в полуконтактном режиме монокристаллическим кремниевым зондом с алюминиевым 
покрытием. 

Результаты исследования. Довольно хорошо изучена [8, 18] сервовитная пленка, самопроизвольно 
возникающая в ходе фрикционного взаимодействия в растворе глицерина в сплаве «медь сталь». Известно, 
что она формируется на поверхности трения в топографически неравноценных условиях. Поэтому в разных 
местах фрикционного контакта состав и свойства сервовитной пленки различны. Ее толщина не превышает 
2 мкм, а механические, физико-химические и трибологические свойства отличаются от свойств обычной ком- 
пактной меди. 

Следует отметить, что одним из продуктов трибоокисления глицерина [19] является карбоновая кис- 
лота. В этой связи интересно рассмотреть влияние состава смазочной среды на физико-механические и трибо- 
логические свойства трибосопряжения «латунь-сталь» при трении в систематическом ряду одноосновных кар- 
боновых кислот. 

Длительные эволюционные трибологические исследования пары трения «латунь — сталь» в водных 
растворах карбоновых кислот выявили следующую зависимость трибологических характеристик от состава 
смазочной среды: длина углеводородной цепи радикала кислоты увеличивается от Су до С; — коэффициент 
трения снижается до 0,007. 

Причем наименьшее значение коэффициента трения, характерное для систем, реализующих эффект 
безызносности [20, 21], достигается в водных растворах валериановой и капроновой кислот. Одновременно с 
этим износ трибопары «латунь — сталь» снижается до 25 раз. В результате избирательного переноса при тре- 
нии на трибосопряженных поверхностях формируется медная пленка с различной шероховатостью и плотно- 
стью покрытия поверхности. 

Поверхность трения сканировалась в условиях перехода в ряду кислот «муравьиная — уксусная — 
пропионовая — масляная — валериановая — капроновая». При этом оптическая профилометрия выявила, во- 
первых, снижение шероховатости сервовитной пленки. Во-вторых, была установлена зависимость шероховато- 
сти от начальной топографии поверхности испытуемого стального диска (его К. равно 118 нм) и состава сма- 


зочной композиции [22, 23] (рис. 1). 


Бурлакова В. Э. и др. Механические свойства сервовитных пленок, формирующихся при трении в водных растворах карбоновых кислот 
ВинмакКота Г. Е. апа йе ойегу. Меспатшса! ргорегйех о} 5егуоуйе тд ютте4 аиниз ртспоп т адиеоих зошпоп$ ор сагБохуйс ас145 








1,7мкм 
0,5 _ 
-0,5 | 
-1,5 


-2,5 
6,0 мкм 
4,0 


0 



































0,5 1,0 15мм 





аа 


0,5 1,0 15мм 


НР АЗ 

















0,5 1,0 1,5мм 


— Профиль сканирования; —Профиль шероховатости; — Профиль волнистости 





Рис. 1. 2)-визуализация, ЗО-визуализация и профиль поверхности из стали (а), (в) и латуни-59 (6), (г) после трения 
в водных растворах уксусной (а), (6) и капроновой (в), (г) кислот 


Шероховатость А„ на базовой длине (Г = 800 мкм) сервовитной пленки снижается при формировании, 
например, в водном растворе капроновой кислоты, до 69 нм. (Для сравнения: в водном растворе муравьиной 
кислоты данный показатель достигает 580 нм, уксусной — 401 нм.) Разница объясняется коррозионной актив- 
ностью среды. 

В растворах муравьиной и уксусной кислот на трибосопряженных поверхностях обнаруживается мно- 
жество царапин, неровностей и пор, а также участки с медной пленкой, образовавшейся по островковому меха- 
низму. При поперечном сечении поверхности трения выявляются глубокие продольные канавки, что свидетель- 
ствует об абразивном характере износа трибосопряжения [23, 24]. 

Более детальное сканирование поверхности с использованием метода АСМ и визуализация изображе- 
ния обнаруживает значительные повреждения при трении в растворах муравьиной и уксусной кислот. На 3ЗО- 
изображении они визуализируются в виде макроскопических темных полос, соответствующих областям интен- 
сивного фрикционного напряжения, сопровождающегося сильным абразивным износом трибопары (рис. 2). 
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Рис. 2. Топография поверхности сервовитной пленки, полученной при трении в водном растворе уксусной (а) 
и капроновой (6) кислоты 
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Сравнительный анализ морфологии поверхности трения при переходе к кислотам с большей длиной 
углеводородного радикала приводит к выводу о смене характера износа, модификации структуры тонкого при- 
поверхностного слоя трибосопряжения. Это является следствием массопереноса компонентов контактирующих 
тел, а также адсорбции нанокластеров меди из рабочей среды на контртело в результате трибоэлектрохимиче- 
ских процессов в зоне трения при реализации избирательного переноса. При этом на стальной поверхности при 
фрикционном взаимодействии в водных растворах капроновой и валериановой кислот образуется медная плен- 
ка из мелкозернистых нанометровых кластеров меди. Сформированный на поверхности слой достаточно плот- 
ный, с малым разбросом частиц по размеру (рис. 1, 2). 

К настоящему времени накоплено значительное количество экспериментальных данных [25—29], сви- 
детельствующих о существенном изменении механических свойств при уменьшении характерных размеров 
элементов структуры объекта менее 100 нм. Это позволяет предположить изменение механических параметров 
сервовитных пленок по сравнению с таковыми для меди в объеме [30]. Для выяснения этого факта проводилось 
инструментальное индентирование. 

Значения твердости Н и модуля Юнга Е определяли по методике Оливера-Фара [10] из характерных Р- 


й-диаграмм (рис. 3) в соответствии со стандартами [9]. 
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Рис. 3. Типичные Р-й-диаграммы для сервовитной пленки, полученной при трении пары «латунь — сталь» в водном 
растворе валерианой кислоты (а), капроновой кислоты (6) 


Проведенные исследования показали увеличение значений Н и ЕЁ пленок по сравнению с аналогичными 
характеристиками медного образца. 

Увеличение твердости Н и модуля упругости ЕЁ пленки, образующейся на поверхности трения, может 
быть связано с изменением механизма пластической деформации. Высокая концентрация вакансий по границам 
зерен при фрикционном взаимодействии способствует зернограничному проскальзыванию структурных эле- 
ментов относительно друг друга. При этом реализуется механизм недислокационной пластичности [31-34], что 
приводит к улучшению прочностных характеристик поверхностного медного слоя и увеличению пластичности 
при переходе в ряду кислот «муравьиная — уксусная — пропионовая — масляная — валериановая — капроно- 
вая». Из анализа зависимости величин Н и Ё сервовитной пленки от глубины пластического отпечатка й. также 
следует, что на твердость влияет масштабный фактор. Значения модуля Юнга остаются практически постоян- 
ными вплоть до глубин, сопоставимых с толщиной сервовитного слоя. Наиболее ярко это заметно на пленках, 


образующихся при трении в валериановой кислоте (рис. 4). 
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Рис. 4. Зависимость твердости (Н) и модуля Юнга (ЕР) от глубины пластического отпечатка (й.) (показано в 
полулогарифмических координатах) при индентировании сервовитной пленки, полученной в водном растворе 
валериановой кислоты (а), капроновой кислоты (6) 


Выявлен следующий размерный эффект: 

— величина Н растет с ростом Й, в области малых глубин; 
— при дальнейшем росте й. быстро уменьшается Н. 

Для других образцов размерные эффекты выражены слабее. Так, например, для пленки, сформирован- 
ной в капроновой кислоте, величины Н и Е остаются практически постоянными во всем исследованном диапа- 
зоне й. (от 23 нм до 1,1 мкм). Затем, с ростом йс, значения Н и Е уменьшаются. Очевидно, это связано с увели- 
чением размера зоны деформирования: на Н и Е начинают влиять свойства переходного слоя и материала под- 
ложки. Таким образом, в области больших глубин отпечатка (й.> 1,1 мкм) величины Н и Ё уменьшают- 
ся (рис. 4). 

Размерный эффект [35, 36] также обнаруживается при изучении зависимости коэффициента трения от 
прилагаемой нагрузки в нано- и микрошкале [11, 16, 17] (рис. 5). 
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Рис. 5. Зависимость коэффициента трения К сервовитной пленки от нагрузки в наношкале 


Вначале он выражается в резком падении коэффициента трения с ростом нормальной нагрузки, плоть 
до некоторого значения ЕР». Затем коэффициент трения плавно увеличивается с ростом Ах. Таким образом, 
зависимости А =Х(ЁЕу) имеют явно выраженный минимум при достижении некоторого значения 
Емопт [11, 16, 17]. 

Полученные экспериментальные данные целесообразно сопоставить с теоретической зависимостью ко- 
эффициента трения от нормальной нагрузки. Таким образом, становится очевидно их хорошее качественное 
соответствие для всех исследованных образцов пленки. При этом на величину Р’м»и влияет тип используемой в 
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качестве смазочной среды карбоновой кислоты. Так, например, для капроновой кислоты Ем» = 30 мН, а для 
валериановой Г» = 100 мН. 

Обсуждение и заключения. Исследованы механические свойства сервовитных пленок, формирую- 
щихся в паре трения «латунь-сталь». В нано- и микрошкале выявлены размерные эффекты в механических 
(твердости) и трибологических (коэффициент трения) свойствах. Показано, что исследованные свойства и ха- 
рактер их изменения зависят от типа используемой карбоновой кислоты. Уменьшение среднего размера кри- 
сталлитов в составе наночастиц пленки приводит к увеличению ее прочности и проявлению сверхпластичности 
при растяжении и сдвиге. 

Полученные в работе результаты позволяют сформулировать следующие выводы. 

1. При фрикционном взаимодействии пары «латунь — сталь» в водных растворах карбоновых кислот 
на поверхностях трения формируется антифрикционная сервовитная пленка, способствующая резкому сниже- 
нию коэффициента трения. 

2. Сервовитная пленка, формирующаяся при трибовзаимодействии пары «латунь-сталь», является 
наноструктурной. 

3. Трибологические параметры системы «латунь — водный раствор карбоновой кислоты — сталь» за- 
висят от длины углеводородного радикала карбоновой кислоты, используемой в качестве смазочной среды. 

4. Состав смазочной среды влияет на параметры шероховатости сервовитной пленки: при переходе от 
муравьиной к капроновой кислоте шероховатость поверхности пленки снижается в 4 раза. 

5. Локальные физико-механические свойства зависят от условий получения сервовитного слоя, нагруз- 
ки и размера зоны деформирования (глубины внедрения индентора). 

6. Обнаруженные размерные эффекты (зависимость определяемых величин от глубины отпечатка ин- 
дентора) подтверждают вывод о наноструктурном характере сервовитной пленки, формирующейся в трибоси- 
стеме «латунь-водный раствор карбоновой кислоты-сталь». 
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Введение. Статья посвящена исследованию энергосиловых 
параметров совмещенного процесса «прокатка — равно- 
канальное угловое прессование (РКУ)». Цель работы — 
определение возникающих усилий прокатки и прессова- 
ния при деформировании совмещенным способом. 
Материалы и методы. Был проведен прочностной расчет 
матрицы и эксперимент по деформированию образцов из 
алюминия марки АГЗ] 6063. В ходе опытов значения уси- 
лий фиксировались с помощью тензометрической стан- 
ЦИИ. 

Результаты исследования. Результаты прочностного ана- 
лиза позволяют утверждать, что данная конструкция мат- 
рицы пригодна для создания экспериментальной установ- 
ки совмещенного процесса «прокатка — РКУ- 
прессование», поскольку рассчитанного запаса прочности 
вполне достаточно для реализации прессования на пре- 
дельных условиях. Усилия прокатки на всех этапах де- 
формирования превышают соответствующие усилия прес- 
сования, что является необходимым условием совмещен- 
ного процесса. 

Обсуждение и заключения. Полученные результаты могут 
быть использованы при проектировании эксперименталь- 
ных установок, реализующих исследуемый совмещенный 
процесс. При этом тензометрический метод исследования 
прочностных характеристик пригоден и при использова- 
нии калиброванных валков. 


Ключевые слова: совмещенный процесс, прокатка, рав- 
ноканальное угловое прессование, усилие, прочностной 
анализ, тензометрия. 
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Введение. Повышение качества деформируемого металла является одной из главных задач при созда- 
нии нового или совершенствовании известного технологического процесса обработки металла давлением 
(ОМД). Вне зависимости от того, к какому типу (прокатному или кузнечному) относится разрабатываемый 
процесс, рекомендуется применять схему деформирования, реализующую чистый сдвиг. В этом случае затра- 
чивается минимальное количество энергии на деформацию металла, достигается максимальная и однородная 
проработка заготовки по сечению, обеспечивается закрытие внутренних дефектов. 

За последние два десятилетия активно разрабатывались способы обработки давлением, позволяющие 
получать заготовки с ультрамелкозернистой структурой. Как показал их анализ, они основаны на принципах 
реализации сдвиговых или знакопеременных деформаций. Отдельно стоит упомянуть процессы, представляю- 
щие собой комбинацию указанных принципов. Такие процессы позволяют реализовать особый вид обработки 
давлением, который называется «интенсивная пластическая деформация» (ИПД). ИПД реализуют такие 
методы, как кручение под высоким давлением [1-3], равноканальное угловое (РКУ) прессование [4—8], 
винтовая экструзия [9-12] и другие. Особого внимания заслуживает ступенчатая РКУ-матрица (или РКУ- 
матрица с параллельными каналами) [13-14]. Она позволяет реализовать сдвиговую деформацию при 
прохождении заготовки через ее каналы и одновременно с этим два знакопеременных очага деформации, при 
условии сонаправленности входного и выходного каналов. В сравнении с обычной угловой матрицей 
ступенчатая РКУ-матрица является энергосберегающей, поскольку позволяет реализовать большую степень 
деформации за один проход при одинаковом усилии. 

Рассмотренные выше способы ОМД относятся к дискретным способам прессования. Эти технологии 
производства пресс-изделий имеют ряд недостатков. Главные из них обусловлены: 

— дискретностью (прерывностью) процесса прессования; 
— наличием реактивных сил трения на контакте металла с контейнером. 

Таким образом, вследствие неравномерности деформации и высокой энергоемкости процесса прессо- 
вания ограничивается длина прессуемых изделий, снижается их качество. Применение схемы непрерывного 
прессования позволяет устранить эти недостатки. Данный способ принципиально отличается от предыдущих 
тем, что для деформации и экструдирования заготовки через отверстие матрицы используются активные силы 
трения. Длина заготовки при этом не ограничена. 

В последние годы разрабатываются так называемые «совмещенные» процессы ОМД, которые пред- 
ставляют собой сочетание двух и более процессов деформирования. 

Следует особо отметить важнейшую особенность совмещенного процесса ОМД: зачастую при его реа- 
лизации недостатки составляющих отдельных процессов снижаются либо полностью устраняются. 

Кроме того, в последнее время все большее внимание уделяется энергосберегающим технологиям, ос- 
нованным на использовании для деформации активных сил трения. 

Авторы данного исследования разработали совмещенный способ прокатки и прессования в равнока- 
нальной ступенчатой матрице, который по сравнению с обычным прессованием в равноканальной ступенчатой 
матрице снимает ограничения по размерам исходных заготовок [15]. Он представляется одним из наиболее 
перспективных и малоизученных в практическом плане процессов прокатки и прессования и по многим пара- 
метрам превосходит известные способы обработки металлов давлением. 

Суть данного способа деформирования заключается в следующем. Предварительно нагретая до темпе- 
ратуры начала деформирования заготовка подается к прокатным валкам. Первая пара валков за счет сил кон- 
тактного трения захватывает заготовку в зев валков, а на выходе из него проталкивает заготовку через каналы 
равноканальной ступенчатой матрицы. После того как заготовка полностью выйдет из каналов матрицы, ее за- 
хватывает вторая пара валков, которые также за счет контактных сил трения полностью вытягивают заготовку 
из матрицы. Таким образом, в данном случае прессование заготовок обеспечивается контактными силами тре- 
ния, возникающими на поверхности контакта металла с вращающимися валками. 

Ранее проводились различные исследования данного совмещенного процесса — как теоретические [16-— 
19], так и лабораторные (на установке с одной парой валков) [20-23]. В настоящей работе приведены результаты 
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исследования энергосиловых параметров, изученных в ходе разработки опытно-промышленной установки с 
двумя парами валков (рис. 1). 





Рис. 1. Экспериментальная установка для реализации совмещенного процесса «прокатка — РКУ-прессование» 


Прочностной анализ РКУ-матрицы. При реализации совмещенного процесса «прокатка — РКУ- 
прессование» использовалась экспериментальная установка, созданная на основе лабораторных прокатных ста- 
нов ДУО-200 и ДУО-250. Их энергосиловые характеристики позволяют деформировать алюминиевые и медные 
сплавы в холодном состоянии. 

При проектировании новой технологии деформирования на этапе создания экспериментальной уста- 
новки особое место занимает прочностной анализ. Он позволяет определить стойкость и жесткость рабочего 
инструмента, оценить его запас прочности, что даст информацию о его долговечности. Теоретические исследо- 
вания совмещенного процесса «прокатка — РКУ-прессование» были выполнены с помощью компьютерного 
моделирования в программе РЕЕОКМ. В результате установлено [17, 18], что самой нагружаемой деталью яв- 
ляется матрица для прессования. Это объясняется тем, что в зоне стыка каналов возникают наиболее жесткие 
условия деформирования. Поэтому для проверки пригодности выбранной конструкции матрицы проводился 
прочностной анализ. С этой целью была построена модель матрицы с углом стыка каналов 150 градусов по 
размерам, представленным в [20]. 

При реализации совмещенного процесса целесообразно подбирать ширину заготовки таким образом, 
чтобы при деформировании в матрице металл не контактировал с ее боковыми стенками. Это позволит суще- 
ственно снизить общее усилие противодавления со стороны матрицы. Однако полностью исключать такой слу- 
чай нельзя. Поэтому в программном комплексе РЕЕОЕМ была создана модель РКУ-прессования, в которой 
высота и ширина заготовки соответствуют высоте и ширине канала матрицы. Данный случай деформирования 
является предельным и почти неосуществим на практике. Поэтому, если матрица окажется достаточно проч- 
ной, она подойдет и для совмещенного процесса, условия которого существенно «мягче». В частности, высота 
и ширина заготовки несколько меньше высоты и ширины канала матрицы, что является необходимым условием 
для прохождения металла по каналу. 

Начальные размеры заготовки — 10х40х90 мм. Условия и допущения, принятые при моделировании: 
— материал заготовки до деформации является изотропным, в нем отсутствуют какие-либо напряжения и де- 
формации; 

— заготовка разбита на 24 000 конечных элементов, средняя длина ребра элемента —1,5 мм; 

— начальная температура заготовки — 700 °С, также учитывались деформационный разогрев и теплопередача 
от заготовки к инструменту и в окружающую среду; 

— инструмент принят как абсолютно жесткое тело; 

— модель заготовки принята как упругопластическая; 

— материал заготовки — АТЗТ 1045, соответствующий стали 45; 
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— коэффициент трения на контакте заготовки и матрицы равен 0,1; 
— скорость движения пуансона равна 10 мм/сек. 
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В результате расчета определялось усилие деформирования на всех элементах матрицы с целью 
нахождения наиболее нагруженных элементов конструкции. Наибольшие усилия возникают в направлении 
действия нормальных напряжений (рис. 2). 
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Рис. 2. Графики усилия РКУ-прессования для боковых стенок (а) и центральных сегментов (6) 


Из рис. 2 видно, что самыми нагруженными элементами матрицы являются центральные сегменты. 
При деформировании на них возникает усилие до 260 кН, в то время как на боковых стенках оно не превышает 
75 кН. На основании этих данных был проведен анализ напряжений на инструменте в программе РЕЕОКМ с 
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помощью специализированного модуля. Принцип данного анализа заключается в том, что абсолютно жесткие 
элементы модели становятся упругими телами, воспринимающими напряжение. От заготовки проводится ин- 
терполяция действующих сил на инструмент. Затем задаются места фиксации, чтобы упругие тела под действи- 
ем сил не разлетались. В нашем случае местами фиксации являются кольцевые отверстия для крепежных эле- 
ментов. Помимо этого в данном расчете необходимо задать материал инструмента. Для изготовления матрицы 
была использована инструментальная сталь 5ХВ2С. Она применяется для изготовления: 
— штампов сложной формы, работающих с повышенными ударными нагрузками, 
— ножей для холодной резки металла, 
— резьбонакатных плашек, 
— пуансонов и обжимных матриц при холодной работе. 

В результате были получены следующие данные (рис. 3). 
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в) 


Рис. 3. Результаты прочностного анализа: напряжение на боковых стенках (а); напряжение в центральных сегментах (6); 
линейное перемещение (в) 


Из рис. 3 можно сделать следующий вывод: наиболее нагруженные зоны в центральных сегментах — 
это поверхностные области возле стыка каналов и дуговые зоны вокруг первых двух пар отверстий. Здесь воз- 
никает напряжение, доходящее до 250 МПа. 

В областях между отверстиями напряжение находится в пределах 125-150 МПа (рис. 3, а). 

На боковых стенках максимальные напряжения возникают в зоне контакта с деформируемым металлом 
и равны 95-100 МПа (рис. 3, 6). 

Для изучения жесткости конструкции были рассмотрены результаты расчета линейных перемещений 
элементов конструкции под действием упругой деформации в специализированной расчетной библиотеке 
«АПМ» [24]. Изучение упругой деформации показывает ее наибольшую величину в поверхностной области 
возле стыка каналов. Следует отметить два факта. Первый: действующие напряжения на всех элементах матри- 
цы значительно ниже предела текучести стали 5ХВ2С. Второй: максимальное значение линейного перемеще- 
ния достигает 0,019 мм (рис. 3, в), что пренебрежимо мало и позволяет говорить о повышенной жесткости 
предлагаемой конструкции. Учитывая эти два обстоятельства, можно сделать вывод о данной конструкции 
матрицы: она пригодна для создания экспериментальной установки совмещенного процесса «прокатка — РКУ- 
прессование», поскольку рассчитанного запаса прочности вполне достаточно для реализации РКУ-прессования 
на предельных условиях. 

Лабораторные эксперименты по изучению усилий. В результате проведенного прочностного расче- 
та определено место установки на матрице тензодатчиков для измерения усилия деформирования. Места креп- 
ления находятся на расстоянии 75 мм от переднего края центрального сегмента, где возникают максимальные 
линейные перемещения. 

Для измерения возникающих усилий деформирования использовалось следующее оборудование: 

— тензостанция ХЕТ-017-Т8 (ЗАО «ЭТМС», Россия); 
— датчики усилия с тензорезисторами ТКФО]1-2-200 (ЗАО «ЭТМС», Россия); 
— ноутбук для управления тензостанцией и записи сигнала. 

Тензостанция ХЕТ 017-Т8 производства ЗАО «ЭТМС» (г. Москва) предназначена для проведения из- 
мерений энергосиловых и иных параметров по нескольким каналам с временным разрешением записи до 
20 кГц [25]. 

При изучении энергосиловых параметров совмещенного процесса «прокатка — РКУ-прессование» 
планируется исследовать усилия деформирования на стадиях прокатки в валках и прессования в матрице. Для 
этого необходимо установить тензорезисторы на матрицу и на месдозы прокатной клети. Суть тензометрии 
заключается в том, что в процессе деформации на валки и матрицу воздействуют усилия прокатки и прессова- 
ния. В результате валки и матрица изгибаются. Вместе с ними деформируются тензорезисторы, прикрепленные 
к матрице и месдозам. Изменяется величина их электрического сопротивления. Тензостанция фиксирует это 
изменение и в соответствии с результатами предварительной тарировки датчиков автоматически преобразует 
его в график изменения усилия. 

Для предварительной тарировки тензорезисторов центральный сегмент матрицы и месдоза испытыва- 
лись на сжатие на гидравлической крутильно-разрывной машине МИ-40КУ. Использование крутильно- 
растяжной машины позволяет развивать значительные усилия (35—40 кН). При этом их значения фиксируются 
точно и надежно. Изменение напряженно-деформированного состояния под нагрузкой вызывает деформацию и 
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линейное изменение сопротивления тензорезисторов. Суть тарировки состоит в построении зависимости, свя- 
зывающей электрическое напряжение в схеме и прикладываемое усилие. Эта зависимость всегда имеет строго 
линейный характер. 

Для проведения тарировки исследуемые элементы последовательно нагружались с шагом 5 кН в диапа- 
зоне от нуля до 35 кН. Фиксировались соответствующие значения напряжения в схеме под нагрузкой и после ее 
снятия. В целях повышения точности выполнено три прохода по указанному диапазону - вверх, вниз, вверх (42 
измерения). В результате для каждого элемента получены зависимости «усилие — перемещение» («перемеще- 
ние» соответствует значению хода верхнего бойка). Данные тарировочных испытаний были статистически об- 
работаны, и по ним получены уравнения линейной регрессии, связывающие прикладываемое к инструменту 
усилие (Р‚, Н) с напряжением ((,, мВ) в схеме. 

Уравнение для измерения усилия прокатки имеет вид: 

Р;=-4852,70;- 23538. 

Коэффициент детерминации А” = 0,998. Стандартная ошибка измерения по результатам 42 тестов со- 
ставила менее 0,2 %. 

Уравнение для измерения усилия прессования имеет вид: 

Р‚=-—4120,20;- 21609. 

Коэффициент детерминации А” = 0,997. Стандартная ошибка измерения по результатам 42 тестов со- 
ставила менее 0,3 %. 

Полученные данные были занесены в программу регистрации и обработки измерений тензостанции 
ГЕТ-017-Т8. Это позволяет записать сигнал в виде графика усилия. 

После проведения тарировочных работ матрица с тензорезистором была размещена на эксперимен- 


тальной установке (рис. 4). 





Рис. 4. Измерение усилия при совмещенном процессе «прокатка — РКУ-прессование» 


Здесь, на месдозах первой клети также были установлены тензорезисторы для измерения усилия про- 
катки, которое выполнялось только на первой клети. Это ключевой параметр, согласно концепции процесса. 
При реализации совмещенного процесса усилие прокатки должно превышать возникающее усилие прессова- 
ния. 

В качестве исходного материала использован алюминиевый сплав 6063 (аналог алюминиевого дефор- 
мируемого сплава АДЗ1) при комнатной температуре. Его химический состав и механические свойства приве- 
дены в табл. 1. 












































Таблица 1 
Химический состав и механические характеристики алюминиевого сплава 6063 
Массовая доля элементов, % 
51 Ее Си Мп Ме Ст 7 т Прочие 
0,2-0,6 0,35 0,1 0,1 0,45—0,9 0,1 0,1 | 0,1 0,15 
Механические свойства при 20 °С 

Предел прочности, МПа Предел текучести, МПа Твердость, НВ Удлинение, % 

90 48 25 20 
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Заготовки одинаковых размеров (15х25х450 мм) прокатывались в первой клети до толщины 9 мм. Во 
второй клети зазор между валками был равен 7 мм. Длину заготовки выбирали таким образом, чтобы ее перед- 
ний конец захватывался второй парой валков, а задний конец в этот момент еще прокатывался в первой кле- 
ти (рис. 5). 


ОВО 
ОВО АВ 


ИИ 
ИОВА 
ОВО 


И 
| 


| 








Рис. 5. Изменение усилия при деформировании алюминиевого сплава 6063: прокатка в первой клети (а); 
прессование в матрице (6) 
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Данные замеров усилий сведены в табл. 2. 

















Таблица 2 
Результаты измерения усилия 
1-я стадия 2-я стадия 3-я стадия 
Этап прокатки, кН 205 217 242 
Этап прессования, кН 59 104 172 
Разница (запас сил трения), % 247 108 40 

















Также была определена разница между значениями усилия прокатки и прессования. В данном случае 
эта разница характеризует уровень запаса активных сил трения: 


А (жк _1).100%, 
ПРЕСС 
где Рирок — усилие прокатки в первой клети; Риресс — усилие противодавления в матрице. 

Сравнительный анализ значений возникающих усилий показал, что усилия прокатки на всех трех эта- 
пах превышают значения соответствующих усилий прессования. Таким образом, выполняется необходимое 
условие для реализации совмещенного процесса. При этом запас активных сил трения в первой клети с про- 
движением заготовки по каналам матрицы существенно снижается. 

Выводы. Результаты исследования свидетельствуют о том, что запаса прочности матрицы вполне до- 
статочно для реализации РКУ-прессования на предельных условиях. Следует отметить также, что усилия про- 
катки на всех этапах деформирования превышают усилия прессования, и таким образом выполняется необхо- 
димое условие реализации совмещенного процесса. При этом запас активных сил трения в первой клети с про- 
движением заготовки по каналам матрицы существенно снижается. 

Полученные результаты могут быть использованы при проектировании экспериментальных установок, 
реализующих исследуемый совмещенный процесс. Тензометрический метод исследования прочностных харак- 
теристик пригоден и при использовании калиброванных валков. 
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Влияние объема древесных опилок на механические свойства композитного материала" 


А. Г. Дьяченко, Т. п. Савостина”, Б. И. Саед 


ь? Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 


3 Университет Алеппо, г. Алеппо, Сирийская Арабская Республика 


ЕНесЕ 0Ё зазу ди уошише оп шесвашса! ргорег@ез оЁ сотрозНе та{ена! 


А. С. Эуаспепко', Т. Р. Зауоз@та?, В. 1. баед *" 


+2 роп Зе ТесЬса! Ошуегзну, Возюу-оп-Ооп, Ваззап Еедеганоп 
3 ОпёуегзИу оЁ Аерро, А!ерро, Зупап АгаБ Вер Ис 


Введение. Рассматриваются вопросы создания композици- 
онных материалов с добавлением древесных опилок. Дре- 
весная составляющая композитов обрабатывается водой 
для усиления шероховатости поверхности. Это увеличива- 
ет площадь контакта, вследствие чего усиливается взаимо- 
действие между опилками и волокнами полимера. Цель 
исследования — изучение возможности упрочнения компо- 
зиционных материалов, полученных из опилок. 

Материалы и методы. Образцы изготовлены из композитных 
материалов на основе ненасыщенных полиэфирных смол, 
армированных древесной стружкой. Обработка водой прово- 
дилась при комнатной температуре в течение 2, 4, 6 и 8 дней. 
Затем образцы проходили испытания на изгиб и сжатие. 
Результаты исследования. В результате испытаний образцов 
зафиксированы изменения их механических свойств. Опре- 
делено, каким образом сопротивления изгибу и прочности 
при сжатии зависят от времени обработки водой. Построе- 
ны графики, отражающие указанные зависимости. 
Обсуждение и заключения. После обработки водой компо- 
зиты, армированные древесными опилками, демонстриру- 
ют более высокую стойкость при изгибе. Это объясняется 
усилением шероховатости поверхности опилок и, как 
следствие, увеличением поверхности сцепления с компо- 
зитной основой. Кроме того, обработка водой повышает 
стойкость образцов при сжатии. Особенно заметно упрочняют- 
ся образцы, созданные на основе крупных опилок. Это объяс- 
няется образованием углублений на поверхности древесных 
опилок, что также усиливает адгезию между ними и поли- 
мерной основой композита. 


Ключевые слова: древесные волокна, опилки, полимер, 
композитный материал. 
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Работа выполнена в рамках инициативной НИР. 
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Введение. В современной производственной практике довольно широко используются композитные 
материалы на основе полимеров, армированных различными примесями. Основные преимущества примесей из 
древесины (опилки, стружка) — их приемлемая стоимость и доступность [1, 2]. 

Для определения соотношений волокон древесных опилок использовалась вибрационная установка с 
набором сит. В итоге были выбраны три объема опилок (5 мм*, 25 мм, 120 мм») [2, 3, 4]. В ходе экспериментов 
образцы подвергались трем видам нагрузки (давлению, сдвигу и изгибу) для определения влияния объема опи- 
лок на прочностные характеристики композитных материалов. 

В 2001 году образцы с древесным наполнителем подвергались опытам, целью которых было устано- 
вить качество сцепления между волокнами и связующим материалом. Результаты позволяют утверждать, что 
армирование древесными волокнами увеличивает сопротивление материала [3]. В том же году проводились 
эксперименты, доказавшие, что обработка глицероловой кислотой повышает прочность на разрыв образцов, 
армированных древесными опилками [3]. 

В ходе опытов 2002 года древесина и древесное волокно подвергались химической обработке с исполь- 
зованием МАРР (та|езце4 ро!ургору!епе — малеиновый полипропилен). Выясилось, что у таких образцов выше 
прочность на разрыв и изгиб [5]. 

В 2006 году проводились эксперименты по обработке поверхности древесины МаОН и ССН2СООН. 
Опыты показали, что такое воздействие позволяет достичь высоких показателей по прочности на разрыв [6]. 

В рамках данного исследования планируется продемонстрировать целесообразность проведения хими- 
ческой обработки для улучшения свойств композиционных материалов. 

Материалы и методы. Древесные опилки (объемом 5 мм?, 25 мм?, 120 мм”) [1] обрабатывались водой 
при комнатной температуре в течение 2, 4, 6 и 8 дней. 

100 граммов древесных опилок помещали в пластиковый контейнер, заливали 1 литром воды и остав- 
ляли на 2, 4, би 8 дней. Затем воду сливали и заливали массу дистиллированной водой на 1 час. По истечении 
этого времени опилки вынимали и сушили в печи в течение 24 часов при температуре 110 °С. Далее все иссле- 
дования проводились также при комнатной температуре. 

Результаты исследования 

Определение потери веса обработанных опилок 

Использовались по 10 г опилок разных объемов (5 мм’, 25 ммз, 120 мм"), обработанных описанным 
выше способом. После обработки образцы рассматривали под оптическим микроскопом. В этом случае потеря 
веса опилок равна разности веса до и после обработки. Соотношение потерь веса рассчитывается следующим 


образом: 


ПИ р, 
и 


где И» — вес образца волокна до обработки (г); И’, — вес образца волокна после обработки (г); Й’* — процент 


потери веса (рис. 1). 
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Рис. 1. Зависимость потери веса древесных опилок от продолжительности обработки ВОДОЙ 


График показывает, что опилки объемом 120 мм* теряют в весе не намного больше, чем опилки объе- 
мом 25 мм (разница — около 4 %). Если же сравнивать опилки объемом 5 мм’ и 120 мм*, то потери веса по- 
следних выше в 8 раз. При этом у всех образцов потеря веса напрямую зависит от времени обработки (чем 
дольше испытания, тем меньше вес). При воздействии воды в течение длительного времени в жидкость перехо- 
дит часть растворимых веществ, поэтому после промывания и высушивания материал становится легче. 

В дальнейших экспериментах рассматривались по три образца из композита. Для его получения опилки 
различного объема (5 мм*, 25 мм”, 120 мм*) пропитывались 10-процентным полиуретановым эфиром, хорошо 
перемешивали, добавляли затвердитель и выливали смесь в пресс-формы. Полиэстер составлял 90 % всей мас- 
сы образца. 

Испытание на изгиб композитных материалов 

Для проведения опытов были подготовлены образцы композитного материала 8х15х160 мм по АТМ 
790 [7, $8]. 

Сопротивление изгибу измерялось в трех точках. На рис. 2 показан образец композитного материала, 





армированного опилками, перед испытанием на изгиб. 





Рис. 2. Образцы композитного материала, армированного опилками, перед испытанием на изгиб 


Сопротивление образца изгибу (МПа) определяется по формуле: 
_ З.Р.Ё 
С > 2.В.Р ? 

где Р — сила нагрузки (Н), Г. — длина образца (мм), В — ширина образца (мм), { — толщина образца (мм). 

Испытание на сжатие композитных материалов 

Для проведения опытов были подготовлены образцы композитного материала @ 12,7х25,4 мм. Сопро- 
тивление измеряли в сжатом состоянии. Нагрузка прикладывалась до тех пор, пока не происходил разрыв. Об- 
разцы армированных опилками композитных материалов, подготовленных к испытанию, представлены на 
рис. 3. Сопротивление сжатию (МПа) определяется по формуле [9]: 


2 
где Р — сила сжатия (Н), Е — площадь поперечного сечения (мм). 
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Рис. 3. Образцы композитного материала, армированного опилками, перед испытанием на сжатие 





Результаты испытаний на изгиб 
3 р 
Из рис. 4 видно, что образцы, изготовленные из опилок объемом 5 мм”, более устойчивы к изгибу. 
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Рис. 4. Изменение сопротивления изгибу в зависимости от периода обработки образцов 


Максимальное сопротивление изгибу (230,88 МПа) продемонстрировал образец с опилками объемом 
5 мм? после обработки в течение двух дней. 

Образец 25 мм* демонстрировал максимальное сопротивление изгибу (98,6 МПа) при обработке в тече- 
ние 8 дней. Образец 120 мм° демонстрировал максимальное сопротивление изгибу (194,55 МПа) при обработке 
в течение 4 дней. 

Увеличение сопротивления изгибу в данном случае объясняется следующим образом. При обработке 
водой усиливается шероховатость поверхности опилок, увеличивается площадь взаимодействия с основой. Та- 
ким образом упрочняется сцепление между опилками и полиэфиром и, как следствие, увеличивается прочность 
образцов при изгибе. 

Результаты испытаний на сжатие 


Образцы из опилок объемом 5 мм* обладают меньшей прочностью на сжатие (в сравнении с образцами 
из опилок объемом 25 мм'и 120 мм? (рис. 5). 
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Рис. 5. Изменение сопротивления сжатию в зависимости от периода обработки образцов 
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Максимальное значение прочности на сжатие (156,67 МПа) демонстрируют образцы с опилками объе- 
мом 120 мм, обработанные в течение четырех дней. Это объясняется тем, что при обработке водой на поверх- 
ности более крупных древесных опилок образуются углубления. Это приводит к усилению адгезии между 
опилками и основным полимером и повышает устойчивость композита к сжатию. 

Обсуждение и заключения. Испытания образцов из композитных материалов, армированных древес- 
ными опилками и обработанных водой, позволяют сделать следующие выводы. 

1. Образцы из опилок объемом 120 мм? теряют в весе больше, чем образцы из опилок объемом 5 мм’ и 
25 мм’. 

2. Потеря веса образцов напрямую зависит от времени обработки водой. 

3. При обработке образцов водой первое резкое увеличение сопротивления изгибу наблюдается в тече- 
ние двух дней. 

4. Максимальное сопротивление изгибу (230,88 МПа) продемонстрировал образец с опилками объемом 
5 мм°. 

5. Максимальное сопротивление сжатию (156,67 МПа) продемонстрировал образец с опилками объе- 
мом 120 мм* 
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Введение. Статья посвящена изучению процесса электро- 
акустического нанесения износостойких покрытий. Целью 
работы является получение математической модели влия- 
ния ультразвука на процессы формирования и развития 
электрической искры, происходящей в процессе электро- 
акустического напыления. 

Материалы и методы. В основе анализа влияния ультра- 
звуковых колебаний на процессы, протекающие при фор- 
мировании и развитии искрового разряда, рассмотрены 
уравнения непрерывности, движения и переноса энергии с 
учетом вклада ультразвукового поля. Учтены факторы, 
влияющие на теплопроводность и электропроводность 
сильно ионизованного газа. 

Результаты исследования. При получении модели были 
сделаны предположения, что отвод тепла из канала осу- 
ществляется «прозрачным излучателем». Тогда для обла- 
сти канала было принято автомодельное решение: давле- 
ние, температура и плотность постоянны по сечению, а 
скорость пропорциональна радиусу. Получена математи- 
ческая модель, описывающая процессы, протекающие в 
искровом канале с учетом влияния энергии ультразвуко- 
вого поля. 

Обсуждение и заключения. На основании построенной 
модели установлено, что под действием ультразвука уве- 
личивается радиус и температура искрового канала, а 
также создаются условия двукратной ионизации при вы- 
соких энергиях ультразвука. 


Ключевые слова: 
упрочнение, упрочняющие покрытия, высококонцентри- 
рованные потоки энергий, ультразвук, электрическая ис- 


электроакустическое напыление, 


кра, математическая модель, токопроводящий канал, тер- 
модинамика, ионизация. 
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Введение. При нанесении износостойких покрытий методом электроакустического напыления образу- 
ется сравнительно узкий токопроводящий канал с высокой температурой и ионизацией, в котором выделяется 
джоулево тепло, что приводит к повышению давления и расширению канала. Расширяющийся канал действует 
подобно поршню на остальной газ и, так как расширение происходит со сверхзвуковой скоростью, вызывает в 
нем ударную волну, которая распространяется впереди этого своеобразного поршня. Температура в области 
ударной волны гораздо выше, чем в невозмущенном газе, а температура в самом канале во много раз больше, 
чем в ударной волне. Соответственно, плотность газа в канале очень мала, подавляющая часть массы движуще- 
гося газа вытесняется из него, что и дает возможность рассматривать границу канала как поршень [1]. 

Сам факт образования узкого канала можно понять так: при действии ультразвука, а также после про- 
боя газа и появления в нем проводимости в местах протекания тока выделяется джоулево тепло. Электропро- 
водность газа, как известно, сильно возрастает с температурой. Так, при высокой степени ионизации, когда су- 
щественны столкновения электронов с ионами, электропроводность пропорциональна 73”. При малой иониза- 
ции эта зависимость более резкая, т.к. с ростом Т быстро возрастает степень ионизации и, следовательно, появ- 
ляется тенденция к концентрации тока в сравнительно узком канале. В местах, где выше температура, прово- 
димость тока больше, значит, там течёт больший ток и выделяется больше тепла. Это приводит к еще большему 
разогреву и так далее [2]. 

Физическими процессами, определяющими ширину канала и предел концентрации тока, является отвод 
тепла из канала и расширение нагретой области под действием давления. Каналом можно считать область до 
точки, где температура и степень ионизации существенно понижается. В канале можно пренебречь инерцией 
газа, но надо учесть выделение и перенос тепла. В области ударной волны нужно учитывать инерцию, но мож- 
но пренебречь электро- и теплопроводностью. Эти две области разделены переходным слоем — «оболочкой» 
канала. В оболочке происходит нагрев и ионизация входящего в канал газа [3, 4]. 

Материалы и методы. Основными уравнениями задачи являются уравнения непрерывности, движе- 
ния и переноса энергии с учетом действия ультразвукового поля. Они имеют вид [5, 6] 

рр, О (0 
91 ди тд 














(2) 
52 Ооо 
91 9") дк 
2 2 д 
ВР” ие - Три’ ава” + (9) ео 3) 
91 2 гдг р 2 тд 


где р — плотность, / — скорость, р — давление, & — внутренняя энергия на единицу массы газа, 4 — поток 
тепла, / — плотность тока, Ё — электрическое поле, И — вклад ультразвукового поля. 
Уравнение состояния имеет вид 
р=(и, +п,)Т =РТ(2+1) т, (4) 


где т, — средняя масса атома, ип п, — числа электронов и ионов в единице объема; 2 средний заряд иона; 


е 
Й,—7П, Температура выражена в энергетических единицах. 
Будем считать, что ионизацию в канале можно подсчитать по формуле Саха. Внутренняя энергия газа 
гв, 1 (5) 
р Е 
2р т р2 (2+1Т 





где 1 — полная энергия ионизации плюс энергия диссоциации, отнесенная на один атом. Формулу (5) удобно 
применять в случае полной ионизации; при неполной ионизации с увеличением Т возрастает и величина 1. При 


этом, как следует из формулы Саха, /ИТАсоп»1, поэтому для & более удобна формула 
1 6 
В (6) 
у=Ьр 
где У— эффективный показатель адиабаты; для воздуха У=1,22. 


Электропроводность ОИ теплопроводность Х сильно ионизированного Газа равны 
с=в, (Т)Т° = 30 (=)Т*? /(4е? Влт».) 


5/2 и 
Хе (ЕТ 
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Здесь е, т — заряд и масса электрона, /, = п (эт к. / (== Али, )) у ов) — безразмерный коэффициент. 


Я 
Для 2=1 значение ©’ =1,95. Величина Хе” | ОТ ‚ согласно закону Видемана-Франца, порядка 1. Для Д=5 © 


1)=3,4 те в, 


х.()=3,9 и вы 


Результаты исследования. Предположим, что зависимость от времени радиуса канала, граница кото- 
К ь 
рого играет роль поршня, вытесняющего газ, имеет вид а(И=А ;‚ движение в области ударнои волны определя- 


ется двумя размерными параметрами А,р. 
В уравнениях (4)-(6) введем безразмерные обозначения (а, — радиус фронта волны) 








х=и/а, (0), р'(х) = р/бь» /’(х)=У/а.› '(х) = р/роа? ® 
Пренебрегая выделением и переносом тепла, запишем систему (1)-(3) в форме 
, а’ ‚„аху’ 
И + =0 
х) Ах р хах 
, , 9 
(-ри+ о. ар 0 (9) 
К ах р’ ах 





1 ар’ аху' 
211--|р’-+(Г’-х +у р" =0 
р 9 УВ хх 


С граничными условиями при х=1 
р’=(у+1)/(у-1, Г’=2у-1), р’=2(у-1) (10) 
Положение поршня определяется точкой, где Г=х Е Давление на поршень р р можно выразить через ско- 


рость поршня. 
Рь =К рой’, (0 


2,2 
где «коэффициент сопротивления» К,=9,9 (К,=р (а)а /а ) находится из численного решения системы (9). 


Будем пренебрегать излучением и примем, что 


ат (12) 
= 
Аг 
Пусть температура Т в канале много больше, чем требуется для полной ионизации, следовательно, на 


краю Т много меньше, чем в центре. Положим Т=0 при г=а. Введем безразмерные обозначения 


т? Т ии’ Г а э5Ии г (13) 
5=— <, 0(5) ==, и= ‚У + у 
а (:) ув 02а\а а 2а [ра 2\а а 


где Т, —_ температура на оси. Давление считаем постоянным по сечению канала. 











Если отвод тепла из канала осуществляется «прозрачным излучателем», то для области канала можно 
указать простое автомодельное решение: давление, температура и плотность постоянны по сечению, а скорость 
пропорциональна радиусу. Все падение температуры сосредоточено в оболочке. Там же поглощается излучение 
и происходит ионизация газа, поступающего в канал. Считая оболочку тонкой, можно получить систему урав- 
нений для основных параметров канала. В общем случае можно пользоваться для оценок этими уравнениями 
как математической моделью, описывающей, хотя и грубо, основные процессы в канале. При этом приближен- 
но учитываются действия ультразвука и теплопроводность [7]. 

Уравнения баланса энергии для канала и оболочки имеют вид 


ай’ ап а* (14) 
а ев +9, , 


15 
м ое, -- 
р] 


где М, И’ — масса и энергия газа в канале. Уравнение (15) получается интегрированием (3) по сечению канала 
(включая оболочку) без предположения о виде распределения величин по сечению. Для однородной модели по- 
ложим 
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Й=М-=, М=п-а’р. 
Уравнения (14), (15) есть следствие закона сохранения энергии. Выражения для выделения тепла О, за 
счет электрического поля и тепла О’, за счет ультразвукового поля, а также для теплоотвода излучением О, и 
теплопроводностью О’, можно взять в виде 
О =Л/пав, О, =п® (16) 
Оетао (рт Обо т (17) 


где ©) — частота ультразвуковых колебаний, "| — размерный коэффициент. 
Сравнивая (14) и (15), получим, что 


О.+0,=н (0,+0,) , (18) 
где И — коэффициент порядка 1. Если Т не зависит от #, то 
42а? ый (19) 
=9|1+(у-1)2а?| -—— 
и=у 1+(у-1) р 


(1114 -1 Е 
Рассмотрим канал в воздухе с проводимостью © =2`10° сек ; Кр=0,9; 7=1,2; 1—6 следовательно, 


Е =4,5. Для радиуса канала получим выражение 


а=0,93(1--0)% р; (20) 


Здесь, а измеряется в мм, / в КА, Ё в мксек; р/=1,29710`3 с при атмосферном давлении. Если же 
ударная волна слабая, то радиус аналогично находится из (19). 
Ниже в таблице 1 приводятся значения радиуса, рассчитанные по формуле (20), при разных значениях 
0 и  (мксек) при напряжении разряда /=30 В и емкости батареи с=0,15 мкФ, а также при индуктивности в цепи 
[=4 нГн (что овен О А/сек). 
Таблица 1 
Радиус канала при различных значениях 6 1 











т 
0 0,3 0,5 1 

0 1,62 
1 1,82 
2 1,95 
3 0,82 1,26 2,04 








Проведем оценки температуры в канале. Считаем, что Д -—1, для того же разряда, что и выше в момент 
-1 мксек при /=4 нГн имеем, что О, +О, = (1+-6)1,7.10° эрг/см сек. Если считать, что все тепло переносится 


электронной теплопроводностью, а излучением пренебречь, то получим, что 
И 


Т=4(1+6)ет. Приняв Т=4еГ найдем, что количество ионов в единице объема в этом случае будет п-9 10 , что 
по порядку соответствует экспериментальным значениям. 

Заключение. На основании построенной приближенной модели можно сделать следующие выводы о 
влиянии ультразвука на развитие искрового канала. 


1. Радиус канала увеличивается в ( 1+0) 6 раз по сравнению со случаем, когда нет ультразвука, где 0 
— отношение энергий электрического и ультразвукового полей. 

2. Температура в канале увеличивается пропорционально (1+0) в предположении, что отток тепла осу- 
ществляется электронной теплопроводностью. 

3. Уже в момент образования ударной волны наступает почти полная ионизация в канале и могут со- 
здаваться условия для двукратной ионизации при высоких энергиях ультразвука. 
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+? Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 


Верайшто Не таш зВай о? дгуег тоа$ег ыы 


У. С. Гучаши$Ку', $. $ Аззашепко””" 


+2 роп Зе Тесымса! Ошуегзну, Возюу-оп-Ооп, Ваззап Еедеганоп 


Введение. Статья посвящена анализу причин повреждений 
и износа главного вала тостера сушилки. Рассмотрены спе- 
цифика ремонтного производства и ограничения в выпол- 
нении сварочных работ, которые необходимо учитывать 
при создании технологии, обеспечивающей восстановление 
эксплуатационных свойств конструкции. 

Цель работы — создание технологии ремонта главного вала 
тостера без его разборки. Для реализации этой задачи была 
установлена расчетная ремонтная конструкция и выполне- 
ны послесварочные операции, обеспечивающие требования 
технической и нормативной документации к данной кон- 
струкции. 

Материалы и методы. В программной среде «Компас 3)» 
разработаны модели: целого вала (конструкция вала в про- 
ектном состоянии); вала, поврежденного в результате дли- 
тельной эксплуатации (более 15 лет); поврежденного вала с 
приваренной ремонтной конструкцией. Проведено числен- 
ное моделирование напряженно-деформированного состоя- 
ния (НДС) моделей. 

Результаты исследования. Разработано программное обес- 
печение для численного моделирования НДС ремонтной 
конструкции. С помощью модели получена конструкция 
ремонтного вала, в которой максимальные напряжения не 
превышают напряжений вала в проектном состоянии. Для 
исключения влияния агрессивной среды на коррозионную 
усталостную прочность вала использован изоляционный 
способ. Для присоединения ремонтной конструкции к валу 
разработана технология выполнения швов, позволяющая 
ограничить биение 0,12 мм. 

Обсуждение и заключения. Рассмотрены работающие при 
кручении валы, поврежденные в результате износа и корро- 
зии на большую глубину. Для восстановления их конструк- 
тивной прочности целесообразно использовать комплекс 
перечисленных далее методов. Конструктивный (заключа- 
ется в установке дополнительных деталей, компенсирую- 
щих недостаточную прочность, и обеспечивает уменьшение 
концентрации напряжений в наиболее нагруженных зонах). 
Технологический (уменьшает остаточные сварочные 
напряжения за счет рациональной последовательности де- 
формирования, способствующего наведению благоприят- 


* Работа выполнена по договору № 130 от 19.04.2014 г. 


* Е-пшай: [уадтизку40@тай.га, аззашетко з@ тайтл 


** ТБе гезеагсВ 15 Чопе ул ш Ше йаше оЁ Сошегасе Мо. 130 оЁ 19.04.2014. 
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ных остаточных напряжений сжатия). Изолирующий (осно- 
ван на нанесении антикоррозионных покрытий). 
Экономическая целесообразность разработанной техноло- 
гии ремонта несомненна. Ремонт вала без его разборки сто- 
ил 180 тысяч рублей, а цена нового вала — 3,8 млн рублей. 


Ключевые слова: ремонт сваркой в производственных 


условиях, ремонт вала, выбор типа ремонта, ЗО- 
моделирование, метод конечных элементов (МКЭ), допол- 
нительные детали, концентрация напряжений, технология 
ремонта, пластическое деформирование, экономическая 


целесообразность. 
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Введение. Выход из строя машин часто бывает обусловлен негативным влиянием конструктивных, 
технологических или эксплуатационных факторов. Чтобы выбрать тип ремонта, необходимо знать: 

— особенности работы конструкции, 

— условия ее эксплуатации; 

— возможные причины выхода из строя основных элементов. 

Специфика ремонтного производства состоит в том, что при выполнении работ приходится иметь дело 
с уже созданными конструкциями. Это требует выбора определенных способов сварки, расположения ремонт- 
ных сварных соединений, ограничивает доступ к местам сварки. В некоторых случаях сужаются возможности 
подогрева металла при сварке и последующей термической обработке сварного соединения [1]. Все это приво- 
дит к необходимости создания специальных технических решений в каждом конкретном случае. 

Маслоэкстракционные заводы для отгонки растворителя (бензина) из шрота (жмыха) используют чан- 
ный испаритель (тостер), снабженный греющими рубашками и трубопроводами подачи пара. Тостер, показан- 
ный на рис. 1, представляет собой колонный аппарат, собранный из чанов (1), установленных один над другим. 











































































































Рис. 1. Схема тостера сушилки 


Людмирский Ю. Г. и др. Ремонт главного вала тостера сушилки 


ГуиатизКу У. С. апа ше ойег5. Кераттз йе тат зйай о 4гуег оаяег 





Через все чаны проходит вал 2, на котором закреплены парные мешалки шрота 3, установленные над 
днищами чанов с небольшим зазором. Привод вала осуществляется через редуктор 4 от двигателя мощностью 
160 кВт и вращается с частотой 10 об/мин. Крутящий момент на валу тостера — 14500 кГ-м. 

Тостер сушилки, исследованный в рамках данной работы, в течение 15 лет эксплуатировался на масло- 
экстракционном предприятии филиала «МЭЗ Юг Руси», что привело к значительному местному износу вала. 

Проектные размеры вала: диаметр 285 мм, длина 6860 мм. Как показала проведенная диагностика, вы- 
работка вала локализовалась в его нижней части, где он контактирует со шротом (действующим как абразив). 
Температура шрота 105 °С. Место контакта обдувается влажным воздухом. В этой части вала на участке 700 мм 
его диаметр в результате коррозионно-механического износа уменьшился до 245 мм. 

Изношенная часть вала показана на рис. 2. 





Рис. 2. Износ вала после 15 лет эксплуатации 


Расчеты, выполненные согласно работе [2], показали, что полярный момент инерции вала в проектном 
состоянии составлял 6,6-10*см“, а в результате эксплуатации он уменьшился до 3,8-10*см“. Дальнейшая экс- 
плуатации вала без восстановления его несущей способности может привести к разрушению и взрыву паров 
бензина в тостере. В связи с этим оценивались различные варианты ремонта главного вала тостера. Замена не 
рассматривалась из-за высокой стоимости индивидуального изготовления вала. Кроме того, его демонтаж тре- 
бовал разборки тостера и помещения, в котором он установлен, что практически невозможно. 

В настоящее время наиболее широко применяют следующие способы ремонта изношенных валов: 

— наплавка металла на поврежденную поверхность [3] с целью восстановления размеров и несущей способно- 
сти деталей; 

— замена поврежденного участка конструкции новым; 

— установка дополнительных деталей или ремонтной конструкции [4]. 

Первый вариант ремонта был отклонен, т. к. в этом случае требовалось наплавить на вертикально рас- 
положенный вал большой объем металла (более 140 кг), что неизбежно привело бы к появлению дефектов и к 
большим деформациям. 

Второй способ трудноосуществим. При его реализации возникнут сложности с вырезкой поврежденной 
части вала, разделкой кромок под сварку, сборкой и сваркой. В результате появятся недопустимые сварочные 
деформации. В этом случае даже при удачной замене поврежденного участка весьма вероятно нарушение экс- 
плуатационных свойств тостера. 

Третий способ ремонта обычно рекомендуют, если повреждение поверхности имеет большие разме- 
ры [4, 5]. Таким образом, было решено ремонт главного вала тостера выполнять путем установки ремонтной 
конструкции (рис. 3). 
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Рис 3. Схема ремонтной конструкции 


Конструкция вала 1 состоит из двух полуколец 2 и двух полумуфт 3. Размеры муфты и кольца назначи- 
ли, руководствуясь прочностными и технологическими соображениями. Длина муфты, равная 800 мм, превы- 
шала длину поврежденного участка вала на 100 мм. Внутренний диаметр муфты равен 285 мм, наружный — 
325 мм. Для восстановления статической прочности вала в его изношенной части требовалось приварить эле- 
мент, полярный момент инерции которого восполнил бы его потерю при износе (6,6-10“см“- 
3,8-10* см“ = 2,8-10* см“). При вышеназванных диаметрах полярный момент инерции полумуфты равен 
3,1.10*см“. Таким образом, размеры ремонтной конструкции, приваренной к валу, не только восстанавливают 
расчетный момент инерции вала, но и увеличивают его на 15 %, что обеспечивает статическую прочность вала. 

В работах [6, 7, 8] показано, что долговечность сварных конструкций, работающих при повторно- 
статических нагрузках во влажной среде, во многом зависит от напряженно-деформированного состояния 
(НДС) в зоне перехода от металла шва к основному металлу. Для оценки влияния геометрии сварных швов на 
НДС использовали метод конечных элементов (МКЭ). Проектируя ремонтную конструкцию, стремились к то- 
му, чтобы напряжения в районе сварных швов не превышали напряжений в районе шпоночных пазов, имею- 
щихся на валу для передачи крутящего момента на мешалки шрота. 

Методика расчета НДС представлена в работе [9]. В программной среде «Компас 3)» были разработа- 
ны ЗО-модели целого вала с первоначальным шпоночным пазом, поврежденного коррозией вала и поврежден- 
ного вала с учетом установки ремонтной конструкции. Для расчета использовали программный продукт 
Апбуб [10, 11]. З)-модели импортировались в Аизу5 Иот’Ьепсй. С помощью этой программной оболочки вы- 
полнены предварительная обработка геометрии и подбор материала, соответствующего материалу вала — 
ТУ7ГТС (табл. 1). 

Таблица 1 


Параметры выбранного материала 














Свойство Единица измерения Значение 

Плотность кг/м? 7850 
Предел текучести МПа 250 
Предел прочности МПа 460 

















Далее был проведен расчет, который включает в себя создание конечноэлементной сетки, а также уста- 
новку граничных условий путем приложения крутящего момента со стороны редуктора и жесткого закрепления 
вала с обратной стороны. 

На рис. 4 представлены результаты рассчитываемых конструкций, а в таблице 2 — максимальные 
напряжения на валу: в местах расположения шпоночных пазов (рис. 4, а), в месте износа вала (рис. 4, 6), в 


наиболее нагруженном месте ремонтной конструкции (рис. 4, в). 


Людмирский Ю. Г. и др. Ремонт главного вала тостера сушилки 
ТуиатизКу У. С. апа йе ойегу. Кератиз йе тат 5йай оГагуег поаяег 





а) Ь) с) 


Рис. 4. Модели расчета НДС в местах расположения: шпоночных пазов (а); износа вала (6); ремонтной конструкции (в) 


Таблица 2 


Максимальные расчетные напряжения 














Место расчета напряжений на валу Максимальное напряжение, МПа 
В расположении шпоночных пазов 74 
В месте износа вала 124 
В ремонтной конструкции 75 














Видно, что напряжения в разработанной ремонтной конструкции практически не превышают проект- 
ных напряжений в местах расположения шпоночных пазов. 

НДС зависит не только от рабочих напряжений, но и от величины и знака остаточных напряжений в 
зоне сопряжения металла шва с основным металлом. Поэтому в технологии ремонта использовали метод пла- 
стического деформирования [12]. Он позволяет уменьшить остаточные напряжения растяжения или даже пере- 
вести их в благоприятные напряжения сжатия. 

Перед началом ремонтных работ была проведена дегазация тостера и цеха, в котором он находится. 
Поврежденную часть вала и места сварки зачистили до металлического блеска вращающимися щетками. 

Сборку ремонтной конструкции начали с установки двух полуколец 2 (рис. 3), их стяжки и прихватки. 
Для сварки полуколец использовали электроды УОНИ 13/55 [13] на режимах, указанных в табл. 3. 




















Таблица 3 
Режимы ручной дуговой сварки покрытыми электродами 
Марка электродов Рекомендованная сила тока при РДС*, А, в положении: 
Диаметр, мм 
нижнем вертикальном горизонтальном и потолочном 

УОНИ 13/45, 3,0 80-100 60-80 70-90 
УОНИ 13/55 4,0 130—160 100-130 120-140 

5,0 170-200 140—160 150—170 

















*РДС — режим дуговой сварки. 











После этого произвели зачистку швов заподлицо с основным металлом. 
Затем кольцо установили перпендикулярно валу, как показано на рис. 3, приварили его и выполнили 


сварку кольцевых швов, присоединяющих кольцо 2к валу, по схеме, показанной на рис 5, а. 
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Рис. 5. Схема выполнения: кольцевых швов (а), продольных швов (6) 


Сварку выполняли одновременно два сварщика. Порядок наплавки валиков (1-18) показан на рис. 3. 

На приваренное кольцо установили две полумуфты — так, что расстояние от продольных швов кольца 
до продольных швов муфты составляло 120-130 мм. Полумуфты обжали вокруг вала с помощью двух центра- 
торов и выполнили в каждом стыке по три прихватки длиной 50-60 мм. Для уменьшения вероятности появле- 
ния дефектов при сварке начало и конец каждой прихватки зачистили абразивным инструментом. После этого 
верхний центратор удалили и приступили к сварке продольных швов муфты. Сварку выполняли обратносту- 
пенчатым способом [14] снизу вверх — в порядке, показанном на рис. 5, 6. После выполнения 70 % длины про- 
дольного шва удалили нижний центратор и закончили выполнение корневого шва. Заполняющие и облицовоч- 
ные валики выполняли в том же порядке. Каждый последующий слой смещали относительно предыдущего на 
25-30 мм. Силу тока, рекомендованную в табл. 2 для сварки вертикальных швов, уменьшили на 10 %. 

Кольцевые швы, присоединяющие муфту к валу и кольцу, выполняли одновременно два сварщика, 
расположенные в диаметрально противоположных местах относительно вала (рис. 5, а). Вначале были наплав- 
лены валики 19—27, а затем 28—37 (рис. 3). После выполнения кольцевых швов произвели их зачистку и пласти- 
ческое деформирование по линии сопряжения металла шва с основным металлом вала. Качество сварных со- 
единений проверяли визуальным контролем послойно. 

Радиальное биение вала измерялось выше ремонтной конструкции на 100 мм индикаторами часового 
типа в двух взаимно перпендикулярных направлениях. Его значение не превысило 0,12 мм. 

С целью защиты от коррозии муфты и вала пространство между ними заполнили подсолнечным мас- 
лом [15], а наружную поверхность вала изолировали эпоксидным клеем марки ЭДП по ТУ 07510508.90-94. 

Целостность отремонтированного вала контролируется ежеквартально. Он безотказно работает в тече- 
ние четырех лет при полном соответствии функциональных показателей требованиям технической и норматив- 
ной документации. 

Экономическая целесообразность разработанной технологии несомненна. Ремонт обошелся в 
180 тысяч рублей, а цена нового вала — 3,8 млн рублей. 

Выводы. Для восстановления конструктивной прочности валов, работающих при кручении, целесооб- 
разно использовать комплекс перечисленных далее методов. Конструктивный предполагает установку допол- 
нительных деталей, компенсирующих недостаточную прочность, и обеспечивает уменьшение концентрации 
напряжений в наиболее нагруженных зонах. Технологический уменьшает остаточные сварочные напряжения за 
счет рациональной последовательности деформирования, способствующего наведению благоприятных оста- 
точных напряжений сжатия. Изолирующий заключается в нанесении антикоррозионных покрытий. 
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Методика функциональной унификации адаптивного модуля гидропривода с функцией 
стабилизации нагрузки на рабочем органе мобильных машин” 


В.А. Першин', Т. А. Хиникадзе?“" 


1? Институт сферы обслуживания и предпринимательства (филиал) г. Шахты Донского государственного технического университета, 


Российская Федерация 
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12 шение оЁ Зегутсе ап Визтезз (О$ТО БгапсВ), ЗВаКму, Визчап Еедеганоп 


Введение. Исследованы вопросы — функциональной 
унификации адаптивного модуля гидропривода. Впервые 
рассмотрены принципы  самоадаптации с учетом 
внутренних адаптивных связей управления нагрузкой и 
согласованности движений на рабочем органе мобильных 
машин. 

Цель работы — 
функциональной унификации адаптивного модуля 


гидропривода. Для достижения поставленной цели решен 


создание И анализ мегодики 


ряд задач. Обоснован выбор технических устройств — 
унифицируемых адаптивных модулей гидропривода 
мобильных машин. Предложены методика и показатели 
функциональной унификации Рассмотрены 
внутренние связи: прямая положительная и обратная 


модуля. 
отрицательная. Показано их влияние на свойство 
функциональной унификации адаптивного модуля. 
Материалы и методы. Для синтеза и анализа показателей 
функциональной унификации адаптивного модуля принят 
метод подобия функционирования технических систем. 
Результаты исследования. Разработана методика 
проведения конструктивно-функциональной унификации 
модулей с самоадаптацией. Показаны возможные варианты 
модификации унифицируемых модулей и соответствующие 
им сочетания гидродвигателей, регулирующей аппаратуры 
и математических моделей адаптивных — связей. 
Предложены критерии и индикаторы подобия. Выполнен 
анализ функциональной унификации внутренних 
адаптивных связей модуля с разными типами 
гидродвигателей и элементов дросселирования жидкости в 
гидросистеме. Сформулированы рекомендации проведения 
функциональной унификации при типизации и 
эксплуатации (наладке) адаптивного модуля. 

Обсуждение и заключения. Методика рекомендована для 
функциональной унификации гидравлического модуля с 
самоадаптацией. Она может использоваться при разработке 
типоразмеров модуля и в процессе его эксплуатации в 
качестве — автономного 
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Введение. Многофункциональные и комбинированные мобильные машины активно используются в 
сложных условиях: при выполнении дорожно-строительных и сельскохозяйственных работ. Как правило, их 
эксплуатация предполагает ручное или жестко запрограммированное управление. Таким образом, актуально 
создание простых функционально унифицированных приводов (модулей) мобильных машин, способных 
адаптироваться к переменным свойствам обрабатываемой среды [1]. 

Следует отметить, что вопросы унификации связаны с решением ряда научно-практических задач. 
Одна из них —Й проблема типоразмеров функционально унифицированных изделий. Для ее решения 
используется соответствующий аппарат фундаментальной [2] и прикладной [3, 4, 5] математики. Проводятся 
исследования адаптивных систем различной структуры и назначения [6], гидроприводов устройств с адаптивно 
согласованными движениями рабочего органа, воспринимающего стохастические переменные нагрузки [7, 8, 9, 
10]. В рамках данной работы рассматривается принцип функциональной унификации в гидроприводе 
модульного типа с функцией самоадаптации по нагрузке. Исследование и практическое испытание такого 
модуля проводятся впервые. 

Материалы и методы. В качестве базовых объектов исследования приняты устройство для бурения 
породы с переменными свойствами [11] и устройство для обработки криволинейных поверхностей [12] (рис. 1, 
2). 





Рис. 1. Принципиальная схема устройства для бурения породы с переменными свойствами [3] 
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Рис. 2. Принципиальная схема устройства для обработки криволинейных поверхностей [2] 


На представленных схемах пронумерованы следующие элементы (для рис. 2 обозначения даны в скобках): 
1 (10) — гидронасос постоянной производительности; 2 (11) — предохранительный клапан; 3 (12) — 
фильтр; 4, 7, 8 — гидрораспределители; 5 (4) — дроссель (регулятор режимов); 6 (1) — гидромотор главного 
движения; 9 — регулятор потока; 10 — подпорный регулируемый дроссель; 11 — гидроцилиндр подачи; 12, 13, 
14 — манометры; 15 — бак; (2) — рабочий орган; (3) — гидросистема; (5) — реверсируемый гидромотор 
подачи; (6) — вал гидромотора подачи; (7) — механизм подачи; (8, 9) — регулируемые дроссели. 

Элементы гидросистемы этих устройств (в том числе их модификации [13]) структурно объединены по 
дифференциальной схеме. 

Между сливными линиями регулятора режимов и гидромотора главного движения дополнительно 
включены гидроцилиндр (рис. 1) или гидромотор (рис. 2). Таким образом устанавливаются две связи: обратная 
отрицательная и прямая положительная. 

На сливной линии гидромотора главного движения рассматриваемого устройства для бурения 
находится регулятор потока (рис. 1). На сливной линии гидромотора главного движения устройства для 
обработки криволинейных поверхностей находится регулируемый дроссель (рис.2). В этом состоит 
принципиальное отличие данных устройств. 

Такие вариации унифицированных модулей позволяют обеспечить сочетание движений рабочего 
органа мобильной машины: вращение — подача; вращение — вращение; подача — подача [14]. 

Конструктивные отличия влияют на характеристики прямой положительной и обратной отрицательной 


связей устройств, одновременно унифицированных и модульных (табл. 1, 2). 

















Таблица 1 
Функциональные зависимости адаптивных связей модуля с регулятором потока 
В Функции адаптивных связей модуля 
Типы гидродвигателей модуля 
Обратная отрицательная Прямая положительная 
Гидромотор главного движения 
«ши = АМо — (ВМед + СМеди) гп = К Эд 
— гидромотор подачи 
Гидромотор главного движения 
«ши = АМо > (ВМед + СМеди) Уцп — К20гд 
— гидроцилиндр подачи 
Гидроцилиндр главного движения 
ЕН й Уца = АМо — (ВМад + СМаш) Уци = 3 Ул 
— гидроцилиндр подачи 
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Таблица 2 
Функциональные зависимости адаптивных связей модуля с регулируемым дросселем 
Функции адаптивных связей модуля 
Типы гидродвигателей модуля 
Обратная отрицательная Прямая положительная 

Гидромотор главного движения ов |2АР;ьв 
© = АМо — (ВМад + СМедп) т = Кд — И АРВ. |7 АР, 

— гидромотор подачи Ты р 
Гидромотор главного движения ов |2АРьв 
и = АМо — (ВМад + СЕ) И п = 20 д — И АВВ: | ееАРЬ: 

— гидроцилиндр подачи а р 
Гидроцилиндр главного движения Ев |2АРльв 
7 = АМ — (ВЕд +СЕш) ип = ЕзИ д — Арес АРВ 

— гидроцилиндр подачи Ра р 














Отметим, что соотношение скоростей выходных элементов гидродвигателей (см. табл. 2) не зависит от 
расхода через регулятор потока, так как этот показатель не зависит от перепада давления на аппарате, т. е. 
Ор = с0п54. 

Адаптивные связи играют особую роль при действии рабочих органов машин в режимах с требуемыми 
показателями скорости выполнения операций и устойчивости [14]. Рассмотренные выше конструктивно- 
функциональные особенности устройств (наличие регулируемого дросселя или регулятора потока) позволяют 
реализовать такие связи. 

В таблице 1 и 2 приняты условные обозначения: 17 — скорость штока; Ли — площадь поперечного 
сечения гидроцилиндра подачи; о;„ —Й угловая скорость гидромотора главного движения; А, ВС — 
постоянные коэффициенты, зависящие от моментов инерции, кинематических параметров и КПД узлов 
привода главного движения и привода подачи модуля; Мо — номинальный и расчетный суммарный момент 
сопротивления, соответствующий процессу обработки без адаптации; Дрв — площадь проходного сечения 
дросселя Дрз; Мал, Ми — действительные моменты сопротивления, воспринимаемые соответственно валом 
привода главного движения и валом подачи; Ем,Е!п —Й действительные усилия сопротивления, 
воспринимаемые соответственно штоком гидроцилиндров главного движения и цилиндра подачи; ДАР — 
перепад давления на дросселе Дрз; р — плотность рабочей среды в гидросистеме (масло); Кл, К2, Кз 
коэффициенты перевода конструктивных особенностей и размерностей скоростей движений. 

Отрицательная обратная связь позволяет в процессе работы модуля стабилизировать действительный 
суммарный момент (усилие) сопротивления на рабочем органе. В этом случае автоматически сравниваются: 

— действительный момент сопротивления со значением, заданным регулятором режимов; 
— изменения скоростей подачи и главного движения со знаком, обратным сигналу рассогласования. 

Ниже приводятся соответствующие уравнения обратной связи по ускорениям вала гидромотора 
подачи. 

При наличии регулятора потока: 
1 +4 (ВАМОд+ САМО) 
| (ВАМ. + САМ0 1)05 |" (1) 





При наличии регулируемого дросселя: 


а 0 0 Г Г 
1 +а(ВАМад+ САМО п) а |2(АРдрз+ 6(АРдрз) 
Фи =- [(А Г] 0 05 + Мдрздрз — )]. (2) 
Ксгп (ВАМегд+ САМсди) ЧЕ р 


Методика проведения функциональной унификации определяется ее целями и особенностями, а также 
этапом жизненного цикла модуля. Вместе с тем следует отметить и общие шаги ее реализации: 





— уточнение функционального назначения и условий функционирования унифицируемого модуля; 

— предварительный сбор данных в рамках конструктивной и функциональной унификации 
гидроаппаратов типового и унифицируемого модуля (в том числе условий однозначности; номинальных и 
граничных значений параметров функционирования, особенно перепадов давлений на гидроаппаратах в 
условиях динамического равновесия гидросистемы); 

— уточнение вида адаптивных связей по конструктивным особенностям модулей (в соответствии с 
табл. 1и2). 
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При соблюдении отмеченных общих положений функциональная унификация адаптивных связей 
выполняется в определенной последовательности. 

— Расчетом или по техническому заданию, паспорту уточняют значения рабочих (или максимальных) 
нагрузок на рабочем органе мобильной машины с унифицируемым модулем. 

— В соотношении с типовым модулем определяют функциональные и конструктивные значения 
параметров элементов гидросистемы унифицируемого модуля. При этом используются номограммы, графики 
или расчеты (например, методом подобия функционирования) [15]. 

Результаты исследования. При выполнении унификации используется принцип суперпозиции 
функций. 

Кратко рассмотрим пример методики проведения функциональной унификации гидродвигателей 
гидросистемы модуля. 

Выходными характеристиками гидромотора и гидроцилиндра являются угловая скорость (частота 
вращения) и момент на валу гидромотора, линейная скорость и усилие на штоке гидроцилиндра. 
Функциональные зависимости, определяющие эти характеристики, имеют вид [16]: 


а) для гидромотора 


М. м Л . я Ч, В, У, а : ВР (ед я и й р)), 


(3) 
Пи = рр (О, Ч, о, В, О АЕ, я. р)); 
6) для гидроцилиндра 
7, = (0,4, „6,р,АР,,Е,И,), а) 
т — ГАР, 2 АР, ы Ту ы т ы и, бы ы В). 
Здесь и„„ — частота вращения вала; 4 — объем масла; О — расход (подача) масла; Е — модуль упругости; и — 


динамическая вязкость масла; р — плотность масла; Р„, — давление на входе; АР,, — перепад давления на 
элементе утечек; Оу, — величина утечек; /„ — площадь щели, через которую проходит утечка; / — момент 
инерции; В — коэффициент сжатия масла; Ри, Е», — приведенные силы страгивания и внутреннего трения в 
гидроцилиндре; / — радиус вращения масс; &, у — конструктивные коэффициенты; М, — момент на валу; М.р 
— момент трения; АР, — перепад давления на гидроцилиндре; Р„ — усилие на штоке; АРи — разность усилий 
на поршне. 

При проведении конструктивно-функциональной унификации используют частные критерии подобия 
унификации [17], полученные методом анализа размерностей [15]. 

Для гидромотора критерии подобия унификации имеют вид: 





пи Ч м5. Е. 
т, о в, й- Е - ДР’ =ВР., (5) 
Мы’ САР ЛР. 
о а о 
а Оси? _ р а”” 
(9. 





с Ра ре Р пы” ° 

Критерии подобия унификации имеют определенный функциональный или конструктивный смысл, а 
именно: Л, —Й критерий расходно-кинематических параметров; л, —Й критерий, характеризующий 
соотношение конструктивных и расходных (функциональных характеристик); п; — критерий, 
характеризующий номинальную жесткость рабочей жидкости; л, — критерий теоретического (индикаторного) 
момента на валу насоса; л, — критерий упругой деформации рабочей жидкости; п„„ — критерий потерь на 
трение в сопряжениях; л, — критерий потерь на инерционное, гидромеханическое сопротивление; ло — 
критерий потерь на гидродинамическое сопротивление рабочей жидкости; л,, —Й критерий потерь на 


преодоление утечек рабочей жидкости; л„ — критерий потерь на внутреннее (вязкое) сопротивление рабочей 


среды. 
Таким образом, каждый критерий отражает сущность одного из свойств функциональной унификации 
модуля, а любая система критериев отражает это свойство обобщенно. 
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Обобщенные критерии подобия унификации, полученные путем 
объединения частных критериев, представляются в виде: 


ПРП, Ра°,ЛУ „Ор | 


тр 


2 
-_ Ч/»ЬЕАЁ, 
М, = ЗИ, = УЕ. 
ППП пот, — ВНЕ) 4 УРРО, 
Соответствующие частные индикаторы подобия функциональной унификации С, получают из 
уравнений (5) для унифицируемого и типового модулей [15]. При этом должны быть равны соответственные 
числа подобия лу; и ли; для типового и унифицируемого модулей: 
пу = пи = ет. (7) 
В рассматриваемом примере частные индикаторы подобия функциональной унификации имеют вид: 





(6) 


С. С СЕ См:рб 3 
1= Сл; 154 Са; 15; 1 = Саб; 1 ЕН (8) 
Са ит Са Л ц САР 2% Рвх/ 
ыы Смгы° . 1 АР. 1 — С9Срозз . — С9 САРгм 
= — — 0,33 Й = 2 ь 
СРьхС]з 2 Вх" Ри””?а Са” Ср 


Так же, как и критерии подобия функциональной унификации, частные индикаторы можно объединять 
в комплексы, получая эффект суперпозиции. 

Изменение масштаба параметров С; унифицируемого модуля в индикаторах (8) по отношению к 
аналогичным параметрам индикаторов типового модуля должно соответствовать условию равенства 
единице [15]. Например, для выходной характеристики ДР.„ по индикатору 


= Со. САРтм 
1-6 ©) 


можно исследовать функциональную унификацию по перепаду давления, характерному (рабочему) объему и 
расходу гидромотора, плотности масла. 

Искомый типоразмер одного из параметров унифицируемого модуля определяется подстановкой в (5) и 
(8) известных значений параметров гидромотора типового модуля и априори известных (установленных 
техническим заданием и т. п.) значений параметров унифицируемого модуля. 

Особенности процедуры исследования функциональной унификации зависят от поставленной задачи. 
При проведении унификации модуля нового типоразмера сначала выполняют конструктивно-функциональную 
унификацию, выбор типоразмеров элементов гидросистемы и затем проводят проверку их функционального 
соответствия внутренним адаптивным связям. Для решения первой задачи используются частные критерии и 
индикаторы подобия функциональной унификации. Для решения второй задачи используются обобщенные 
критерии, индикаторы, а также уравнения прямой и обратной адаптивных связей. 

Другими задачами унификации адаптивного модуля являются: исследование причин нарушения 
функциональной унификации в процессе эксплуатации одного и того же адаптивного модуля, а также 
исследование эффективности функциональной унификации при модернизации адаптивного модуля 
гидропривода. При этом также может быть использован метод подобия функционирования технических 
систем [18], [19]. 

Обсуждение и заключения. Обоснована актуальность разработки и применения в мобильных 
машинах унифицированного гидропривода модульного типа. Данная конструкция позволяет адаптировать 
силовые и кинематические параметры при воздействии на рабочий орган переменных нагрузок. Методика 
проведения конструктивно-функциональной унификации модулей с самоадаптацией предполагает 
использование математических моделей адаптивных связей. Показаны возможные варианты модификации 
унифицируемых модулей с самоадаптацией и соответствующие им сочетания гидродвигателей, регулирующей 
аппаратуры и математических моделей адаптивных связей. Приводится пример методики проведения 
функциональной унификации с использованием метода подобия функционирования технических систем. 
Методика рекомендована для исследования функциональной унификации при разработке типоразмеров и в 
процессе эксплуатации предложенного типа модуля. 
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Прогнозирование уровня подземных вод месторождения цементного сырья на основе 


динамических окрестностных моделей” 


И.А. Седых"”" 


' Липецкий государственный технический университет, г. Липецк, Российская Федерация 
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1. А. бедукв"” 
'ТдрезК Зе Тесыиса! ОшуегзИу, ГлрезК, Виззап Еедегайот 


Введение. Статья посвящена разработке математической 
модели уровня подземных вод месторождения цементного 
сырья, расположенного в задонско-елецком водоносном 
горизонте, являющимся основным источником 
хозяйственно-питьевого водоснабжения города Липецка. 
Поэтому на стадии разработки месторождения 
необходимо проводить постоянный мониторинг и иметь 
возможность прогнозирования уровня подземных вод. 
Цель работы — идентификация и исследование 
динамической окрестностной модели с переменными 
иерархическими окрестностями уровня подземных вод, 
позволяющей с достаточной точностью прогнозировать 
значение уровня вод в обследуемых скважинах. 
Материалы и методы. Приведено определение 
динамической окрестностной модели с переменными 
отличающейся 
изменяющимися во времени двухуровневыми 


окрестностными связями между узлами первого и второго 


иерархическими окрестностями, 


уровня. В каждый следующий дискретный момент 
времени узлы окрестностной модели меняют свое 
состояние под воздействием текущих управлений и 
состояний узлов, входящих в их окрестности. В качестве 
частного случая рассмотрена модель с линейными 
функциями пересчета состояний. Параметрическая 
идентификация динамической окрестностной модели 
заключается в нахождении параметров системы для 
каждого узла второго уровня и основана на методе 
наименьших квадратов. 

Результаты — исследования. Разработана линейная 
динамическая окрестностная модель с переменными 
иерархическими окрестностями для прогнозирования 
уровня подземных вод месторождения цементного сырья, 
расположенного в задонско-елецком водоносном горизонте. 
Для параметрической идентификации и моделирования 
функционирования рассматриваемой динамической 
окрестностной модели разработано программное 
обеспечение на языке С++, позволяющее для заданной 
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структуры находить параметры функций пересчета 
состояний узлов, а также прогнозировать поведение модели 
в процессе функционирования. Приведена иерархическая 
структура и проведена параметрическая идентификация 
линейной динамической окрестностной модели уровня 
подземных вод. После выполнения параметрической 
идентификации на обучающей 
математическая модель проверена на контрольной выборке. 
Обсуждение и средние 
относительные ошибки идентификации и прогноза 


выборке данных 


заключения. Полученные 
свидетельствуют об адекватности разработанной модели и 
позволяют рекомендовать ее для прогнозирования уровня 
подземных вод месторождения цементного сырья. 


Ключевые 
месторождение 
окрестностная модель с переменными иерархическими 


слова: уровень 
цементного 


подземных вод, 
сырья, динамическая 
окрестностями, параметрическая идентификация. 
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Введение. Задонско-елецкий горизонт является основным продуктивным комплексом известняков, 
пригодным для цементной и металлургической промышленности [1]. Известняки светлого и светло-серого 
цвета, с желтоватым оттенком, средней крепости, трещиноватые, мелкозернистые. Часто встречаются крепкие 
известняки, представленные окремнёнными разностями. Известняки нередко пористые и кавернозные. 
Глинистые породы, пригодные для цементного производства, залегают среди четвертичных отложений, 
покрывающих известняки. 

Вскрышные породы представлены плодородным слоем почвы, суглинками, некондиционной глиной, 
песком и слоем щебенки в верхней разрушенной части известняков. 

Гидрогеологические условия месторождения Сокольско-Ситовское месторождение 
приурочено в геоморфологическом отношении к присклоновой части долины реки Воронеж, чем обусловлена 
гидрогеологическая обстановка в районе месторождения [2]. 


простые. 


Питание водоносного горизонта осуществляется за счет инфильтрации атмосферных осадков и за счет 
поглощения паводкового стока по долинам балок и оврагов. Разгрузка подземных вод происходит в реке 
Воронеж. 

Современное положение уровня подземных вод изучено по режимным наблюдениям в 7 скважинах, 
равномерно расположенных на площади месторождения. Уровень подземных вод в меженный период меняется 
от 102,5 м до 109,7 м. 

В настоящее время задонско-елецкий водоносный горизонт является основным источником 
хозяйственно-питьевого водоснабжения города Липецка. В связи с этим разработка известняков может 
производиться только в необводненной их части с оставлением предохранительного целика не менее 2,0 м, что 
рекомендовано санитарной службой Липецкой области. 

Водоприток в разрабатываемый карьер месторождения возможен только за счет атмосферных осадков. 
В этой связи на стадии разработки месторождения необходимо проводить постоянный мониторинг состояния 
подземных вод задонско-елецкого водоносного горизонта и иметь возможность прогнозирования уровня 
подземных вод. 

В работе для моделирования уровня подземных вод рассматриваемого месторождения применяются 
динамические окрестностные модели с переменными иерархическими окрестностями, позволяющие 
моделировать сложные пространственно-распределенные процессы и объекты, изменяющие свое состояние во 
времени [3-8]. 

Окрестностные модели впервые предложены в конце 90-х гг. ХХ века [9]. В [9-12] приведены 
основные определения и алгоритмы теории окрестностного моделирования. Понятия «окрестности», 
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«окрестностные связи» рассматривается и в [13-14]. В них используются агенты, перемещающиеся по 
окрестностям и взаимодействующие друг с другом по определенным правилам. Сегодня теория окрестностного 
моделирования активно развивается. Появились динамические [6-8], недетерминированные модели [11], 
модели с переменными окрестностями [15]. 

Материалы и методы. В данной работе используются динамические окрестностные модели, 
применяющиеся для прогнозирования состояний пространственно распределенных систем. В них каждый узел 
которых является самостоятельным объектом, функционирующим во времени и связанным окрестностными 
связями с другими объектами системы. Это отличает рассматриваемые модели от широко применяющихся 
сегодня нейронных сетей, которые можно использовать для моделирования процесса функционирования 
каждого объекта или узла в отдельности. Методика окрестностного моделирования предназначена для 
одновременного совместного моделирования и прогнозирования поведения всех элементов распределенной 
системы. 

Далее рассмотрим динамические окрестностные модели «вход-состояние» с переменными 
иерархическими окрестностями. Они отличаются изменяющимися во времени двухуровневыми 
окрестностными связями между узлами первого и второго уровня и линейными функциями пересчета 
состояний. Они позволяют, по сравнению с известными одноуровневыми окрестностными моделями, 
выполнять прогнозирование с более высокой точностью. 

Динамическую окрестностную модель «вход-состояние» с переменными [15] иерархическими 


окрестностями можно задать набором №. = (№ , Г, Ч,Х [0],#) ‚ где: 


1) №=(4,0.,0,,О.„) —щ двухуровневая структура окрестностной модели; даный, — 
множество узлов первого уровня; О, и О, — окрестности связей узлов по состояниям и по управлениям 
соответственно; О„ — иерархические окрестностные связи между узлами. Для каждого узла аеА 


п 
определена своя окрестность по состояниям О [а,]<А и управлениям О[а]<А; О, =чО, [а, ] Ь 
а Е 


п 
О, = ЗО, [а,]. 
Некоторым узлам первого уровня а, е А поставлено в соответствие множество узлов второго уровня 
О. „[а,= {а!,...а*}. Все узлы а, Е О[а,], а, е ОДа, | ‚ оказывают влияние на узлы второго уровня а?’ ЕО. [а]. 
В каждый момент времени определен только единственный активный узел Б=1.....с такой, что 
а’ЕО.„[а,,. 


Узел а’ является активным в момент времени , если для него выполняется заданное условие 


активации /*[!] = ‘тае. Для всех узлов второго уровня а’ О.[а’]=О[а, 1; О[а’]=ОДа,]. 


п 


2) Хе =” — блочный вектор состояний окрестностной модели в текущий момент времени, каждый 
блок которого Х[а,| = Хе А” — вектор состояний в узле а, системы, #=1..... п. 

3) ГЕК" — — блочный вектор управлений в текущий момент времени, каждый блок которого 
Иа]=ИШЙе А” — вектор управлений в узле а, системы, 1=1,...и. 

4) С:Хо хИ, >Х — вектор-функция пересчета состояний окрестностной модели, где Хо, 
множество состояний узлов первого уровня, входящих в окрестность О,; 7, — множество управлений узлов 


первого уровня, входящих в окрестность О,. 
Для узлов первого уровня а, е А функция С. будет иметь вид: 


хи+ьй=ви=У хи+ь]=У 6, (1) 


= 


где 6’ Ао ыН® Г а ХР] — функция пересчета состояний для узла второго уровня а’ : 


Для каждого узла второго уровня а’ еО..[а,,й] функция (6 в текущий момент времени в линейном 


случае имеет вид: 
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ХИ+ЬР]= У &Р,ЛХЬЛ+ У ГР, МИЬ+ 8, (2) 

а, ЕО, [а ] а =О,[аь ] 
где аа, е А (К № ый — Узлы первого уровня модели; ХЬ?]е А” — состояние в узле а’ в момент 
времени # ИЬР]еЕ К” — управление в узле а’ в момент времени & 8’, Ле К”””, э[Р,Ке В" , 


8'[Р]е А” — матрицы-параметры модели. 
Ь Ь > 
Для каждого узла второго уровня а’ ®О.,[а,,Ё] функция С) в текущий момент времени равна нулю, то 


есть Х[+ЬР]= [И =0. 
› Р 

5) [0] ЕР" ° — начальное состояние модели. 

6) { — текущий дискретный момент времени функционирования модели. 

Структура окрестностной модели может быть представлена в виде одного двухуровневого графа 
структуры окрестностной модели — ориентированного графа с двумя видами дуг: по состояниям и 
управляющим воздействиям — или двух ориентированных графов — внешней и внутренней структур. 

Параметрическая идентификация динамической окрестностной модели [16-17] заключается в 
нахождении параметров системы для каждого узла второго уровня и основана на методе наименьших 
квадратов: 


Е=У|ХИ+ЬЙ- 6 па. 


Результаты исследования. Рассмотрим линейную динамическую окрестностную модель уровня 
подземных вод месторождения цементного сырья. Граф внешней структуры модели представлен на рис. 1. Узел 
а’ — внешняя среда. В настоящее время для целей мониторинга уровня подземных вод используются 7 
скважин, которые на рис. 1 соответствуют узлам а,-а,. Управляющие воздействия узлов И[Ё,Й Е Е® состоят 


из количества атмосферных осадков в мм и среднесуточной температуры воздуха в °С за последние 5 суток 
перед измерением уровня воды в скважинах в момент времени Е, 1=1,...,8. Состояния узлов Х|ЬЙеЕ ВР равны 


уровню воды в скважине Г в момент времени #, 1=2,,...,8. 





Рис. 1. Граф внешней структуры окрестностной модели 


Узлы первого уровня окрестностной модели а, (1=2,...,8 ) являются иерархическими: О.„[а,]= {а',а’} 


и соответствуют положительной и отрицательной среднесуточной температуре воздуха. Их структура 
представлена на рис. 2. 
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7Й 





Рис. 2. Граф внутренней структуры узлов первого уровня 


Ь КИ 
Система (2) для каждого узла второго уровня а’ окрестностной модели уровня подземных вод будет 
иметь вид: 
.Б БГ; :6 Ы) БГ; .6 БГ +В 
А а + 8, (3) 
реа Е 0. ЬГ: Е $ = 
где в [1ПеЮ; #[1ПеЕК”; & ПЕК; 1=2.....8; Б=12. 

Для параметрической идентификации и моделирования функционирования рассматриваемой 
динамической окрестностной модели разработана программа на языке С++. Начальными данными для 
идентификации являются структура и обучающая выборка. Программа позволяет находить параметры функций 
пересчета состояний узлов, а также прогнозировать поведение математической модели в процессе 


функционирования. 
После проведения параметрической идентификации на обучающей выборке данных полученная модель 


была проверена на контрольной выборке. Нормированные исходные и модельные значения уровня подземных 
вод для узла а, для обучающей и контрольной выборок приведены соответственно на рис. 3 ирис. 4. 
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Рис. 3. Нормированные исходные и модельные значения уровня подземных вод для узла а, 
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Рис. 4. Нормированные исходные и предсказанные значения уровня подземных вод для узла а, 
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Седых И. А. Прогнозирование уровня подземных вод месторождения цементного сырья на основе динамических окрестностных моделей 
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Средняя относительная ошибка идентификации (прогноза) окрестностной модели вычисляется по 
формуле: 





ма ХИ-+ЬЙ-Х, + 
А= 1 уз = [+ ‚Й 1 + 5Й 100%, 
Мп А,Е+ЪЙ 
где Х„[+ЬЙ — состояние узла а, в т-ом кортеже обучающей (контрольной) выборки; Х„[#+Й — 


модельные значения состояния узла а; М — объем обучающей (контрольной) выборки. 
Средняя относительная ошибка прогноза приведена на рис. 5. 
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Рис. 5. Средняя относительная ошибка прогноза для каждого кортежа контрольной выборки 


Средняя относительная ошибка идентификации составила 0,19%, средняя относительная ошибка 
прогноза — 0,23%, что свидетельствует об адекватности разработанной модели и позволяет рекомендовать ее 
для прогнозирования уровня подземных вод месторождения цементного сырья. 

Обсуждение и заключения. Разработана линейная динамическая окрестностная модель с 
переменными иерархическими окрестностями для прогнозирования уровня подземных вод месторождения 
цементного сырья, расположенного в задонско-елецком водоносном горизонте. 

Для параметрической идентификации и моделирования функционирования рассматриваемой 
динамической окрестностной модели разработано программное обеспечение на языке программирования С+-. 

Реализованная модель с достаточной точностью предсказывает уровень подземных вод в обследуемых 
скважинах и может быть эффективно использована для прогноза уровня подземных вод месторождения 
цементного сырья. 
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Выбор элементов комплекса защиты информационной системы предприятия на основе 
требований нормативно-правовых документов” 
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Е. А. УНепиго ', А. А. Геубоуа 7" 
+? Уоеортаа Зе Ошуегзиу, Уо1еоета, Визз1ап ЕедегаНоп 


Введение. Качество производственных процессов во многом 
зависит от инфраструктуры управления — в частности, от 
эффективности информационной системы (ИС). Менеджмент 
компаний уделяет все большее внимание обеспечению без- 
опасности этой сферы, на ее поддержку регулярно направля- 
ются финансовые, материальные и другие ресурсы. В пред- 
ставленной работе рассмотрены вопросы построения 
комплекса защиты информационной системы предприятия. 
Материалы и методы. Охрана ИС предприятия учитывает 
особенности объекта защиты и актуальные угрозы информа- 
ционной безопасности. В рамках данного исследования при- 
нято, что ИС представляет собой комплекс информационных 
ресурсов. По результатам специального анализа определены 
категории угроз информационной безопасности предприятия: 
взлом; утечка; искажение; утрата; блокирование; злоупотреб- 
ление. Выявлена связь данных угроз, компонентов ИС и эле- 
ментов комплекса защиты. Рассмотрены требования норма- 
тивно-правовых актов Российской Федерации И 
международных стандартов, регулирующих данную сферу. 
Показано, каким образом результаты данного анализа позво- 
ляют обосновать выбор элементов комплекса защиты ИС. 
Результаты исследования. Сравнительный анализ регламен- 
тирующей литературы, относящейся к данному вопросу, поз- 
волил выявить следующее. Разные документы предлагают 
разный набор элементов (подсистем) комплекса защиты ИС 
предприятия. Разрабатывая программу защиты ИС, следует 
руководствоваться Приказом ФСТЭК № 239 и стандартом 
800-82 Веу1510п 2 Се ю ГС$ Зесигйу. 

Обсуждение и заключения. Результаты представленного ис- 
следования являются основой для формирования программно- 
го комплекса интеллектуальной поддержки принятия решений 
при проектировании системы защиты информации на пред- 
приятии. В частности, можно разрабатывать гибкие комплек- 
сы, позволяющие расширять состав элементов (подсистем). 
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Введение. Информационные системы (ИС) все активнее используются в производственных и управленче- 
ских процессах. В связи с этим обостряется проблема информационной безопасности (ИБ) ИС. В частности, недо- 
статочная изолированность ИС упрощает несанкционированный доступ к ним [1, 2, 3]. Последствиями вредонос- 
ного воздействия на ИС могут быть простои производства, финансовые потери, а при реализации худшего 
сценария — даже техногенные катастрофы [4]. Таким образом, актуальной задачей является формирование ком- 
плекса защиты промышленных ИС, эффективно препятствующего злоумышленным действиям. 

Материалы и методы. Создание комплекса защиты информации основывается на результатах предпро- 
ектного обследования, в ходе которого определяются состав объекта защиты и актуальные для него угрозы. 

Объект защиты представляется как множество информационных ресурсов: 

ОБесЕ. ‚ = {МЕ,СС, 15,55, И/5, РЕ, О$,55, 45 „1Р,5и, КМ, 5М 14}. 

Здесь МЕ — множество сетевого оборудования; СС — множество каналов связи; /5 — множество инфраструк- 
турных серверов; 545 — множество систем хранения данных; И — множество рабочих станций пользователей; 
РЕ — множество периферийного оборудования; 05 — множество операционных систем; 55 — множество си- 
стемного программного обеспечения (ПО); 45 — множество прикладного ПО; /Р — множество информацион- 
ных процессов, протекающих в ИС предприятия; $и — подсети; Ю5М — множество съемных носителей инфор- 
мации; °М — электронные носители информации; [А — информационные активы. 

Множество актуальных угроз ИБ Тйгеа! определяется [5]: 

Тигесй = {Вгеаюте, Геак, Олюгпоп, [.055, Вост, АБизе}. 


Здесь Вгеате — угрозы взлома; Геак — угрозы утечки информации; Гу5ют’йо" — угрозы искажения; [055 — 
угрозы утраты; ВюсЯпие — угрозы блокирования информационных ресурсов ИС предприятия; АБизе — угрозы 


злоупотреблений. 
Комплекс, противодействующий данным угрозам, представляет собой систему защиты ИС предприятия 
(СЗИС) (рис. 1). 
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Рис. 1. Взаимосвязь объектов защиты и угроз в схеме СЗИС 





Система защиты информации 5Р/ (5ует оф ргыесйоп ор тюгтайоп) — двухуровневая и включает подси- 
стемы (компоненты) [6]: 
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— множество подсистемы (5и65узет)защиты информации; 
— множество средств защиты (МР, теап5 офргоесноп) информации. 
В общем виде структура СЗИС представлена на рис. 2. 
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Рис. 2. Обобщенная структура СЗИС предприятия 
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При определении компонентов комплекса защиты информации специалисты исходят из анализа имею- 
щейся нормативно-правовой документации и стандартов, действующих на предприятии. Необходимо также учи- 
тывать международный опыт. Довольно широко в мировой и отечественной практике применяется стандарт 800-82 
Веу1з1оп 2 Сливе 1ю шдизла! Сопго| Зузетаз (1С$) Зесигцу [1]. Он разработан Национальным институтом стандар- 
тов и технологий США. В нем, в частности содержатся рекомендации по повышению безопасности в системах 
промышленного контроля, включая системы диспетчерского управления и сбора данных. Показано, каким угрозам 
подвергаются организационные процессы и бизнес-функции, описаны типичные уязвимости. Особое внимание 
уделяется мерам безопасности и контрдействиям, которые следует предпринять в угрожающей ситуации. 

Результаты исследования. Отечественные нормативно-правовые акты (НПА), регламентирующие вопро- 
сы защиты ИС предприятия, условно можно разделить на две категории [6]: 

— НПА по обеспечению информационной безопасности автоматизированных систем управления технологическим 
процессом (АСУ ТП); 
— НПА по защите критической информационной инфраструктуры (КИИ). 

Следует особо отметить, что уязвимости в защите КИИ могут повлечь значительный материальный и эко- 
логический ущерб. Недостаточная охрана КИИ чревата социальными и военно-политическими проблемами. 

Проектирование системы защиты информации (в частности, при создании модели интеллектуальной под- 
держки принятия решений) предполагает предварительное проведение сравнительного анализа профильных нор- 
мативно-правовых актов Российской Федерации. Следует рассмотреть, например, следующие документы: 

— Приказ Федеральной службы по техническому и экспортному контролю (ФСТЭК России) от 14 марта 2014 г. 
№ 31 «Об утверждении требований к обеспечению защиты информации в автоматизированных системах управле- 
ния производственными и технологическими процессами на критически важных объектах, потенциально опасных 
объектах, а также объектах, представляющих повышенную опасность для жизни и здоровья людей и для окружа- 
ющей природной среды» [7]; 

— Приказ ФСТЭК России от 25 декабря 2017 г. № 239 «Об утверждении требований по обеспечению безопасности 
значимых объектов информационной инфраструктуры Российской Федерации» (проект) [8]; 

— международный стандарт 800-82 Веу1з1оп 2 Оли4е ю шиза! Сопго] Зузетв (1С$) Зесигиу [9]. 


Сравнительный анализ Приказов ФСТЭК России № 31 и № 239 представлен в табл. 1. 
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Таблица 1 
Элементы (подсистемы) комплекса защиты ИС предприятия в Приказах ФСТЭК России № 31 и № 239 





Подсистемы СЗИС 


Приказ ФСТЭК от 14.03.14 № 31 Приказ ФСТЭК от 25.12 2017 г. № 239 








Идентификация и аутентификация субъектов доступа и объектов доступа (ИАФ) 





Управление доступом субъектов доступа и объектов доступа (УПД) 





Ограничение программной среды (ОПС) 





Защита машинных носителей информации (ЗНИ) 











Регистрация событий безопасности (РСБ) Аудит безопасности (АУД) 
Антивирусная защита (АВЗ) 
Обнаружение вторжений (ОВ) Предотвращение вторжений (ПВ) 





Контроль (анализ) защищенности информации (АНЗ) |Защита информационной (автоматизированной) 
системы и ее компонентов (ЗИС) 








Обеспечение целостности (ОЦЛ) 





Обеспечение доступности (ОДТ) 





Планирование мероприятий по обеспечению защиты информации (ПЛН) 





Защита технических средств и систем (ЗТС) 





Обеспечение безопасности разработки программного | Реагирование на инциденты информационной без- 
обеспечения (ОБР) опасности (ИНЦ) 








Защита среды виртуализации (ЗСВ) Информирование и обучение персонала (ИПО) 





Управление обновлениями программного обеспечения (ОПО) 





Обеспечение действий в нештатных ситуациях (ДНС) 





Анализ угроз безопасности информации и рисков от| — 
их реализации (УБИ) 











Управление конфигурацией автоматизированной системы управления и ее системы защиты (УКФ) 








Примечание. Для большей наглядности отличия не только разнесены по разным ячейкам, но и выделены серым 
фоном. 





Итак, Приказ ФСТЭК № 239 предусматривает наличие в комплексе защиты ИС предприятия следующих 


подсистем: 

— аудит безопасности (АУД); 

— защита информационной (автоматизированной) системы и ее компонентов (3ЗИС); 
— реагирование на инциденты информационной безопасности (ИНЦ); 

— информирование и обучение персонала (ИПО). 


Следует отметить, что решение о комплектности СЗИС до известной степени зависит от финансовых во3з- 


можностей предприятия. Однако если стоимость защищаемых ресурсов и потенциальный ущерб от злоумышлен- 
ных действий выше, чем стоимость СЗИС, то целесообразно внедрить АУД и ЗИС. 


Сравнительный анализ Приказа ФСТЭК от 25 декабря 2017 г. № 239 и стандарта 800-82 Веу1з1оп 2 Сшае ю 


паизела| Сопёго| Зузепа$ (1С$) Зесигиу приведен в табл. 2. 


Таблица 2 


Элементы (подсистемы) комплекса защиты ИС предприятия в Приказе ФСТЭК от 25 декабря 2017 г. № 239 и 


стандарте 800-82 Веу151оп 2 Слмае ю шаиза! Сопёго] Зууетз (1С$) ЗесигИу) 








Подсистемы СЗИС 


Приказ ФСТЭК от 25.12 2017 г. № 239 800-82 Кеу1510п 2 Сливе ю 1С$ Зесигиу 
Идентификация и аутентификация (ИАФ) 
[епийЯсайоп ап аяфеписайоп 











Управление доступом (УПД) 
Зузет ап4 соттишсайоп$ рто{есйоп, 
ЗесигИу аззеззтеп{ ап4 аиогханоп 
Ограничение программной среды (ОПС) Зузеш ап4 шоппайоп ищеэтйу 











Защита машинных носителей информации (ЗНИ) 
Мела ргоесНоп 





Аудит безопасности (АУД) 
Ачат$ ап ассоита Пиу 











Антивирусная защита (АВЗ) Зузет ап шогтпайоп пиезтйу 
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Гйепвиге Е. А. ап4 ше ойег$. З@есйииз зареу расказе сотропеп5 ор ешегрие трогтайоп зуяет рюПоилиз гедшгетет$ 








Я а Приказ ФСТЭК от 25.12 2017 г. № 239 800-82 Веу1510оп 2 Слшае ю ГС$ Зесигиу 
Предотвращение вторжений (компьютерных атак) | Зузет ап4 шРоплайоп пиеэтиу 

(СОВ) 
Защита информационной (автоматизированной) | Зузет ап4 шРоплайоп ищеэтиу 
системы и ее компонентов (ЗИС) 














Обеспечение целостности (ОЦЛ) Зузет ап шогтайоп пиезтйу 





Обеспечение доступности (ОДТ) 
бузет апа зегу!сез асаи1иоп 





Планирование мероприятий по обеспечению безопасности (ПЛН) 
Р!аппипе, сопипзепсу р!аппте 





Защита технических средств и систем (ЗТС) 
Маиепапсе 





Реагирование на инциденты информационной безопасности (ИНЦ) 
шсеи гезропзе 





Информирование и обучение персонала (ИПО) 
Регзоппе! зесигИу 





Управление обновлениями программного обеспе- | Огоашхайоп — Ул4е шЮппайоп зесигйу ргоэтат 
чения (ОПО) тапазетепе сопёго]5 








Обеспечение действий в нештатных ситуациях (ДНС) 
РрБузса| ап епупоптеша! ргоесйоп 
Ау’агепезз ап4 тайпте 





Управление конфигурацией (УКФ) 
Сопйгигайоп тапазетег{ 





Ев В15К аззеззтепе 











— Зузет ап4 сотилишсайоп$ ргобесНоп 








Примечание. Для большей наглядности отличия не только разнесены по разным ячейкам, но и выделены серым 
фоном. 








В данном случае наиболее очевидны следующие различия: 

— 800-82 Веутз1юп 2 Слиде ю 1С$ Зесигиу объединяет в подсистеме Зу%ет ап шЮюгтайоп пиеэтиу функционал 
подсистем ОПС, АВЗ, СОВ, ОЦЛ, ЗИС, определенных в Приказе ФСТЭК; 

— 800-82 Кеу1з1юп 2 Сшае ю 1С5 Зесигиу предусматривает наличие подсистемы защиты систем связи — Зузет 
ап соттишсайоп$ ргоесНоп; 

— Приказ ФСТЭК объединяет в подсистеме управления доступом (УПД) функционалы подсистем Зуу%ет апа 
соттишсайоп$ ргоесйоп и Зесигйу аззеззтеп! ап4 аафоптайоп; 

— Приказ ФСТЭК объединяет в подсистеме планирования мероприятий по обеспечению безопасности (ПЛН) 
функционал подсистем Р1аппше и Сопйпзепсу р1апппз; 

— Приказ ФСТЭК объединяет в подсистеме обеспечения действий в нештатных ситуациях (ДНС) функционал 
подсистем РКузтса| ап4 епупоптета! ргофесйоп и Ау’агепез$ ап4 гашпие. 

Следует особо указать на подсистемы оценки рисков и защиты систем связи [10]. Это наиболее важные 
элементы комплекса защиты ИС предприятия, из тех, которые не предусмотрены отечественной нормативно- 
правовой документацией. Их внедрение позволит усилить защиту, оперативно реагировать на инциденты, возни- 
кающие в ИС предприятия, своевременно и точно противодействовать атакам. 

Выводы. Результаты анализа элементов комплекса защиты ИС будут использованы для построения моде- 
ли интеллектуальной поддержки принятия решений при проектировании СЗИС. В частности, планируется преду- 
смотреть возможность расширения состава подсистем СЗИС. Выбор элементов такого комплекса будет зависеть от 
оценки рисков, размера потенциального ущерба от вредоносного воздействия, стоимости компонентов СЗИС. 
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Введение. Полиномы нескольких переменных над полями 
Галуа лежат в основе теории кодов Рида-Маллера, а также 
используются в ряде криптографических задач. В работе 
изучаются свойства таких полиномов, заданных над 
произвольными полями Галуа нечетной мощности. Для 
полученных результатов предложены два практических 
приложения: схема разделения данных и декодер кодов 
Рида-Маллера. 

Материалы и методы. С использованием линейной 
алгебры, теории полей Галуа и общей теории полиномов 
нескольких переменных получены результаты, связанные 
с дифференцированием и интегрированием полиномов 
нескольких переменных над полями Галуа нечетной 
мощности. Для векторов построен и изучен аналог 
оператора дифференцирования. 

Результаты исследования. На основе полученных 
результатов о дифференцировании и интегрировании 
полиномов предложен новый декодер для кодов Рида- 
Маллера второго порядка и предложена схема 
организации разделенной передачи конфиденциальных 
данных, т.е. такой системы связи, в которой исходные 
данные на стороне отправителя разделяются на несколько 
частей и, независимо друг от друга, передаются по 
различным каналам связи, а на стороне получателя из 
принятых частей восстанавливаются исходные данные. 
Особенностью предлагаемой схемы является то, что она 
позволяет защищать данные, как от нелегитимного 
доступа, так и от непреднамеренных ошибок, при этом в 
обоих случаях используется один и тот же 
математический аппарат. Разработанный декодер для 
кодов Рида-Маллера второго порядка, заданных над 
произвольным нечетным полем Галуа, может иметь 
некоторое ограничение по числу исправляемых ошибок, 
однако, его использование целесообразно для ряда 
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Введение. Полиномы нескольких переменных над полями Галуа и их производные применяются в 
различных областях защиты информации. Некоторые вопросы, связанные с интегрированием и 
дифференцированием полиномов нескольких переменных, рассмотрены в ряде работ. Например, в [1] 
исследуются полиномы, заданные над простыми полями Галуа, в [2-4] получены результаты для булевых 
функций, ав [5—6] получены результаты для полиномов, заданных над троичными полями Галуа. 

В работе рассматриваются полиномы нескольких переменных, заданные над произвольными полями 
Галуа нечетной мощности. Для таких полиномов получены результаты, связанные с вычислением производных 
по направлению, а также с восстановлением полинома по набору его производных, вычисленных в базисных 
направлениях. Для полученных результатов предложены два возможных практических приложения: схема 
разделения данных и декодер кодов Рида-Маллера (РМ-коды). 

Предложенная схема разделения данных может быть использована для организации разделенной 
передачи конфиденциальных данных, т.е. такой системы связи, в которой исходные данные на стороне 
отправителя разделяются на несколько частей и, независимо друг от друга, передаются по различным каналам 
связи, а на стороне получателя из принятых частей восстанавливаются исходные данные. Особенностью 
предлагаемой схемы является то, что она позволяет защищать данные, как от нелегитимного доступа, так и от 
непреднамеренных ошибок. При этом в обоих случаях используется один и тот же математический аппарат, 
связанный с РМ-кодами и дифференцированием полиномов. Разделенная передача может быть использована 
как для повышения скорости связи, так и для обеспечения конфиденциальности данных за счет усложнения 
задачи перехвата из нескольких линий связи. Некоторые вопросы разделения данных рассмотрены в работах 
[7-11]. 

Для РМ-кодов второго порядка детерминированные декодеры известны только для некоторых 
значений мощности 4 полей Галуа. Например, довольно много известно декодеров для случая 4=2 , например 


[12-13], для случая 4 =3 и использования полунепрерывного канала связи сконструирован декодер [5]. В [14] 


предложен декодер кодов Рида-Маллера второго порядка, заданных над полями Галуа мощности 2, 4 и 8. 
Предлагаемый в данной работе декодер РМ-кодов второго порядка, заданных над произвольным нечетным 
полем Галуа, основан на редукции к кодам Рида-Маллера первого порядка, кодовые слова которых можно 
декодировать любым подходящим декодером. В случае РМ-кодов, заданных над полями мощности больше 3, 
предлагаемый декодер имеет некоторое ограничение по количеству исправляемых ошибок. Следует отметить, 
что использование предлагаемой схемы декодирования в случае полей мощности больше трех, несмотря на 
имеющееся ограничение, может быть целесообразным при невысоком уровне зашумления используемых 
каналов связи. 
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Дифференцирование полиномов нескольких переменных. Пусть 4=р’, р — простое нечетное 


число, 5= М, Е — поле Галуа мощности 4. Рассмотрим кольцо полиномов от т переменных Е/[х,...х„] над 
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ы (^) 
полем Е,. Линейное пространство полиномов из Е/[х,...х„| степени не выше г обозначим Е”[х,,...,х„]. 
Пусть Е’ — т-мерное линейное пространство над Е. 


ОИЗвОДНОЙ полинома Е [хь...х по направлению Е называется результат действия 
п Ех, ... БеЕ" 


т 


оператора дифференцирования [3]: 





(Бу) = Л@)-Л@), хеЕ’, (0 
где /;(х) = /(х +Ь). Легко показать, что О; Е и ], а оператор 
ВЕРЕ Е аьным, | (2) 
является линейным. 
Сумму координат вектора а“еЕЕ,, где р — простое, как натуральных чисел обозначим р(0.). 


Полиномы Г еЕ”[х,,...х„] будем записывать в каноническом виде 


Я) = У, ди? =ах + У пл У ам, (3) 


аеЕ,' р(@)=1 р(а)=2 


где при записи монома х“ = хр р показатели а; будем отождествлять с элементами поля | ‚› а слагаемые в 


каждой сумме будем располагать в лексикографическом порядке по возрастанию. Если последняя сумма в (3) 
(0 
равна нулю, то получаем полином из Е, [х...,х, ]. 


Лемма 1. Пусть 4=р’, р — простое нечетное число, Г(х)е Е”[х,...х,] — полином в 


каноническом виде (3), Б = (Ь,,...,6„)=ЕЕ,/. Тогда 


1 @) = Лю. + оби Лоо. ) +хАх" ы (4) 
(25) = Б (Лоо лю»: Ди), +2ХАБТ +БАБ" = 2ХАБ" + У(Б)— Хи, (5) 





где 
Лю. Ло. 12 Ла. 12 в Лоо 12 Ло. 12 
Ло. 12 Лоо. 00 Лон. /2 р Лою.ло 12 Лоо. о /2 
Лил. 12 Лон. 00 12 Лод. 00 яя Ло. л0 12 Лол. 12 


Лю. 1о 12 Лоо. 10 12 Лол. 10 12 в: Лию.20 ооо. 12 


Ло. 12 Лоо. 01 12 Лоо. 12 ие ооо. 12 Лоо. 02 


Доказательство. В случае простого поля Галуа доказательство содержится в [1]. Используя (1), (4) и 
симметричность матрицы А получаем: 


у (х) В С +5) = Лого + (х Ва }. + (х +5) А(х +5)" > 
(Р.Г) = 6) -1@)=Ь (К, Ли оз-- Ло} +Х4БТ +БАХ" +ЬАБ' , 


(Р.Г) = Б (Фо лв». +2АБТ +ЬАБ! = 2ЖАБТ + ГБ) Го 


Докажем теорему, которая определяет способ восстановления с точностью до постоянного слагаемого 
полинома из Е | по набору его производных, вычисленных в базисных направлениях. 





Рассмотрим полином ЛеЕ”’[х,,х,,...х„| вида (4): 


Уз = Лю. +х ры ) +хАХ" . 


Если 
(2, Хе) = ах, нах, +0] (6) 
. =... т 
то 
о. а в. м. 
А-1 В а . с ь | о 
2 0 В В В 
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—1 


Ло. 0 ок -ВАЬ! РА Ь, ыы Ь, 
Лю _ б СВАЬ . е ” у р, (8) 
Лоо. о, > Ь, АБ, " [2 7 т. 
Доказательство. ИЗз (5), (6) получаем: 
У 1=Ь..., т: АБТ = ее ЛВ, - Ло = о (9) 
Тогда 
ыы о 
ое ь |9 4... 9 
В. 6, 2 Вы)’ Ка 6, ют 


Следовательно, формула (7) верна. 


Из (4) следует, что для любого Бе В | 


(6) — 6.00 — м у +ЬАБ" . 


Возьмем в качестве р векторы Б ‚=В и воспользуемся равенством (В) Гри =@, из (9). Тогда 


У 1=Ь.., т: 6% 8 ‚АБР ао и > 


Следовательно, 
Гот ДАТ 1 1 1 
0% — БАБ Ь Ь, ры Ь, Ло. 
2_т ШТ 2 2 2 
0, —Б, АБ, , В = В Ха. 00 
т т тт т т т 
0% Ь, АБ, Ь Ь, Ь, Уи 


и формула (8) доказана. ® 
4-ичные коды Рида-Маллера АМ, (г,т). Рассмотрим РМ-коды над конечным полем Е, где а=р’, р 


— простое нечетное число, 5= М [15-16]. Элементы еж] являются информационными полиномами 
кода КМ. (т’,т); будем полагаль, что т>г>0, т>2. Вектор т ‚ составленный из коэффициентов 
информационного полинома Х(х,...,х„), называется информационным вектором. 

В векторном пространстве Е” зафиксируем некоторое упорядочение 

{04 , ..., 0, } (9%, = (д, 2,-... 0 и) п=а". (10) 
Произвольный информационный полином /(х) Ев м] кодируется путем вычисления его значений в 


т. 
точках упорядоченного пространства Е, : 


СОР) = (@и),...,Л(@,)), (р 


и тем самым определяется оператор кодирования 
С: в Е 
Коды Рида-Маллера определяются натуральными параметрами г и т (г<т) 
ВМ, (г, т) = {С (© ),...,Л(а,)) |1) ЕЕ [х,..„х„],Чев(Р) <} Е,, 


параметр и’ называется порядком кода. Они образуют семейство линейных [и,^,4]4-кодов, длина и и 





размерность К которых определяются по формулам 
„ [| 
т й / т-—1 
п= Ч > | = У уз (-10/ С. ны 2 


1=0 7=0 

где | | — округление до меньшего целого, а минимальное кодовое расстояние 4 кода КМ, (г,т) удобно 
2% Я 

вычислять, используя параметры дуального кода АМ, (7 ‚т), где ’ =т(а--г-1. Пусть р — остаток от 


деления 7+1 на 9-1: /'+1=0(а-Ю+р, где р<а-1, тогда параметр 4 кода КМ (г,т) задается 


выражением 


Деундяк В. М. и др. Дифференцирование полиномов нескольких переменных над полями Галуа 
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4=(р+1а°. (12) 
Отметим, что произвольный [п,К,4]«-код позволяет исправить {= | (а —1)/ 2] ошибок в одном кодовом слове 


[17]. 

Далее будем рассматривать РМ-коды порядков 1 и 2, заданные над полями Галуа нечетной мощности, 
соответствующие им информационные полиномы записывать в виде (3), а для нумерации координат 
информационного вектора использовать упорядочение (10). 


Лемма 2. Пусть ге {1;2}, 4>3, тогда минимальное кодовое расстояние кода КМ, (",т) вычисляется 
по формуле: 
4, =(а-г)а”', (13) 


а значения гарантировано исправляемых ошибок 1, = | (а —1)/ 2 кодами КМ,(г,т), ге\1;2}, связаны 
следующим образом: 
аа: (14) 

Доказательство. Воспользуемся тем, что ’<а4, и вычислим с и р — неполное частное и остаток от 

деления /` +1 на а—1 соответственно: 
о =(т(а-П-г)ау(а-=т-1; 
р=(т(а-Ю-г)то4(а-1 =а-1-г. 

Тогда из формулы (12) получаем (13). Из равенств 4 =(а-1а”", а, =(а-2)4"' получаем, что искомое 


неравенство (14) имеет вид 





1 | (9-Е «| (9-2) "1 
2 2 г 2 | 


Отметим, что 4, —Щ_ нечетное, а а, —_ четное, поэтому 


(а-2)а" -1|_(а-2”""-1 
2 . р 








2 








2 





(9-1а"'-1|_ (а-Ва""-2 
2 | О 


Следовательно, неравенство (14) приобретает вид 


1 (а-14”"' о р (а-2)а”"' —1 
2 2 в 2 





2 


легко видеть, что оно эквивалентно неравенству а>3.® 

Следствие. Если 4>3, то в (14) — строгое неравенство, то 4=3, то в (14) — равенство. 

В таблице 1 приведены параметры некоторых РМ-кодов. В трех верхних строках указаны такие 
параметры рассматриваемого кода АМ,(г,т) как 4, т, п. В следующих трех строках для кодов АМ, (1, т) 
содержатся значения: К, — размерность кода, 4, — минимальное кодовое расстояние и й — число 


исправляемых ошибок. В следующих трех строках представлены аналогичные значения А, 45, ь ДЛЯ КОДОВ 
ВМ, (2,т). 
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Таблица 1 


Значения параметров некоторых РМ-кодов 


2187 Вр= 3125 | 78125 СЕ: 2401 | 16807 


рр ев ве ава 
СИ О Е ВЕСИ И Е ВЕ И СО И ЕВ 
И И ЕЕ 


Теперь введем аналог оператора дифференцирования 2; ‚ действующего в пространстве полиномов (2), 





для пространства Е, ‚ где п=4”. Координаты векторов из Е, будем нумеровать векторами из упорядоченного 
множества Е" (см.(10)). Рассмотрим оператор сдвига т, : Е, —> Е, ‚ действующий по формуле 
т;(а)= (Я 5-Я +5), 

где а = (а ,...,аз, )= Е Ь = (6, ...,6„)Е Е’. Отметим, что оператор сдвига т; является перемешивающим 
биективным отображением. Линейный оператор Л;:Е, >Е,, являющийся аналогом Ш;, определим 
формулой: 

А;(а)=т,;(а)-а, а=(а, ,...а, )ЕЕ,. (15) 
Будем называть А‚(а) производным вектором вектора @ по направлению Ь. 


Лемма 3. Рассмотрим полином ГеЕ® [х,,х,,...х„|, вектор Ь = (В, ....Б,) Е”, операторы А;, О; и 


Хи 
С. Тогда 
(©) =С0), С.Л = АС. (16) 

Доказательство проводится прямыми выкладками и для 4=3 имеется в [6].® 

Отметим, что из (2) и (16) вытекает, что если С(у) е КМ, (2,т) , то А, (С(и)) е АМ, (т). 

Ниже рассмотрим примеры практических приложений полученных теоретических результатов. 

Схема разделения данных. Для разделения и восстановления данных в предлагаемой схеме 
используем [и,А,, 4], -код АМ, (1, т) и [п,К,,4.],-код КМ ,(2,т) ‚ заданные над произвольным полем Галуа Е, 
нечетной мощности. Значения 4 и т являются параметрами этой схемы. 

Алгоритм разделения данных 

Вход: информационный вектор %Е Е? кода АМ, (2,т) и упорядоченный набор базисных векторов 


В= {6, =(6,Ы,...6,) ЕЕ" и, (17) 


который является секретным ключом рассматриваемой схемы. 





Выход: векторы 5, Е =. 
Шаг 1. Сопоставим входному вектору % информационный полином = их) и закодируем его с 
Иа 
использованием (11) в вектор С(и’) е Е, кода КМ. (2,т). 


Шаг 2. Сформируем т производных векторов (см. (15)): 





А; (С(и)) = СО; (№) ЕЕ, , 1=Ьт, В ЕВ. 
Отметим, что С(Р; (у)) Е АМ. (т). 





В р://уезииК.Чопзеа.ги 


Шаг 3. Каждый вектор С(Б; (и)) ЕЕ, ,‚ {= т конкатенируем с коэффициентом Л» := "(0) кодового 


вектора С(\): 
5, =С(Ф; | Лоь Е". 


[9% 
> 


4 
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Далее векторы 5, е Е* ‚ =Ьт, передаются по т различным линиям связи. Очевидно, что во время 


передачи векторы 5,, 1=1,т , могут быть искажены. Таким образом, из канала связи будут получены векторы 


ю. 


1 


п-+1 ; 
Е =: 


00..00 


$," = (С(Ф, (м) Г. 
где (С(Р; (\)))" — возможно искаженный вектор С(Р, (и), а скаляр Л, ГИ — возможно искаженное значение 


Лю.ю. Скаляр /„ соответствующий 5,' обозначим / 


00.005 ° 
Алгоритм восстановления данных 
Вход: векторы 5,', {=1 т и секретный ключ В (см. (17)). 


Выход: вектор \'е Е . 


Шаг 1. Из каждого вектора 5, = тт ‚ выделяем две компоненты: (С(; ())))'Е Е и Л №: = тт : 


00..00,Г 
Шаг 2. Векторы (СБ; (и’)))" направляем в декодеры кода АМ, (1, т). Отметим, что декодеры могут 
быть использованы произвольные, например, [16], [18]. На выходе рассматриваемых декодеров формируются 
полиномы р; (ис 1 ры. ‚=1,....тТ. 


Шаг 3. Из коэффициентов / формируем вектор (Х’°, 


00..00.7 00..00,1 00..00,2 2 


т) 


00..00.т 


‚) и подаем на вход декодера 
кода КМ, (0,т), фактически совпадающего с кодом т-кратного повторения. Результатом работы этого 
декодера является скаляр Гуд. 

Шаг 4. Значения коэффициентов полиномов В; (СЕ? быв ‚ 1=Ьт, и ключ В (см. (17)) 
подставляем в формулы (7) и (8) и из полученных результатов строим полином /(х). Затем вычисляем 
искомый полином у”(х) = Г (%)+ Год. 

Шаг 5. Получателю сообщений выдаем информационный вектор и'е Е ‚ соответствующий полиному 
и"(х). 

Замечание 1. Корректность алгоритма восстановления данных зависит от числа ошибок, повредивших 
векторы 5, Е во время их передачи по линиям связи, а также от корректирующей способности 


используемых РМ-кодов первого порядка. Очевидно, что если используемые декодеры корректно восстановят 
векторы С(р; (»))и значение Гу, ТО и восстановление исходных данных с использованием результатов 


теоремы 1 будет корректным, и, следовательно, вектор и”’х), полученный на выходе алгоритма 
восстановления данных, совпадает с исходным информационным вектором \Е вы Отметим, что работа 
декодеров по восстановлению С(р; (у) корректна, если 

У=Ьт: 4, (ССР, (м), (СС, (и) < (а, -1/2 |, 
где 4,(х,У) — расстояние Хемминга между векторами х,у.Скаляр Г» восстанавливается корректно, если 


вектор О п и т Л, 


00..00 т 


), сформированный на шаге 3 алгоритма, содержит менее п/2 координат, 


отличных от значения Г», = и{0). Еслив 5, ЕЕ"", {=1, т, произошло ошибок больше, чем могут исправить 
используемые декодеры, то восстановление те Е” не гарантируется. 


Замечание 2. В предложенных алгоритмах разбиения и восстановления используются коды Рида- 
Маллера как первого, так и второго порядков, однако декодеры применяются только для кодов первого 
порядка. 

Замечание 3. Конфиденциальность передаваемых данных обеспечивается не только необходимостью 
знания ключа, но и за счет использования нескольких ЛИНИЙ связи, т.к. в этом случае перехват данных является 
более затруднительным для злоумышленника, чем нелегитимное получение данных из одной ЛИНИИ СВЯЗИ. 


Декодер кодов Рида-Маллера второго порядка. Сначала рассмотрим идею организации алгоритма 
декодирования, а затем опишем алгоритм по шагам. 


Зафиксируем некоторый базис В ={6, =(Б,№,...М)еЕ”}._„ в пространстве Е” 


ы „ где 4=р`, р — простое 


нечетное число. Предположим, что на вход декодера поступает У = С(и)-+е(е Е,), где » — информационный 


полином, С(\) — кодовый вектор [и,^,, 4, |, -кода КМ,(2,т), ее Е„ — вектор ошибок, для которого 
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и, (е)<ь, (18) 
где и/,(.)— вес Хемминга, #, = | (а, —10/ 2| . По вектору У с помощью оператора А; построим т 
производных векторов, вычисленных в базисных направлениях: 
А; (У) = А; (Сс) +е) = А; (СоХ)+А; (е), 1=1,..,т 
каждый из которых представляет собой вектор А; (С(»)) е КМ ,(1,т), искаженный вектором ошибок А; (е) 


Е Е" 


а? 
работающим до половины кодового расстояния (см., напр., [16], [18]), в случае, если число ошибок не 
превосходит 1 =| (4-1 /2 |, т.е. когда 


и может быть безошибочно декодирован произвольным декодером [и,А,,4,],-кода КМ. (т), 


ий, (А, (е)) <. (19) 
Если векторы А; (У) декодируются правильно, то искомый информационный полином и кода АМ (2,т) 


может быть восстановлен с использованием теоремы 1 с точностью до одной координаты, которая затем может 
быть найдена, например, декодированием по максимуму правдоподобия. Таким образом, для правильного 
декодирования У по предложенной схеме требуется выполнение условия (19). 

Алгоритм декодирования кода КМ, (2,т) 


Вход: параметры [и,К,,4, |, -кода КМ. (2,т), У =(У,,У,,,..У )ЕЁЕ,. 
Выход: восстановленный информационный вектор %. 


Шаг 1. Зафиксируем некоторый базис В = {5 =(&,В,...6.) ЕЕ пространства Е” и вычислим 


ыы т 
производные векторы по всем направлениям 5, ЕВ: 
А; (Г) =, (Р-Я. 
Шаг 2. Декодируем А; (7), 1=1,....т ‚ используя произвольный декодер АМ (т) -кодов, работающий 
до половины кодового расстояния, и в результате получаем векторы р; и их полиномиальные представления 
РС) = сих, + 2х, +... + х„ +0 С Е [х,,х,,..5 „|, = Ъ..ьт. 
Шаг 3. Используя полиномы Ру (х), 1=1.....т, по формулам (7) и (8) найдем полином Х(х) с 


нулевым свободным членом. 
Шаг 4. Для всех 2Е В. ВЫЧИСЛИМ 


2 


4 (2) = >С 18) +2). -У, 


где С(/<)+2), — @,-тая координата вектора С(/(х)+2) (см. (11)). Обозначим 2,значение 2 на котором 





функция (2) достигла своего минимума. 

Шаг 5. Результатом декодирования является вектор ® взаимно однозначно соответствующий 
информационному полиному их) = /(х) +2, .® 

Теорема 2. Для того, чтобы построенный алгоритм декодирования кода КМ, (2,т) исправил все 
ошибки в У=С(и)+е достаточно, чтобы выполнялись условия (18) и 

и, (е)<&/2, (20) 

где й = | (а, —0/ 2] ‚ 4, — минимальное кодовое расстояние кода КМ, (т). 

Доказательство. На шаге 2 на вход декодера поступают векторы 

А; (У) =т (7) -У =т, (С(®)) +1; (2)-С(®)-2 =А, (СС) +А, (©). 


Напомним, что А; (С(у))е КМ, (Ьт), а декодеры кода КМ, (1,т), работающие до половины кодового 





расстояния, исправляют до 1 ошибок в кодовом слове. Из условия (20) вытекает, что 
У, (А; (е)) < ит, (с; (е)) +ил, (е) =2мл,(е) <& 
Следовательно, векторы Р; ‚ которые формируются на шаге 2, совпадают с А; (С(С\)). Из этого вытекает, что 
на шаге 3, в силу теоремы 1, формируется полином Г(х)=и(х)-их0). 
Из условия (18) вытекает, что вычисленная на шаге 4 величина 2, равна свободному члену "х0) 
искомого информационного полинома их). Таким образом, на шаге 5 алгоритма получается искомый 


информационный вектор %. ® 
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Отметим, что для правильного декодирования по предложенной схеме требуется выполнение условия 
(20), из которого вытекает (19), хотя более естественным является условие (18). Рассмотрим связь между этими 
условиями. 


Лемма 4. Для кодов КМ, (2,т) в случае 4 =3 условия (18) и (20) равносильны, а в случае 4>3 при 


выполнении условия (18) выполнение условия (20) не гарантируется. 
Доказательство. Из следствия леммы 2 получаем, что при 4=3 справедливо равенство &/2=Ь,, т.е. 


правые части неравенств (18) и (20) совпадают, следовательно, из выполнения одного из них вытекает 
выполнение другого. При а4>3 из следствия леммы 2 получаем, что 1/2 <Ь, т.е. из выполнения (18) не 


вытекает выполнение (20). ® 
Замечание 1. В [5—6] описан декодер АМ, (2,т)-кода, где, как и в предложенном алгоритме, для 


зашумленного кодового слова строятся производные векторы, которые декодируются по алгоритму 
максимального правдоподобия, а затем по полученным значениям восстанавливается искомое информационное 
слово. Однако, в декодере из [5—6] строятся производные векторы во всех 3” возможных направлениях, а не 
только базисных, и используется иной механизм построения искомого информационного вектора. 


Замечание 2. Для кодов КМ й (2,т), а=3, предлагаемый декодер работает до половины кодового 
расстояния. Для кодов КМ,(2,т), а>3, предлагаемый декодер не гарантирует исправление всех ошибок, 


число которых не превосходит {[, ‚ но декодер будет работать корректно при соблюдении более слабого условия 


(20). Отметим, что использование предлагаемой схемы декодирования в случае полей мощности больше трех, 
несмотря на указанное ограничение, может быть целесообразно по следующим причинам. Во-первых, теория 
декодеров РМ-кодов второго порядка проработана недостаточно, но, если имеется декодер первого порядка, то 
предлагаемый декодер, являющийся надстройкой над ним, восполняет этот пробел. Во-вторых, при 
использовании каналов связи, вероятность ошибок в которых такова, что (20) выполняется, переход от РМ- 
кодов первого порядка к кодам второго порядка уменьшает избыточность (см. таблицу 1). 

Заключение. Получены теоретические результаты, связанные с восстановлением полиномов 
нескольких переменных над полями Галуа нечетной мощности по их производным. В качестве практических 
приложений полученных результатов предложены схема разделения данных и декодер РМ-кодов второго 
порядка. В дальнейшем представляется полезным исследовать процесс восстановления полинома по 
искаженным производным и разработать соответствующие модификации предложенных в настоящей работе 
практических приложений. 
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